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SECTION I 


INTRODUCTION 


FORTRAN is an automatic coding language especially suited to solving 
complex scientific and engineering problems. This capability is greatly enhanced 


by the use of subroutines from the Input/Output (I/0) and math libraries to 


perform the basic elements of the FORTRAN language - arithmetic, control, and 
. input/output. 


This manual describes the subroutines that make up the FORTRAN Library. 
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SECTION II 


INPUT/OUTPUT LIBRARY 


The FORTRAN Input/Output (I/0) subroutines perform the functions associated 
with the execution of and input/output requirements of the user's program. The 
descriptions in this section define the general mode of operation of the 
subroutines that constitute the FORTRAN I/O Library. In these descriptions, it 
is assumed that the reader is familiar with the manner in which the General 
Loader processes General Comprehensive Operating Supervisor (GCOS) and $ FFILE 
control cards. 


The I/O subprograms require that program execution be started with the 
FORTRAN initialization subprogram .FSETU, because one of its functions is the 
initialization of fault vector cell 25 (octal). This location contains the 
address where the "logical file/file control block" table begins. The user’ can 
also accomplish this initialization by calling SETLGT when he has created his 
own table. The library subroutines also depend on the limits of unused memory 
being expressed in fault vector cell 37 (octal). This is always done by the 
General Loader. When the $ OPTION FORTRAN card is encountered in the batch 
mode, the General Loader ensures that the FORTRAN initialization subprogram is 
loaded. In the time sharing mode, the RUN subsystem ensures that the proper 
initialization subprogram is loaded. 


LIBRARY CALLS 


A call to any I/O Library subroutine from a FORTRAN language program 
contains, as one of the arguments, the logical file code expressed as an 
integer. This integer is placed in character position 5 (bits 30-35) of cell 
-FBAD. (defined in subroutine FOPEN) by the called I/O subprogram. The called 
subroutine then calls FOPEN which searches the “logical file - file control 
block" table defined as follows: 

Ls Fault vector cell 25 (octal) is of the form 

ZERO TAB,O 
2% The actual logical unit (LGU) table has the form 
ZERO ENDTAB, 0 
TAB VFD 18/FCB1,6/LGU1,6/LGU2,6/LGU3 


VFD 18/FCB2,6/LGU4,6/LGU5 ,6/LGU6 


ENDTAB ZERO 0,0 
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wheres a. TAB - 1 contains the address of the last available location in 
the table. 


b. FCBl contains the address of cell LOCSYM of file control block 
ate 


Cs LGUL, LGU2, LGU3 are the FORTRAN logical files that reference 
file control block 1. Missing files are filled in with zeros. 


If more than three logical files reference the same file control block, the 
FCBl1 pointer and the additional files can occur at any other place in the 
table. 


There are as many entries in the table as needed to express the various 
file control blocks and logical files referencing them. After the last 
entry in the table, zeros fill out the table. 


FOPEN places the address of the file control block for the referenced file 
in bits 0-18 of cell .FBAD., but does not destroy the logical file code in bits 
30-35 (character position 5) of that cell. FOPEN then proceeds to open the file 
and return. 


The subroutine that called FOPEN now has the information necessary to 
perform calls to the proper File and Record Control subprograms. In the case of 
an output file, bits 30-35 of .FBAD. are used as the report code of the output 
record. Thus, if many logical files are connected to SYSOUT, they are separated 
automatically at printing time according to the calling code originally 
specified in the FORTRAN program calling sequence. 
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INPUT/OUTPUT LIBRARY SUBROUTINES 


gives the input/output function and the library 
function in the different modes -- batch and time sharing and using the BCD 
ASCII character sets. 


The input/output library subroutines are summarized in Table 2-1. Table 2-1 


and those that are user callable. 


Table 2-1. 


Function 


Linked Binary I/O Interface 

Short List Binary I/O Interface 

Short List I/O Processor 

Random Binary I/O 

Format Controlled Sequential 1/0 
(Common Procedure) 

Format Controlled Sequential I/O 
(BCD/ASCII Procedure and Data) 
Pormat Controlled Short List I/O 

Processor 
Namelist Input 
Namelist Input 

and Data) 
Namelist and Dump Output 

(Common Procedure) 
Namelist and Dump Output 

(BCD/ASCII Procedure and Data) 
Character String Assignment 
Character String Compare 
Output Stop and Pause Information 
Console Communications 
Object Time Debug Processor 
Memory Dump (Common Procedure) 
Memory Dump (BCD/ASCII Procedure 

and Data) 
Pre-execution Initializer 
Arithmetic Fault Processor 
Define Buffer(s) for File Control Block 
Define File Control Block 
Define Logical File Table 
Backspace Record 
File Closing 
Rewind and Endfile Processor 
End-of-File (On Input) Processor 
Initialization of End-of-File Processing 
Initialization of Data Error 

Processing 
File Opening 
Job Termination 
File Forward and Backspace 


(Common Procedure) 
(BCD/ASCII Procedure 


Carriage Control Simulator 


Call TSS Subsystem 

Create TSS Temporary File 

Specify Record Size, Random Binary File 
Set or Reset Some I/O Parameters 

of Run-Time Library 


User callable 


subroutine that performs’ the 
and 
Subroutines can be categorized as those implicitly called 


Input/Output Library Subroutines 
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SUBROUTINES IMPLICITLY CALLED 
Linked Binary Input/Output Interface 


FUNCTION 


Linked Binary I/O Interface consists of two subroutines (FRDB and _  FBLO) 
called for linked binary input/output interface; for processing linked 
binary elements; linked binary arrays; and for end-of-list processing for 
input/output statements of the following forms: 


READ (fc, END=S1, ERR=S2) list 


WRITE (fc, END=S1, ERR=S2) list 


fc = file code 

END = optional transfer location for end-of-file 

ERR = optional transfer for program error 

list = integer, real, complex, double precision, logical, and 


character elements and arrays 


The FBLO subroutine is called to process short list items; that is, arrays 
referenced without a subscript. 


The following sequence results for a READ statement where the list is an 
integer, real, logical, double precision, or complex element. 


TSX1 =F RDB. Initialization for 
TRA * + offset 

ZERO eBaLisey Esl READ 

ARG ' pointer to file code 

TRA ERR-clause 

TRA END-clause 

TSX1 -FBxT. List processing 

STy element 

TSX1 ~FRLR. Signals end of record 
x L for integer, real, complex, or logical elements 


D for double precision elements 
A for integer, real, complex, or logical elements 
AQ for double precision elements 


The following list processing sequence results for a READ statement where 
the list is a character element. 


EAA list 
TSX1 e-FBBC. (or .FBBCA in ASCII mode) 
NOP element size, DL 
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The following list processing sequence results for a READ statement where 
the list is an integer, real, complex, double precision, character, or 
logical array. 


TSxLl -FBxI. (.FXDI. for double precision arrays) 
ARG array locator 
ARG array size, DL (total number of words in array) 


L for integer, real, complex, and logical arrays 
S for character arrays 


The following sequence results for a WRITE statement where the list is an 
integer, real, complex, double precision, or logical element. 


TSX1 .FWRB. Initialization for 
TRA * + offset 

ZERO -E.L..,E.1. WRITE 

ARG pointer to file code 

TRA ERR-clause 

TRA END-clause 

LDx element List processing 

TSX1 .FByT. 

TSX1 ~FWLR. Signals end of output 


x = A for integer, real, logical, or complex elements 
= AQ for double precision elements 

y = L for integer, real, complex, or logical elements 
= D for double precision elements 


The following list processing sequence results for a WRITE statement where 
the list is a character element. 


EAA list 
TSX1 «-FBBC. (or .FBBCA in ASCII mode) 
NOP element size, DL 


The following list processing sequence results for a WRITE statement where 
the list is an integer, real, complex, double precision, character, or 
logical arrays. 


TSX1l -FBxO. (.FXDO. for double precision arrays) 
TRA * + offset 

ZERO -E.L..,5.1. 

ARG array locator 

ARG array size, DL (total number of words in array) 
x = L for integer, real, complex, or logical arrays 


S for character arrays 


The following table, Table 2-2, contains a summary of the linked binary 
subroutines and entry points. 


t 
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Table 2-2. Linked Binary Subroutines and Entry Points 


Type Short List Short List 
Initialization Element Input Outout End of Input Fnd of Output 


Integer . . A FRDB 
Real , . 7 PRL FRDB 
Comp lex : : : : : RL. FRDB 
Double-Precision : * FRDB 
Logical : 3 | : : FRDB 
Character PRL FRDB 


Subroutine FBLO calls FSLI to set up indexing for the processing of the 
input/output of nonsubscripted arrays. 


CALLING SEQUENCE 


See Function. 


METHOD 


Binary I/O is standard system format I/O unless otherwise specified and 
reads or writes media code 1 records. The binary I/O routine uses File and 
Record Control buffers, eliminating the need for an internal work area. 
Records are assumed to be variable length records and can be up to 4095 
words in length. Logical, variable length records too large to fit in one 
buffer are partitioned and follow the rules for partitioned records as 
defined below. 


Ls Only logical, variable length (VLR) records too large to fit in 
one buffer are partitioned. 


2. A partitioned record begins in a new buffer and ends ina buffer 
with no other records in that final buffer. 


ce Partition indicators: 


Initial Block: 


RCW bits 24-25 = 01 
Intermediate Blocks: 
RCW bits 24-25 = 10 


Bits 26-35 


Partition counter starting with 
2 and advancing by 1 for each 
new block in the record. 
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Final Block: 


RCW bits 24-25 


i 


ll 


Bits 26-35 


Partition counter described above. 


Valid batch mode loader $ FFILE options include NSTDLB, NBUFFS, BUFSIZ, 
FIXLNG, and NOSRLS. Additional information on FRDB can be obtained from 
documentation at the front of the assembly listing. 


RETURN 


Normal return is to the next executable statement in the calling program. 
Error conditions are handled as follows: 


Ls Error code 39 when an end-of-file condition, other than 17 or 23 
(octal numbers ), is detected. 

2% Error code 33 because Read after Write is illegal. 

ce Error code 40 when list exceeds logical record length. All 


remaining list items are set to zero. 


4. Error code 41 when SYSOUT or fixed length records are not smaller 
than block size. 
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Short List Input/Output Processor 


FUNCTION 


FSLI provides list processing for a nonsubscripted array in conjunction 
| with subprograms FSLO, FRDB, and FRRD. 


CALLING SEQUENCE 


The subprograms FSLO, FRDB, and FRRD have calling sequences’ to the 
subprogram FSLI to accomplish the data transfer. 


For Formatted I/O (FSLO) 


Output Input 
E STX2 -FSLIT E STX2 -FSLITI 
TSX2 - FSLI (m) TSX2 . FSLI (m) 
E LD(r) ** 2 E fTSxX1l c 
TSX1 Cc ST (r) **,2 
E ZERO p,0 E ZERO p,l 
™m O,1 


Kh 


A for single precision 

AQ for double precision 

-FCNV. for single and double precision 
-FCNVI for integer 

~FCNVR for real 

-FCNVD for double precision 

~-FCNVL for logical 

~-FCNVC for character 

1 for single precision 

2 for double precision 


For Binary I/O (FRDB, FRRD) 


LDA n,DU 

(ADLA 1,DL) double precision only 
XED . FSLIB, AU* 

ARG ENDTL 

ARG CRTAL 

ARG SVX1 
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METHOD 


For formatted I/O, FSLI initializes the I/O loop and processes the entire 
array, starting with the first element. For binary I/0, an EIS move is 
used to process the entire array. 


RETURN 


For formatted I/0, FSLI restores C(X2) = C(.FSLII)0-17 and exits to FSLO. 
For binary I/O, return is directly to the object code. 
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Random Binary Input/Output 


$ 
£ 


FUNCTION 


FRRD is called for random binary source statements of the following form: 


READ (f'n,optl,opt2) list 


WRITE(f£'n,opt2) list 


f = a file reference 

n = sequence number of the logical record 

optl = end of file transfer (optional) 

opt2 = error transfer (optional) 

list = input/output list that can consist of single precision, 


double precision, character, or array items 


This subroutine also includes the RANSIZ (Random Record Size) subroutine 
that is called by the following sequence: 


CALL RANSIZ (argl,arg2,arg3) 


argl = a file code 
arg2 = the logical record size 
arg3 = standard system format 


When arg3 is missing or zero, the file is assumed to be in standard system 
format. 


RANSIZ needs to be called only once for each file. 


CALLING SEQUENCES 


FRRD contains different entry points for I/O list and end-of-list 


processing. The sequences generated are similar to those generated for 
linked binary. 


TSX1 -FRRD. For READ (n',f) list 

TSX1 ~FRWD. For WRITE(n',f) list 

TSX1 -FRLT. For single precision I/O list entry 
TSX1 ~FRDT. For double precision I/O list entry 
TSX1 -FRST. For binary character (BCD) 

TSX1 -FRSTA. For binary character (ASCII) 

TSX1 -FRLI. -For single precision short list input 
TSX1 eERDE . For double precision short list input 
TSX1 -FRSI. For binary character short list input 
TSX1 ~FRLO. For single precision short list output 
TSX1 ~FRDO. For double precision short list output 
TSX1 -FRSO. For binary character short list output 
TSX1 ~FRRR. For binary input end-of-list 

TSX1 ~FWRR. For binary output end-of-list 

CALL RANSIZ Random Record Size Routine 
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METHOD 


Random binary files are standard system format (unless specified 
otherwise), 320-word blocks, with fixed record size. The module calculates 
logical record location based on physical device block size (64 words), 
logical record size (given by CALL RANSIZ command), buffer block size 
(320), and record number. The record size is compared to the buffer block 
Size and if greater the file is considered partitioned and follows’ the 
rules for partitioned files of standard system format (see the File and 
Record Control manual). 


The user should express record size as the sum of items in the I/O list in 
computer words, remembering that storage allocation for character variables 
and arrays differs when compiling BCD or ASCII (six characters versus four 
characters per word). The record size specified via RANSIZ should not 
include the record control word for the logical record in the size 
argument. An abnormal termination with an associated error message occurs 
when records are read that fall outside the file limits. 


RANSIZ needs to be called once prior to the first file reference. FRRD 
calls FSLI for the processing of the input/output of nonsubscripted arrays. 


RETURN 


3/77 


Normal return is to the next executable statement of the calling program, 
except as specified below: 
nie Error code 72 if the list exceeds logical record length. 
2. Error code 73 for any of the following conditions: 
a. Zero block count on random read. 


b. Block serial number error on random read. 


Cx Record count of zero on random read. 
Be Error code 74 if the file is not present. 
4. Error code 75 for a zero or negative record number on a= random 


read or write. 


5s Error code 76 where the record size is not given. Random files 
have fixed length records. Record size must be specified by a 
CALL RANSIZ for the referenced file. This CALL should precede 
the first I/O to a linked file. 


6. Error code 77 for an attempted random I/O to a linked file. 


7. Error code 78 when the record number in the random read or write 
statement is outside file limits. However, if "END=" is given 
(for READ only), that exit is taken if the record number is 
outside the file limits. A write outside file limits causes the 
file space to grow. 


Ss Error code 79 when the list is greater than specified record 
size. 


9. Error code 80 when file space is exhausted (on output, the file 
can no longer be "grown"). 


2-11 DD20A 


Format Controlled Sequential Input/Output 


FUNCTION 
Formatted sequential input/output is contained in five subroutines: 


FDIO contains the procedures that are common to BCD and ASCII 
formatted sequential input/output. 
FRDD contains the BCD peculiar procedures and data. 
FRDA contains the ASCII peculiar procedures and data. 
FSLO contains the short list input/output interface. 
_FSLI contains the short list input/output processor. 


Formatted sequential input/output is accomplished by calls to these 


subroutines to initialize the input/output; to process the input/output 
list;-and for end-of-file processing. 


Subroutines FRDD and FRDA also contain the processing for the DECODE and 
ENCODE statements. 


CALLING SEQUENCE 


The calling sequences to initialize the input/output are as follows: 


TSX1 entry point in FRDD or FRDA 

TRA * + offset 

ZERO sBeliss, Eels 

ARG file code for formatted statements or 41,42,43 for 
READ, PRINT/WRITE, or PUNCH respectively. 

ARG format statement reference or FORMAT(V) for list 

directed input/output. 
TRA ERR-clause 
TRA END-clause 
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The subroutines and entry points for initialization of input/output are as 


listed in Table 2-3. 


Table 2-3. Initialization of I/0 


Entry Entry 
Module { Point Module; Point 


WRITE (fc,format) list 
WRITE format, list 
PRINT format, list 
PRINT, list 

PUNCH format, list 
READ (fc,format) list 
READ format, list 
READ, list 


The calling sequence for list processing where list is integer, 
double precision, logical, or complex element is as follows: 


LDx element if output 
TSX1 -FCNVy 
STx element if input 


x = A for single precision items 

= AQ for double precision items 

= I for integer, R for real and complex, 
D for double precision, and 
L for logical 


The calling sequence for list processing of character elements 
follows: 


EAA element 
TSX1l ~FCNVC 
NOP element size, DL 


The calling sequence for all short list routines (except character) 
follows: 


TSX1 short list routine 
ARG array locator 
_ ARG array size 
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real, 
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The calling sequence for character short list input/output is as follows: 


TSX1 ~FSCX. 

TRA * + offset 

ZERO elolicigg Bs Ex 

ARG array locator 
ARG array size, DL 
NOP element size, DL 
x = O for output 


I for input 


The following table, Table 2-4, gives the subroutines and entry points for 
list processing and end-of-file processing. 


Table 2-4. List and End-of-File Processing 


Type Short List Short) List End-of -File End-of-File 
Element Out put Input Out put Input 
Ent ry Entry Entry Entry Entry 
Module Point Modu le Point Module Point Module Point Module Point 


BCD 
Integer FCNVI FSLO FS1O FSLO FSII FRDD FFIL FRDD FRIN 
Real FCNVR FSLO FSRO FSLO FSRI FRDID FEIL FRDD FRTN 
Comp lex FCNVR FSLO FSRO FSLO FSRI FRDD FRIL FRDD PRTN 
Double-Precision FCNVD FSLO FS20 FSLO FS21 FRDD FFIL FRDD FRTN 
Logical FCNVL FSLO FS BO FSLO FSBI FRDD FFIL FRDD FRTN 
Character FCNVC FSLO FSCO FSLO FSCI FRDD FFIL FRDD FRTN 

ASCII 
Integer FCNVI FSLO FSIO FSLO FSIL FRDA FPILA FRDA FRTNA 
Real FCNVR FSLO FSRO FSLO FSRI FRDA FFILA FRDA FRTNA 
Comp lex FCNVR FSLO FSRO FSLO FSRI FRDA FFILA FRDA FRTNA 
Double-Precision FCNVD FSLO FS20 FSLO FS2T FRDA FFILA FRDA FRTNA 
Logical FCNVL FSLO FSBO FSLO FSBI FRDA FFILA FRDA FRTNA 


Character FCNVC FSLO FSCO FSLO FSCI FRDA FFILA FRDA FRTWA 


The calling sequences for the DECODE and ENCODE statements are: 


TSX1 . «FEF DEC for DECODE 

TSX1 eFENC for ENCODE 

TRA * +4 or 6 

ZERO -E.L..,E.I1. 

ARG buffer 

ARG format 

TRA ERR~-clause 

TRA END-clause 

NOP record size in characters, DL 
NOP buffer size in words, DL 
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The following calling sequences are used for internal conversion routines 
with format: 


CALL »~BDCNV (buf,form,wda,lines) memory to buffer 


CALL ~DBCNV (buf,form,wda,lines) buffer to memory 


lines = number of lines specified by FORMAT statement. 
wda = word count of each line. 

form = starting location of format. 

buf = starting location of buffer. 


METHOD 


On the entry for initialization, the file code is determined and the file 
is opened and read if input. An appropriate tally word for data transfer is 
built. For output, the address of a 4l-word buffer is used in the tally 
word. (This provides for a maximum length record of 160 characters). For 
input, the address of the File and Record Control buffer is used directly 
so no intermediate storage is necessary; this allows records up to 4095 
characters in length to be processed. After initialization, the format scan 
begins and when a field separator or a right parenthesis is encountered, 
return is made to the user's generated code. The generated code causes 
control to transfer to the appropriate entry point for list processing. 
Thus, each list item corresponds to a format specification. This continues 
until all list items are processed, rescanning the format if necessary. The 
end-of-list entry determines that processing of the record is complete. 


The processing for ENCODE and DECODE is the same except that the address of 


the in-memory buffer provided by the user is used in the data transfer 
tally word. 


The routine FDIO contains code that is common for ASCII/BCD character sets 
while FRDA and FRDD contain ASCII and BCD unique code and symbols. 


FSLI is called by FSLO to set up indexing for the processing of short list 
input/output. . 


RETURN 


Normal return is to the next executable statement of the calling program, 
except as specified for the following error conditions: 


Le Error code 31 for an illegal format statement. Format scan 
proceeds as for end of format. 


2% Error code 32 for an illegal character in data or bad format. 
Data scan treats illegal character as zero. 


3s Error code 33 if user attempts to read an output file. Execution 
is terminated, 


4, Error code 34 for illegal end-of-file mark. FRWD calls FEOF for 
error recovery. 
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3/77 


: 
t 


‘Error code 57 for illegal character for L conversion. Data _ scan 


treats illegal character as space. 
Error code 81 for excessive line length. 


Error code 82 for illegal character as first nonblank character. 
Execution is terminated. 


Error code 86 for a recursive entry to the I/O routine. 
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Namelist Input 


FUNCTION 


The NAMELIST input function is contained in the following subroutines: 
FNLI contains the procedures that are common to BCD and ASCII NAMELIST 
input. 
FVFI contains the BCD peculiar procedures and data. 
FIFA contains the ASCII peculiar procedures and data. 
NAMELIST input is accomplished using File and Record Control variable 
length records for source statements of the following forms: 
NAMELIST /naml/v,w,x 


READ (fc,naml,END=S1, ERR=S2) 


fc = logical file reference. 

naml = NAMELIST name. 

Sl = (optional) is the statement label pointing to the statement to 
which control is transferred upon encountering an EOF for a 
READ. 

S2 = (optional) is the statement label pointing to the statement to 


‘which control is transferred upon encountering an error. 
v,W,X = NAMELIST variables. 


CALLING SEQUENCE 


TSX1 -FVFI. (BCD) 

TSx. -FVFIA. (ASCIT) 

TRA * + offset 

ZERO she bie peel. 

ARG pointer to location of DEC FC 

ARG pointer to location of BCI 1,NAM1 

ARG pointer to location containing address of the EOF exit 
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METHOD 


The input file is scanned for the proper NAMELIST name. When the name is 
found, it is verified that the variables are included in the NAMELIST and 
the input values are stored according to the type specified in the NAMELIST 
table. The entries for the above NAMELIST example (NAM1) are: 


BCI 1,NAM1 

BCI LV 

TALLYD LV;.C, Kk 

ZERO N1,N2 Present only when K = 6 (character data) 
ZERO PROD, di. if c=1 

ZERO 0/-1/d, *d*d3 ,d} *d2 if C = 2 

ZERO -1/d) *d2...ds5.,d1 *d2...*d4 if Cc = 3 

ZERO #1, di. 4do wes. Gg if c= 4 

BCI 1,W 

BCI 1X 


NAM] is the NAMELIST name 
LV is the location of the variable 
V is the variable 


Cc 07777 if V is not dimensioned 

if V has 1 or 2 dimensions 

if V has 3 or 4 dimensions 

if V has 5 or 6 dimensions 
7 


if V has dimensions 


| 
mW Ne 


K is the data type as follows: 
for integer 

for real 

for double precision 

for complex 

for logical 

for character 


x 
fou ie 
DU WN FE 


PROD is the product of all dimensions 
N1,N2 are present when K=6 (character data) 

Nl = number of words per element 

N2 = number of characters per element (maximum of 120) 
ad, ,dz, ...,dm are the dimensions for the variable 
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Bits of index register 6 are set (=l) to indicate the following: 


Bit number Meaning 
0 Logical variable 
1 Integer variable 
2-4 Free 
5 Slash preceded current 


variable name 

6 Complex mode not allowed 

7 Complex mode 

8 Neither star nor character 
field allowed 


9 Decimal point present 
10 Negative number 
11 Negative decimal exponent 
12 Decimal exponent allowed 
13 D field allowed 
14 E field allowed 
15 No decimal point allowed 
16 Two words per element 
17 One word per element 


EXAMPLE 


The following example shows the expansion using variable V with 0,1,2,3,4 
dimensions (these are designated as V0,V1,V2,V3,V4). The dimensions (dj) 
are designated as L,M,N,O. 


VO = undimensioned 


v1 (L) 
V2 (M,N) 
V3 (0,P,Q) 


v4 (R,S,T,U) 


L through U are integers 


Ls With V not dimensioned 
BCI - 1,V0 
TALLYD V0,07777,K 

2. V with a dimension of 1 
BCI 1,Vl 


TALLYD V1,1,K 
ZERO L,L 


3. V with a dimension of 2 (V2) 


BCI 1,V2 
TALLYD V2,1,K 
ZERO M*N,L 


2e19 DD20 


4. V with a dimension of 3 (V3) 


BCI 1,V3 

TALLYD V3,2,K 

ZERO O*P*Q,0 
ZERO 0 or -1,0*P 


5s V with a dimension of 4 (V4) 
BCI 1,V4 
TALLYD V4,2,K 


ZERO R*S*T*U,R 
ZERO R*¥S*T,R*S 


6. Where K = 6 (character data), there is a slight variation in the 
expansion for ASCII and BCD character data. For example, 


CH*10 will cause the following expansion for BCD: 


BCI 1,CH 
TALLYD CH,07777,6 
ZERO 2,10 


CH*10 will cause the following expansion for ASCII: 


BCI 1,CH 
TALLYD CH,0O7777,6 
ZERO 3,10 


RETURN 


Normal return is to the next executable statement of the calling program. 
Under the error conditions listed below, execution is continued only if the 
user initialized FLGERR, which causes a normal return with an indication 
that bad data was encountered, or if the ERR = option was used. FLGERR and 
ERR = can be used in any combination. 


Ls Error code 42 for illegal heading card. FNLI continues as_ for 
end-of-data. 


as Error code 43 for illegal variable name. FNLI continues as_ for 
end-of~-data. 


3. Error code 44 for illegal subscript or when the array size is 
exceeded. FNLI continues as for end-of-data. 


4. Error code 45 for illegal character after right parenthesis. Data 
scan assumes a comma between right parenthesis and the next 
character. 

Ds Error code 46 for illegal character in the data. Data scan treats 


the illegal character as zero. 


6. Error code 48 for illegal logical constant. Data scan treats an 
illegal constant as .FALSE. 


vee Error code 52 for an illegal character field. FNLI continues as 
for end-of-data. 


8. Error code 59 for an empty character field. 
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Namelist and Dump Output 


FUNCTION 


Three subroutines process NAMELIST, DEBUG, DUMP, and PDUMP output using 


standard File 


and Record Control variable length records for source 
statements of the following forms: 


NAMELIST /naml/v,w,x 


WRITE (fc,naml,ERR=S2) 


or 
CALL DUMP (arguments) BCD 
CALL DUMPA (arguments) ASCII 
CALL PDUMP (arguments) BCD 
CALL PDUMPA (arguments) ASCII 


DEBUG is invoked via General Loader control cards. 


The subroutines that process NAMELIST and DUMP output are as follows: 


FNLO contains the procedures that are common to BCD and ASCII NAMELIST 
output. 


FVFO contains the BCD peculiar procedures and data. 


FOFA contains the ASCII peculiar procedures and data. 


fc 
S2 


arguments 


V,W,X 


CALLING SEQUENCE 


TSX1 
TSX1 
TRA 
ZERO 
ARG 
ARG 
ARG 
ARG 


-FVDO 
-FVFO 


Hy 


* + offset 


eBehsc, 
pointer 
pointer 
pointer 
pointer 


E. 
to 
to 
to 
to 


logical file reference 

(optional) is the statement label pointing to the 
statement to © which control is transferred upon 
encountering an error 

see "Memory Dump" for the memory dump subroutines, in this 
section 

NAMELIST variables 


(For DEBUG, DUMP, or PDUMP output) 
(For NAMELIST output) 


I. 

location of DEC FC 

location of BCI 1,NAM1 

location containing address of error exit 
location containing address of EOF exit 
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METHOD 


These subroutines scan a NAMELIST table and print the current value of each 
NAMELIST variable in the format specified by its entry in the table. Refer 
to the description of NAMELIST input for a SEEcE ARE LON of the entries for 
NAMELIST, DUMP (DUMPA), and PDUMP (PDUMPA). 


For NAMELIST, DUMP, and PDUMP output or for NAMELIST input, the entry for a 
NAMELIST veriable, V, is of the following form: 


BCI 1,V always present 

TALLYD LV,D,K always present 

ZERO PROD,D1. present if D is 1 or 2 t 
ZERO O,D1*D2 present if D is 2 


LV is the location of variable V (must not be zero) 
D is octal 7777 if V has no dimensions 
1 if V has 1 or 2 dimensions 
2 if Vv has 3 dimensions 
Dl and D2 are the first and second dimensions (if any) 
PROD is the product of the dimensions 
K for octal 
for integer 
for real 
for double precision 
for complex 
for logical 
for character 


DUR WHF O 


For DEBUG output, the entry for variable V is: 


BCI 1,V always present 

TALLYD LV,D,K+16 always present 

ZERO r2,12 present if D is not octal 7777 

Il and 12 are the initial and final subscripts 


respectively, in increments of D (or 
1 if D is zero) 


If V is octal 777777777777, the output is a memory dump from Il to 12 
in increments of D elements (1 or 2 words). For memory dumps, LV 
is ignored and can be octal 7777. 

K is the same as for NAMELIST 

If V is not octal 777777777777, LV and K are the same as for NAMELIST 
and the octal value 7777 for D indicates that LV has. no 
dimensions. 


All arrays in the DEBUG mode are considered as one-dimensional arrays. 


RETURN 


Normal return is to the next executable statement in the calling program. 


2-22 DD20 


Character String Assignment 


FUNCTION 


FCHA moves character strings from one character variable to another when 
the character variable is used in an assignment statement. 


CALLING SEQUENCE 


The following calling sequence is for A = B where A and B are character 


variables. 
TSX1 ~FCHM. (or .FCHMA for ASCITI) 
TRA * + 6 
ZERO eBiticcc. eke 
ARG B,I 
ARG J,I 
ARG A,I 
ARG eet 
B = location of sending field 
J = location of variable containing size (in characters) of B 
A = location of receiving field 
I = location of variable containing size (in characters) of A 


METHOD 


The character string identified by character variable B is assigned to 
character variable A. If size B>size A, B is truncated at the size of A. 


If size B< size A, the remaining characters in A after the assignment are 
blank filled. 


RETURN 


Return is to the next executable statement in the calling program. 


2-23 DD20 


Character String Compare 


FUNCTION 


FCOM compares BCD or ASCII character strings for character variables 


appearing in logical or relational expressions of the form IF (A .EQ. B) 
where A and B are character variables. 

CALLING SEQUENCE 
TSX1 ~-FCOM (or .FCOMA for ASCIT) 
TRA * + 6 
ZERO -E.L..,E.I1. 
ARG A,I 
ARG a a 
ARG B,I 
ARG J,1 
A = address of character string A 
I = address containing size of A in characters 
B = address of character string B 
J = address containing size of B in characters 

METHOD 
Le The length of both strings in words is determined. If the lengths are 

different, the difference is saved. 
2% The two strings are compared on a word-for-word basis over equal 
length. Indicator register is set as follows: 
a. If ACB, the carry indicator is set off 
(At the first non-compare, the comparison stops; > or < is 


determined by the collating sequence of the character set.) 


b. If A=B, the zero indicator is set 
Cc. If A>B, the carry indicator is set 
3. If the strings are equal over their equal length but one 


string 


longer than the other, the difference is compared against blanks. 


the difference is blanks, the strings are equal. 


RETURN 


Return is to the next executable statement in the calling program. 


is 
If 
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Output Stop and Pause Information 


FUNCTION 


FPAW outputs a message to the system console (user terminal for time 
sharing) when a PAUSE or STOP statement is encountered in the _ source 
program. 


CALLING SEQUENCES 


-FPAW. (or .FPAWA for ASCII) 
*+9 
-E.L..,E.I. 


Generated for a PAUSE statement with no message. 


-FPAI. (or .FPAIA for ASCII) 
¥43 

-E.L.. rol. 

pointer to integer 


Generated for a PAUSE statement with integer message. 


~-FPAC. (or .FPACA for ASCIT) 
+4 

~E.L.+,E.L. 

pointer to character message 
number of characters, DL 


Generated for a PAUSE statement with character string message. 


~FIXT. (or .FIXTA for ASCII) 
#43 

Beles ey Bs Ls 

pointer to integer 


Generated for a STOP statement with integer message. 


~-FCXT. (or .FCXTA for ASCII) 
*+4 

sBeliseehe Ls 

pointer to character message 
number of characters, DL 


Generated for a STOP statement with a character string message. 
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METHOD | 


This routine returns the following information to the user terminal (time 
sharing user) or the system console (batch user) when the following 
statements are encountered: 


STATEMENT | | MESSAGE 
PAUSE PAUSE LINE # nnnn 

SNUMB xxxx-xx (batch only) 
PAUSE integer PAUSE nnnnnn 

SNUMB xxxxx-xx (batch only) 
PAUSE character string PAUSE user message 
STOP integer | STOP AT LINE nnnn 
STOP character string STOP user message 


RETURN 


Normal return is to the next executable statement in the calling program. 
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Object Time Debug Processor 


FUNCTION 


FDBG decides whether or not to produce DEBUG output based on the contents 
of the IF and FOR statements. 


CALLING SEQUENCE 


FDBG is entered by a DRL (Derail) instruction. 


ME'THOD 


The location of the DRL instruction by which FDBG was entered is checked to 
see if it was inserted by the General Loader (a legal DEBUG request), or if 
it was originally present in the interrupted program. If it was a DEBUG 
request, the IF and FOR statements are examined in the order in which they 
were specified. If these statements are satisfied, FDBG writes information 
describing the particular input it is interrupting. It then calls FVFO or 
FOFA for the list output, using the NAMELIST table supplied in the DEBUG 
table by the General Loader. The instruction replaced by DRL is’ then 
executed and control is returned to the next instruction in the interrupted 
program. If DRL was not inserted by the General Loader, it is ignored and 
control is returned to the next instruction in the interrupted program. For 
additional details on DEBUG, see the General Loader manual. 


FDBG assumes that the General Loader has placed a debug table in memory. 


RETURN 


Normal return is to the next statement in the interrupted program. 
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Pre-execution Initializer (Batch) 


FUNCTION 


FSTU performs certain initialization functions prior to the execution of 
the user program. The subroutine used depends on the system installation. 


CALLING SEQUENCE 


FSTU is entered by the General Loader only. The entry point is .SETU.; the 
General Loader assumes a six-position storage block (cells .SETU.-1 to 
-SETU.-6 inclusive). 


CELL DEFINITION 
~SETU.-6 Batch/time sharing flag. Set to zero. 
~SETU.=-5 Upper half contains lowest address of memory used by the 


program and the LABELED COMMON region; lower half contains 
the highest memory address used in BLANK COMMON region. 


-SETU.~4 Logical unit table pointer in address field. 


.SETU.-3 Upper half = lowest cell used by the program. Lower half # 
0 = address of pointer to DEBUG subroutine in link 0. 


~SETU.-2 Memory reset constant. 


> SETU.-1l Upper half = entry point address. Lower half # 0 indicates 
low-load job. 


METHOD 


FSTU clears unused memory and sets it to either the constant specified in 
the $ OPTION control card or to zeros if the constant is not defined. It 
places the address of the "logical file - file control block" table in 
fault vector location 25 (octal). In a link job with DEBUG requested in 
link 0, it places a transfer to the DEBUG table in fault vector location 15 
(octal). It places the entry point and a bit indicating a low-load job in 
fault vector location 24 (octal). It then calls the subroutine .FLTPR to 
initialize for fault processing. A secondary SYMDEF, .FLTPR, is imbedded in 
the Execution Error Monitor subroutine to satisfy this SYMREF in a FORTRAN 
execution. This routine places transfers to fault processing routines in 
fault vector locations 7 and 11 (both octal numbers). These are also 
imbedded in the Execution Error Monitor subroutine. There also exists in 
the library a separate subroutine with SYMDEF .FLTPR which zeros fault 
vector locations 7 and 11 (both octal numbers) (divide check and overflow). 
The default option, if either condition occurs, is to abort. This 
subroutine is used if the job uses no FORTRAN library subroutines. Either 
of the .FLTPR subroutines returns to FSTU, which zeros out all index 
registers and performs a TSX1 to the real entry point of the user's 
program. FSTU can be easily changed by an installation to perform the fault 
processing, accounting techniques, etc., for the particular installation. 
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RETURN 


Normal return is to the next executable statement of the calling program. 


RESTRICTIONS 


If the user requires initialization of the cells specified in FSTU, he must 
either use this subroutine or supply one of his own to perform this 
initialization. 
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Pre-execution Initializer (Time Sharing) 


FUNCTION 


FTSU performs the same function for time sharing that FSTU performs for 
batch programs. In addition, FTSU sets up the file control block and 
logical unit table (performed by the General Loader in a batch environment) 
for all files specified in the time sharing RUN command. Cell .SETU.-6 is 
set to nonzero so that it can be used as a test for batch/time sharing. 
Refer to the description of FSTU (batch) for additional details on FTSU. 
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Arithmetic Fault Processor 


FUNCTION 


FFLT zeros fault vectors for overflow, divide check, and underflow faults; 


called by FTSU and FSTU if FXEM is not loaded. (FXEM has its own version of 
-FLTPR that initializes fault vectors.) 


CALLING SEQUENCE 
CALL .FLTPR 


RETURN 


After the execution of the subroutine, control is returned to the calling 
program. 
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Backspace Record 


FUNCTION 


_FBST backspaces one logical record of a file on magnetic tape or sequential 


disk. ' 


CALLING SEQUENCE 


In FORTRAN (indirect) : 
BACKSPACE n 

For GMAP: 

CALL .FBST. (n) 


n = logical file to be backspaced 


METHOD 


3/77 


If the file has buffers and the device is magnetic tape or disk, the 
current record index is examined. If the current record index is for the 
first logical record in a block of logical records, a further check is made 
to determine if partitioned (segmented) records are present. When 
partitioned records are found, the control word of the current record is 
checked to determine if the segment now in memory is the last segment of a 
partitioned record. If so, the file is backspaced until it is positioned 
at the beginning of the first segment of the record. At exit, the "next 
record not there" bit is on, and the file is in the input mode. 


If the current record index is not for the first logical record in the 
block, or if the record is not partitioned, the pointer in the file control 
block and buffer control word are altered. This indicates that the logical 
record preceding the desired logical record has just been read. 


If positioned at end-of-file, the execution of one BACKSPACE positions’ the 
file so that a WRITE statement allows appending to the file, while a READ 
repeats the end-of-file exit if mass store, or an abnormal termination if 
magnetic tape. 


The execution of two BACKSPACE commands positions the file so that a READ 
statement will obtain the last record prior to the end of file. 


2-32 DD20A 


"RETURN 


3/77 


Normal return is to the next executable statement of the calling program, 
except under the following error conditions: 


l. 


Error code 47 if the logical file is assigned to SYSOUT or to a 
device other than magnetic tape or disk. The Execution Error 
Monitor terminates this execution. 


Error code 49 if an erroneous end-of-file appears. The Execution 
Error Monitor terminates this execution. 


Error code 50 if the backspace request is refused because the 
block count in the file control block is zero; i.e., inhibits 
backing out of the file. A message is written on SYSOUT to 
indicate the refusal, and the operation continues. 
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Rewind and Endfile Processor 


FUNCTION 


FEFT, for FORTRAN statement ENDFILE(i), writes an end-of-file record on the 
file. If the file is in input mode prior to executing ENDFILE, FEFT changes 
the mode to output (assuming there is. no error condition). FEFT, for 
FORTRAN statement REWIND(i), rewinds magnetic tape or disk. If the file is 
in output mode prior to executing REWIND, FRWT writes an end-of-file record 
prior to the rewind operation. If the exit from FRWT is without an error 
condition, the file is in the input mode. 


CALLING SEQUENCES 


CALL .FEFT. (arg) compiled for FORTRAN statement ENDFILE (i) 
CALL .FRWT. (arg) compiled for FORTRAN statement REWIND (i). 
arg is the location of 


DEC i 


METHOD 


de ENDFILE - .FEFT. sets the file designator word options to provide the 
following responses: 


a. File initially open and in input mode: FEFT calls GSTOT to change 
the mode to output, and then calls GCLSE to write an end-of-file 
record and close without rewind. 


b. File initially open and in output mode - FEFT calls GCLSE to 
write an end-of-file record and close without rewind. 


on File initially closed - FEFT calls .GOPEN to open the file as 
output without rewind, and then calls .GCLSE to write an 
end-of-file record and close without rewind. 


es REWIND - FRWT sets the file designator word options to provide the 
following responses: 


a. File initially open and in input mode - FRWT calls GCLSE to close 
the file with rewind. 


by File initially open and in output mode - FRWIT calls GCLSE to 
write an end-of-file record and close the file with rewind. 
Following the rewind, the file is set to the input mode. 


Cs File initially closed - FRWT calls GOPEN to open the file as 
input with rewind, and then calls GCLSE to close the file. 
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RETURN 


Normal return is to the next executable statement in the calling program, 
except for error condition 35. Error code 35 applies if there is an attempt 
to rewind or write an end-of-file record if the device assigned to the file 


is SYSOUT or is not magnetic tape or disk. Execution is continued; however 
the endfile/rewind request is ignored. 
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End-of-File (On Input) Processor 


' FUNCTION 


FEOF writes an end-of-file message and either terminates execution or 
returns to the calling program when an end-of-file on input is encountered, 


CALLING . SEQUENCE 


CALL .FEOF. 


It is assumed that the address of the proper file control block is in the 
FOPE subprogram cell .FBAD. 


METHOD 


If a library subroutine detects an end-of-file on an input file, it calls 
FEOF. FEOF places the file number in the message and calls the Execution 
Error Monitor to print the message. 


RETURN 
FEOF selects one of the following: 


1. Terminates execution with error code 34 writing end-of-file 
message. 


2s If the user has provided for end-of-file condition by previously 


calling FLGEOF for the current file, returns to the calling 
program, indicating an end-of-file was encountered. 
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File Opening 


FUNCTION 


FOPE selects and assures that the. physical file associated with the 
specified logical file is open. 


CALLING SEQUENCE 


Ll. 


CALL .FGTFB obtains the file control block address that is in the 
field of cell .FBAD. only. 


CALL .FOPEN(s) obtains the information described in CALL .FGTFB and 
assures that the logical file is open. S indicates the mode in which 
the file is to be opened: 


-1,DL, open the file in its previous mode. 
Odd, open the file as output. 

Even, open the file as input. 

-~1,DL, only bit 35 is examined. 


rh 
13) 
ee 


In both these calling sequences, it is assumed that upon entry to FOPE 
the logical file references are contained in character position 5 of 
cell .FBAD. 


NOTE: When a file is opened, a table is used to find the available 
buffers that can be assigned to the file being opened. The 
table of standard length reusable buffers (321 words) is 
defined as in .FBFTB. 


An equivalent entry to .FOPEN is .FXOP., which is used by the 
Execution Error Monitor to prevent destroying calling sequences in 
case of recursive entry. This entry is used to open P* in order to 
write an error message. 


CALL FLGFRC (lgu, ptr) allows the user to set his error routine 
address in a file control block in the case of a File and Record 
Control error. 


lgu = numeric file code 


ptr address of recovery routine 


See "File and Record Control I/O Error Recovery" in Section V. 
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METHOD 


FOPE is performed in three phases: (1) locating the physical file, (2) 
assigning buffers to the file, and (3) assuring that the file is open. FOPE 
examines the logical file table for a logical file identical to the one in 
cell .FBAD. The file control block. is examined to determine if buffers are 
required that have not previously been assigned. If buffer assignment is 
necessary, table .FBFTB is examined to see if any buffers have been 
released that were previously assigned to another file. These buffers are 
assigned first. If none of these buffers are available, buffers are 
assigned from available unused memory. The file control block is again 
examined to see if the file is open. If it is not open, the File and Record 
Control subprogram OPEN is called for the proper file control block. 
Control is returned to the calling program. 


By the use of FLGFRC, the user can cause word -15 of the file control block 
to point to his error routine. If word -15 is zero when a file is first 
referenced, no change is made to the file control block. However, if word 
-15 is not zero, then a pointer to a translation routine within FOPE is 
stored in word -5, so that a File and Record Control error results in a 
trap to FOPE, an error code 85 message, and a transfer to the user error 
routine in word -15. When running in batch mode, a user error routine can 
be specified by use of a $ FFILE card. This causes the address of the error 
routine to be placed in word-5 of the file control block, thus avoiding the 
trap to FOPE on a File and Record control error. 


RETURN 


Normal return from .FGTFB and .FOPEN is to the next executable instruction 
in the calling program, except for the following error conditions: 


l. Error code 37 if the logical file requested is not in the logical 
unit table. FXER terminates execution. 


2. Error code 38 if there is not enough memory available for I/0 
- buffer assignment. FXER terminates execution. 


36 Error code 56 if there is an attempt to read SYSOUT. FXER 
terminates execution. 


4. Error code 54 if there is an attempt to write GIN. FXER 
terminates execution. 


5. Abort code Q2 if no logical file table exists. 


6. Abort code Q3 if logical file 06 does not exist in the logical 
unit table. A message from FXER cannot be written. 
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Carriage Control Simulator 


FUNCTION 


FSLW formats FORTRAN generated print lines for the printer. 


CALLING SEQUENCE 
CALL .FSLEW(p1) 


pl = location of the print line. It is assumed that location .FBAD. in 
subroutine FOPE contains the address of the file control block for the 
output file, and that word +1 of the file control block contains’ the 
size of the print line. 


METHOD 


FSLW is called by the FORTRAN I/O routine FRWD (I/O Interface by Format 
Control). Control is passed to FSLW after each print line has been prepared 
according to the format specification. Recognized carriage control 
characters are 0, 1, +, and }#, and FSLW looks to see if the first character 
of the prepared print line is one of these. 


If the first character is not a recognized carriage control character, FSLW 
assumes the normal case: single space carriage positioning (J). FSLW 
therefore appends one word of single space slew information to the current 
print line. 


If the first character of the prepared print line is a recognized carriage 
control character, FSLW proceeds as follows: 


l. % - Single Space - This is the normal’case, single space carriage 
positioning (}). 


2. — Space Suppress - The single space slew information appended 
to the previous print line is replaced by space suppress’ slew 
information. 


NOTE: For space suppression control, two conditions are assumed: 

7 (1) that the last record written on the current file is 

the line on which overprinting is desired and (2) that its 

slew information can be changed. Both conditions must be 
satisfied for proper operation. 


Se 1 - Eject Before Printing - A one-word print line (that is, a 
one-word record) consisting of slew-to-top-of-page information is 
generated. This causes a slew to top of page to follow 


immediately after the single space resulting from the information 
appended to the previous line. 
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4. 0 - Double Space ~ A one-word print line (that is, a one-word 
record) consisting of single space slew information is generated. 
This information and the single space information appended to the 
previous print line result in a double space operation. 


Besides taking the actions described above for the recognized carriage 
control characters, FSLW sets any such control character-to a blank (#), if 
it is not already a blank. 


The NOSLEW option on the $ FFILE card causes bit 23 of FCB word -6 to be 
set to l. 


FSLW recognizes this option and changes its normal operation as follows: 
l. The addition of a slew word at the end of a data record is 
inhibited. 


2s The generation of one-word print lines containing only slew 
information is inhibited. 


3. The substitution of a blank character for the carriage control 
character (first character of data record) is inhibited. 


4. A media code 0 is stored in the record control word in place of 
media code 3. 


or The "blank line" records generated by consecutive slashes in 


FORMAT statements are represented by one-word records consisting 


only of blanks. (In the absence of the NOSLEW option, consecutive 
slashes in FORMAT statements cause one-word records containing 


slew characters for single-line slews to be generated. 


RETURN 


Normal return is to the next executable statement of the calling program. 
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SUBROUTINES THAT ARE USER CALLABLE 


Console Communication 


FUNCTION 


FCSL permits operator-program communication via the console; restricted to 
batch mode, not executable in time sharing. 


CALLING SEQUENCE 


CALL CNSLIO (console,message,nwords,nreply,nrepws) 


console = Any of - BCI 1,0000T/ for master console, 
BCI 1,0000T* for tape console, 
BCI 1,0000*T for unit record console, 
BCI 1,0000/T for special purposes. 
If CONSOLE is none of these, BCI 1,0000/T is used. 
message = an array containing one line of the message (in BCI) to be 
output. The message, as received, is prefixed by the SNUMB: 
(SLEW) SSSSS-AA; it will be suffixed by "07701". 
nwords = number of words to be output. Any value greater than 11 is set to 
i oe 
nreply = optional and used if a reply is desired (it need not be any 
particular type). The reply is limited to a maximum of six 
characters unless "nrepws" is supplied. 
nrepws = contains the maximum length in words (up to 11) of the reply. 
Characters are leftjustified in the reply word. 
METHOD 


FCSL outputs the message to the console specified by the calling sequence. 
If the console specification is not one of the four valid specifications, 
the message is output to the master console. This subroutine takes nwords 
(1l maximum) starting at message, prefixes the message with SNUMB and 
suffixes the message with a one line slew, and sends the message to _ the 
console. If a reply is called for (nreply used), the reply is limited to 
Six characters, and is left justified. (Nreply is optional.) 


The program relinquishes control until the console I/O is completed. The 
registers and indicators are restored and a return is made to the next 
executable statement in the calling program. . 
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RETURN 


Return is to the next executable statement of the calling program and is 
done after the console input/output is completed. All registers and 
indicators are restored prior to the return. : 


RESTRICTIONS 


Message may not exceed 11 words. Anything in excess of 11 words is ignored. 
Nwords = 0 results in only the SNUMB being output. 


The reply buffer area (nreply) is not cleared by FCSL. If required, the 
user must clear this area prior to the call to CNSLIO. 
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Memory Dump 


FUNCTION 


FDMP dumps registers and all of memory or designated areas of memory (that 
has been allocated to variables) in a specified format. 


CALLING SEQUENCE 


CALL DUMP (a 1,5,,f£1,,.-+++AnrDnr fn) for BCD 
CALL DUMPA (a1,b1,f1,-..-an,Dnyfn) for ASCII 
CALL PDUMP (a;,6,,f1,-+-+-An-Dns fn) for BCD 
CALL PDUMPA (a, ,b ),£,,+++An-bn, fn) for ASCII 


a and b are variables at the beginning and end of the area to be dumped. a 
or b may represent the first and last variables in the program unit, in 
which case all memory allocated to the variables is dumped. 


f is an integer specifying the dump format as follows: 


te Octal 

Integer 

Real 

Double precision 
Complex 

Logical 
‘Character 


Hout @ on ou 
NU mWNEF © 


If fj is omitted, it is assumed to be zero. If no arguments are given, 
all of memory, including the program object code, is dumped in octal. 


METHOD 


An appropriate NAMELIST table is created, using the parameters specified in 
the calling sequence. FDMP calls FNLO for the actual NAMELIST output 
processing. The panel is dumped, followed by the blocks of memory 
requested. 


RETURN 


If DUMP (or DUMPA) is called, execution is terminated by a call to EXIT. If 
PDUMP (or PDUMPA) is called, the panel is restored and control is returned 
to the calling program. 
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File Control Block and Logical Unit Table Routines 


FUNCTION 


FSET contains three subroutines to allow the user to assign space in memory 
for use as an input/output buffer; to define a file control block for use 
by the input/output subroutines; and to define a logical unit table for use 
by the input/output library subroutines. Each of these three subroutines 
are described separately. 


- DEFINE FILE CONTROL BLOCK 


FUNCTION 


This entry point in FSET, SETFCB, allows the user to define a file control 
block for use by the I/O subroutines. SETLGT must be called first if the 
General Loader has not created any file control blocks. 


CALLING SEQUENCE 


CALL SETFCB(a,i,j...) 


a = location of LOCSYM in the user created file control block 
i,j... = logical files that refer to this file control block 
METHOD 


SETFCB searches the previously defined logical file table for an open space 
to insert the reference to the file control block. It accepts the file 
control block address and appends characters 3,4,5 to the various logical 


file codes referring to this file control block. SETFCB makes as many 
entries as necessary in the “logical file - file control block" table. 


RETURN 


Normal return is to the next executable statement of the calling program 


except when error conditions are encountered. Possible error conditions 
are: 


a Abort code Q2 if there is no logical file table. 


2. Abort code Ql if there is no space available in the logical file table 
for inserting a specified file control block. 
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DEFINE BUFFER(S) FOR A SPECIFIED FILE CONTROL BLOCK 


FUNCTION 


This entry point in FSET, SETBUF, allows the user to assign space in memory 
for use as an input/output buffer. 


CALLING SEQUENCE 


CALL SETBUF (1,a) 


CALL SETBUF (i,a,b) 


i = logical file designator 

a = location of first buffer 

b = location of a second buffer if necessary 
METHOD 


SETBUF searches the logical file table for the specified file and its 
associated file control block. It then attaches the buffers defined to the 
file control block. No check is made to verify that the buffers are of 
sufficient size, as this is the user's responsibility. 


RETURN 
Normal return is to the next executable statement in the calling program. 
RESTRICTIONS 


The size of the total buffer must be one location greater than the area to 


be used for actual record storage. Therefore standard buffer size is 321 
words. 
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DEFINE LOGICAL UNIT TABLE 


FUNCTION 


This entry point in FSET, SETLGT, allows the user to define a logical unit 
table for use by the I/O library subroutine. 


CALLING SEQUENCE 


CALL SETLGT (a,i) 


a = location of logical unit table to be used 
i = number of cells in table a 
METHOD 


SETLGT accepts the array specified by the user as the logical unit table. 
It changes its first location to be a pointer to the last usable position 
of the array and places the address of the array + 1 in fault vector 
location 25 (octal). 


~ 


RETURN 


Normal return is to the next executable statement of the calling program. 


RESTRICTIONS 


SETLGT must be called before any input/output is requested. SETLGT is 
called only when the user wants to suppress the logical file table 
generated by the General Loader and to place the table in his own portion 
of memory. The user should use the NOFCB option on the S$ OPTION General 
Loader control card. 
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File Closing 


FUNCTION 


FCLO closes a file and releases the buffer assigned to that file. The 
buffer is released only if it is standard size (321 words). 


CALLING SEQUENCE 

CALL FCLOSE (u) 

u = logical file number 
METHOD 


FCLO opens the file in its previous mode and calls the File and Record 
Control routine to close the file without rewind. FCLO examines’ and 
releases any standard size buffers assigned to the file. It places’ the 


memory address of these buffers in a table of available buffers (location 
-FBFTB in the FOPE subroutine) for possible reassignment to a newly opened 
file. 

RETURN 


Normal return is to the next executable statement in the calling program. 


RESTRICTIONS 


If more than one logical file refers to one physical file, the physical 
file must be closed, using FCLO only once. 


2-47 DD20 


Initialization of End-of-File Processing _ 


FUNCTION 


FFEE provides a signal to FEOF requesting a return to the calling 
subroutine if an end-of-file condition occurs. 


CALLING SEQUENCE 


CALL FLGEOF (u,v) 


u = logical file number 
v = address of the variable used to indicate an end-of-file condition (user 
must test v when an end-of-file condition could have occurred) 
METHOD 


The address of the variable to be used for end-of-file processing is placed 
in word -15 (upper) of the file control block. The value of the variable is 
set to 0 (set to nonzero if an end-of-file is encountered on logical file 
u). 


RETURN 


Normal return is to the next executable statement in the calling program. 


RESTRICTIONS 


The user must observe the following restrictions on FFEE operation: 


ie If more than one logical file refers to one physical file, the same 


variable must be used for all logical files referring to that physical 
file. 


2% This call must be made prior to any reference to file u. 
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Initialization of Data-Error Processing 


FUNCTION 


FFER provides a variable used in detecting the occurrence of erroneous 
data. 


CALLING SEQUENCE 


CALL FLGERR (u,v) 


u = logical file number 
v = variable used to indicate an input data error (user should test v 
before processing the data because this subroutine gives a normal 
return) 
METHOD 


FFER places the address of the variable in word -16 (upper) of the file 
control block. The value of the variable is initialized to zero (an error 
detected in the input data sets the variable to nonzero). 


spe 


RETURN 


Normal return is to the next executable statement of the calling program. 


’ RESTRICTIONS 


If more than one logical file refers to one physical file, the same 
variable must be used for all logical files referring to that physical 
file. 
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Job Termination. 


FUNCTION 


FXIT contains an entry point .FTERM that is used to terminate the current 
activity. 


CALLING SEQUENCE 


CALL .FTERM 


METHOD 


3/77 


FXIT transforms the logical file table created by the General Loader or by 
FTSU into a file designator word list for closing all files. FXIT calls 
CLOSE which purges all buffers, writes an end-of-file on an output file, 
and notes the closing of an input file. Execution is terminated by a MME 
GEFINI or DRL RETURN. Control is returned to the General Comprehensive 
Operating Supervisor. 


This call will terminate the current activity without checking the wrapup 
list that may have been generated via a CALL ATCALL or CALL NOCALL 
statement. (Refer to Appendix F, FORTRAN Debugging System, in the FORTRAN 
manual.) A CALL EXIT or CALL .FEXIT statement may also be used to 
terminate the current activity and the wrapup list will be inspected. 
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File Forwardspace and Backspace 


FUNCTION 


FFFB allows users who generate multifile tapes to space from one file to 
another (valid only with tape files). 


CALLING SEQUENCE 


CALL FILBSP(xx,n) Backspace n files 
CALL FILFSP(xx,n) Forwardspace n files 


xx = filecode, integer constant or variable 
= number of files to skip; integer variable or constant 


METHOD 


This subroutine is called directly by the user. To ensure proper 


positioning, the current file, if output, should be closed with an ENDFILE 
xx statement and counted as one of the files to be backspaced over. 


There are some file restrictions based on the File and Record Control. 


l. File must be declared multifile. 


2. File must be unlabeled. 
This can be accomplished by using a $ FFILE card as follows: 


$ FFILE XX,MLTFIL,NSTDLB 


RETURN 


Normal return is to the next executable statement of the calling program. 
However, there is one error exit from the subroutine, error code 36. The 
user has requested a backspace larger than the number of files currently 
processed on the tape. The error message is printed and the tape is 
positioned at the first file; execution continues. 
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Call TSS Subsystem 


FUNCTION 


FCAL is used to provide the user with a method to call or 


sharing system subsystem. 


CALLING SEQUENCE 


CALL CALLSS (string) 
or 


CALL CALLSS (string, name) 


go to a time 


time sharing 


a terminating 


string = ASCII character constant or variable that is 
command to invoke the subsystem. The string must contain a 
carriage return or backward slash (octal 134) 
character. 
name = 4-character ASCII name of subsystem to be called 


Return is to the next executable statement in the FORTRAN program. 


CALLGT (string) 
or 


CALLGT (string,name) 


String and name are as for CALLSS. The program terminates rather than 


returning. 


Both calls are ignored in batch. 


METHOD 


DRL CALLSS or DRL T.GOTO are used. 
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Create TSS Temporary File 


FUNCTION 


FDEF is used to create a named temporary file and to access the file in a 
user's available file table (AFT). 


CALLING SEQUENCE 


CALL DEFIL (name, links, mode, istat) 


name = up to 8-character variable or constant containing the ASCII name of 
the temporary file to be created 
links = size in links of file to be created 
mode = 0, sequential file is created 
# 0, random file is created 
istat = status return word as follows: 


0, successful 
3, no room in AFT 


4, temporary file not available 
duplicate file name 


Sy 
6, no room in PAT 


This call is ignored in batch. 


METHOD 


DRL DEFIL is used. 
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Specify Record Size, Random Binary File 


FUNCTION 


FRRD permits the user to specify the record size for a random binary file. 
Normal return is to the next executable statement of the calling program. 
If the record size for a given random file is not provided at load time via 
the $ FFILE card, a call to this subroutine before opening (first I/O to) 
the file is mandatory. 


CALLING SEQUENCE 
CALL RANSIZ (u,n,m) 


logical file designator 
record size 
file format indicator 


How 


u 
n 
m 


u, n and m must be of type integer. They can be any legal arithmetic 
expression. 


Note that a call to RANSIZ can also be used to override a S$ FFILE size 
specification and that this is the preferred method of specification since 
its function works for both batch and time sharing. 


The third argument (m) is optional. When not supplied, file u will be 
processed in standard system format (blocked, variable length records, 
etc.). When supplied, zero indicates standard system format; non-zero 
indicates that block and record control words are not to be processed. This 
latter format provides compatibility with random files generated in the 
time sharing mode. The total file space is available for data; records are 
not blocked, can begin anywhere in a sector and can span device boundaries. 
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Set or Reset Some I/O Parameters of Run-Time Library 


FUNCTION 


This subroutine (FSTU for batch, FTGF for time sharing) permitS-the user to 
set or reset some of the I/O parameters of the run-time library. 
Specifically, it may be used to: 


ds Set the line length (modulo 4) for formatted output directed to a 
terminal. The default setting for this parameter is 72. 


2s Set the media code for unformatted file output. The default setting of 
this parameter is l. 


Bie Set the reflexive read characters that are sent to a terminal to 
request input. The default setting of this parameter is the ASCII 
CHARACTER constant ‘carriage return', ‘line feed', ‘equal sign', X-ON. 


CALLING SEQUENCE 


CALL FPARAM (i,j)> 


i = integer, with a value of 1, 2, or 3 corresponding to one of the three 
functions above 

j = integer, providing the line length or media code for i values of 1 and 
2, Or providing the octal value of four ASCII characters for an i value 
of 3 . 
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SECTION III 


MISCELLANEOUS LIBRARY SUBROUTINES 


This section contains descriptions of library subroutines other than 
input/output. These subroutines are listed in Table 3-1. These subroutines can 
be catagorized as those implicitly called and those that are user callable. 
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Table 3-1. Miscellaneous Library Subroutines 


Function 
as _ 


Double Precision Powers of Ten Table 
Exponent Register Overflow and Divide 
Check Tests 

Sense Light Simulator 

Sense Switch Test 

Restore Links During Execution 

Execution Error Monitor (Common 
Procedure) 

Execution Error Monitor (Batch/Time 
Sharing Procedures and Data) 
Restore Link - H* 

File Transliteration 

Terminal Input Recovery 

ASCII/BCD Indicator 

Date and Time 

Access a Permanent File 

Close File, Detach Buffers, From AFT 

Attach a Temporary Mass Storage or File 
Core Allocator (FDS) 

Special Entry Point (FDS) 

FDEBUG Bootstrap (FDS) 

FDUMP Bootstrap (FDS) 

LINK/LLINK Interface (FDS) 

Delete From Wrapup List (FDS) 
Release Unused Memory (FDS) 

Dummy Setup (FDS) 

Add to Wrapup List (FDS) 

Timing Facility (FDS) 

Wrapup and Loader (FDS) 

Linking Subroutine (FDS) 
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SUBROUTINES IMPLICITLY CALLED 
Double Precision Powers of Ten Table 
FUNCTION 


FDPT stores a table of double precision powers of ten for quick reference 
by all decimal radix conversion routines. 


CALLING SEQUENCE 


-F1DO. is the only SYMDEF symbol. 
METHOD 


FIDO. +2*n is the location of DEC 1.0Dn = 10**n, for n = -34, -37,...,-l1, 
0, 1,..., 37, 38. There are no executable instructions. 
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Restore Link ~- H* 


FUNCTION 


FLHS reloads a program from an H* file (tape) generated in a previous 
General Loader activity. The H* file was generated by a $ TAPE H* control 
card at execution time. 


CALLING SEQUENCE 


This program is called directly from the subroutine library and requires no 
other subprograms. The entire job could be set up as follows: 


SNUMB 

IDENT 

USE -LHSF 

ENTRY .» LHSF 
EXECUTE 

LIMITS 

TAPE js Lara 

DATA AB (Optional) 


MNNNNnMNNNMN 


$ ENDJOB 
*kKEOF 


If the NOFCB option was in effect when the General Loader generated the H* 
file, an entry card of the form: 


$ ENTRY - LHSNF 


must replace the card following the $ USE card. 


METHOD 


The H* file generated by the General Loader contains a link identified as 
//////, Which is the main or common subprogram of the job. If the FCB 
option was in effect during loading (generation of H* file), a second link 
identified as /////1 containing all file blocks generated by the General 
Loader will also be present on. H*. LHSF searches the file for these 
identifiers (/////l is optional), restores them, and enters the main 
subprogram at the entry location specified during the General Loader 
activity. 


Under time sharing execution, the YFORTRAN RUN command restores the main 
link. 


RETURN 


Normal return is to the next executable statement of the calling program. 
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RESTRICTIONS 


The following restrictions apply to FLHS: 


dig 


A $ LOWLOAD card (see General Loader manual) must be included if H* 
file was generated under this option. 


The memory limits in effect when H* file was generated must be 
requested. 


One of the setup subroutines must have been used when H* file was 
generated. Entry to the main link is made through those subroutines to 
initialize fault vectors. 
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Terminal Input Recovery 


- FUNCTION 


FRCV (ASCII) and FRBC (BCD) permit the FORTRAN user to correct a string of 
characters input from a terminal when a character is illegal for the 
current format conversion. 


CALLING SEQUENCE 


CALL .RCOV (curtal,initl)TRA-OK 


curtal = current input tally 
initl = initial input tally 
TRA-OK = recovery transfer point 


If no recovery is possible, the return is 0,1 


METHOD 


When a character is illegal for the current format conversion, the current 
record (current input line) is output with a pointer to the illegal 
character. The user can then input a correction or change (several 
characters) that will replace the corresponding characters previously 
input. The input/output routine will resume with the new string. If the 
user responds with a carriage return, the usual execution error monitor 
message will be output. This routine is not callable by the user. It is 
called by FDIO when appropriate. 
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ASCII/BCD Indicators 


FUNCTION 
FASC is a module that is selected when the ASCII option is used. 


Secondary SYMDEF Entry .ASCB = nonzero 
eASCB+1 = ASCII blanks (040040040040) 
Other ASCII constants are included. 


FBCD is a module that is selected when the BCD option is used. 


Secondary SYMDEF Entry .ASCB = 0 
e-ASCB+1 BCD blanks (202020202020) 
Other BCD constants are included. 
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SUBROUTINES THAT ARE USER CALLABLE 


Exponent Register Overflow and Divide Check Tests 


FUNCTION 


FDCK tests the Fault Status Word for a previous exponent register overflow 
or divide check. 


CALLING SEQUENCE 


CALL DVCHK(j) for divide check 
CALL OVERFL(j) for exponent register overflow 


j} = location of integer variable 


METHOD 


This routine checks the fault vector area and sets an integer variable to l 
if a fault has occurred, or 2 if not. The Fault Status Word (31 octal) is 


not maintained in the time sharing environment. 


RETURN 


Normal return is to the next executable statement in the calling program. 
This routine assumes normal recovery from exponent register overflow and 
divide check. 
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Sense Switch Test 


FUNCTION 


FSWI tests the GCOS switch word for the status of a sense switch. 


CALLING SEQUENCE 


CALL SSWTCH(i,j) to test sense switch i. 
i integer variable or constant, value between 1 and 6 
j integer variable to be set to 1 if sense 

switch i is ON, or 2 if it is OFF 


METHOD 


Bits 6-11 of the GCOS switch word correspond to the sense switches 1-6, 
with 0 denoting OFF and 1 denoting ON. Sense switch specified by i is 
tested and the integer in j is set to 1 if i was ON or to a 2 if i was OFF. 


RETURN 


Normal return is to the next executable statement in the calling program. A 
possible error condition is error code 53. If i is not 1-6, sense switch i 
is declared OFF. 
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FUNCTION 


FLNK enables the programmer to call program overlays in batch mode. 


CALLING SEQUENCES 


CALL LINK (name) 


CALL LLINK (name) 


name = link identifier specified as a character variable or character 
constant 
METHOD 
FLNK assumes that the General Loader has created a file (file code H*) 


containing the user's program segmented into links as specified by $ LINK 
control cards. (H* file is generated by GESAVE.) Both LINK and LLINK use 
the GERSTR function of GCOS. The procedure for restoring a link depends on 
the entry used as follows: 


LINK - Restore the link and transfer to its entry point as specified 
at load time. 


LLINK - Restore the link and return to the next statement or 
instruction in the calling subroutine, 


If DEBUG is requested at load time in any or all of the links, these 
subroutines join the respective DEBUG tables, enabling the user to take 
snap dumps of any links in memory at the time of his request. 


After restoring a link, the link is tested to determine if it contained a 
DEBUG table. If it contained a DEBUG table, the address of this table is 
chained to existing tables. If there was no DEBUG table in memory at this 


time, the address of this table is placed in the DRL cell (cell 13 of the 
fault vector). DEBUG tables, corresponding to links which are overlayed in 
the process, are deleted from the chain. 


The link overlay names must be five to eight characters for ASCII mode 
executions. 


RETURN 


Returns are defined under Method. 
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Restore Links During Execution (Time Sharing) 


FUNCTION 


FTLK performs the same functions for time sharing programs as FLNK performs 
for batch programs. Refer to FORTRAN manual for specific details. 


i ee The YFORTRAN time sharing subsystem RUN command is used to load the 
main program link from an H* permanent file for execution. 


2% The overlay link names on the H* must be five to eight characters when 
running in the time sharing ASCII mode. 


3 The DEBUG option is not available under time sharing execution. 
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Execution Error Monitor 


FUNCTION 
The Execution Error Monitor performs the following functions: 


bs Prints a trace of subroutine calls, if applicable. 


ds Prints execution error messages. 

a Terminates executitn with a Q6 abort or does one of the following: 
a. Continues with execution of the program. 
De Transfers to an alternate error routine. 

4. Allows the user to determine if an error has been processed by 


Execution Error Monitor. 


the 


The Execution Error Monitor functions are performed by the following 


subroutines: 


FXER contains the procedures that are common to both batch and 
sharing operations. 


time 


FXEM contains the procedures and data peculiar to batch operations. 


FXMA contains the procedures and data peculiar to time sharing 


Operations. 


Execution error monitor functions are optional and are determined by three 
groups of switch words (presently there are four words per group) in which 
each bit corresponds to an error code. The first group (.FXSW1) controls 


termination; the second (.FXSW2) controls message printing and trace; 


and 


the third (.FXSW3) controls the alternate error returns. Table 3--2 contains 


the error codes and default procedures. 


SWITCH WORD GROUPS 


1, -FXSW1 (termination). Table 3-2 shows the standard bit settings 


and 


the routines that use the corresponding error codes. The bit settings 


are defined as: 


C (Continue execution) 
A (Terminate with a Q6 abort) 


Termination may be overridden by the corresponding bit of .FXSW3. 


DD20 


-FXSW2 (message printing. and trace) - The meaning of the bit setting 
is 


1 - Suppress printing 
0 - Print 


This group is initialized to zero. Settings may be changed by a 
program call to FXOPT. 


-FXSW3 (alternate error return) - The meaning of the bit setting is 
1 - Alternate error return (overrides termination set in .FXSW1). 
0 - Use normal return. 


This group is initialized to zero. Settings may be changed by a 
program call to FXOPT. 


GMAP CALLING SEQUENCE 


CALL 


* 
Heol 


.FXEM. (x,y) 


address of error description controls 
address of tally word used to indicate card column found in error; y 
is optional and is used only if a card image is to be printed 
Instruction sequence at x 
X ZERO A,B 

ZERO C,D 

ZERO E,F 


A = Address of card image to be printed, or zero when card image is 
not to be printed. 


B = Error code expressed as an integer(n) in the range 1_n 143. 


C = Address of message l. 


D = Word count of message l. 
E = Address of message 2. 
F = Word count of message 2. 
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FORTRAN CALLING SEQUENCES 


ues 


CALL FXOPT (ncode,il,i2,1i2) 


FXOPT is an entry to .FXER anc may be called to alter the standard 
switch word settings. In the statement FXOPT (ncode,il,i2,i13), ncode is 
an error code, and il, i2, and i3 provide the settings for the 
corresponding bits in the three switch word groups. 


Examples: 


Le. GABE -FROPT(32,07.15.0) 
2. CALL FXOPT (32,1,6,0) 
3. CALL FXOPT (32,0,0,1) 


NOTE: Error code 32 denotes an illegal character in input data. 


Example 1 causes a Q6 abort when the error occurs, and no message or 
trace is printed. 


Example 2 causes execution to continue after message and trace are 
printed. 


Example 3 indicates that return is to an alternate error routine after 
trace and message are printed, since the alternate return option takes 
precedence over termination. 


CALL FXALT (SR) 


FXALT is an entry to .FXER and may be called to set the alternate 
error return location. The statement CALL FXALT(SR) communicates the 
name SR of the alternate error routine to the execution error monitor. 
An EXTERNAL SR must be included in the calling routine. If the 
alternate return option for an error code is indicated but no call to 
FXALT has been made, a Q5 abort follows when the error occurs. A 


RETURN statement in the alternate routine continues execution at the 
instruction immediately following the one where the error occurred. 


The statement CALL FXALT($n) designates statement n in the calling 
program as the alternate error return. 


NOTE: If the same error occurs in the alternate error routine, an 
interminable loop results. 
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3% Overflow and divide check fault test 


The fault processor processes divide check, overflow, exponent 
overflow, and exponent underflow faults. A message is output on file 
06 stating the type of fault and the location at which the fault 
occurred. Execution continues in the normal manner, although the EAQ 
registers may have been reset as depicted in the following table. 


FAULT EAQ REGISTERS 
Exponent overflow Large floating-point valuel 
Divide check (FP) Large floating-point valuel 
Exponent underflow Floating-point zero 
Overflow (Integer) No change 
Divide check (Integer) No change 


To have another value returned in the EAQ registers after a divide 
check, CALL FXDVCK(r,m) should be executed prior to the occurrence of 
the fault. This statement causes the value of r to be returned in the 
EAQ registers after a real divide check and the value of m to be 
returned in the Q register after an integer divide check. The first 
argument must be double precision. The second argument may be omitted. 


4, CALL FXEM (ncode,msg,n) 


A FORTRAN-callable entry has been provided so that it may be called 
when the program detects an error condition. 


This statement causes the printing of an error trace and the Hollerith 
message contained in the msg array. The number of words (n) to _ be 
printed must be within the limits 0< ng¢ 20. If only the first argument 
is given, only the trace is printed. 


Se CALL ANYERR (v) 


The user may desire to detect some error that has occurred in an 
input/output routine or a mathematical FUNCTION routine. 


v is a variable into which .FXER will initially place the value zero. 
If an error occurs, the error code is placed in v. The logical IF 
statement provides a suitable means of testing if v is typed LOGICAL. 


lallows further computations without another immediate fault. This value is set 
to approximately 10%**36, 
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METHOD 


Ls 


Error linkage for tracing calls is generated by the Macro Assembler 
Program (GMAP) and by the FORTRAN compiler. Tracing stops when the 
address of the CALL instruction in the error linkage word is zero, or 
when the number of traces exceeds a constant. 


The error trace prints in reverse order. It includes the name of each 
calling routine, identifying number of the CALL instruction, absolute 
location of the CALL instruction, and up to five calling arguments. 


The functions of this routine are optional. The options are controlled 
by the following switch word groups: 


-FXSW1 - Termination 
~-FXSW2 - Message printing 
-FXSW3 ~- Alternate error returns 


Each of the bits (1-143) in a switch word group corresponds to an 
error code. 


Special processing applies to error code 55. When this error is 
encountered, the following message is written: 


ILLEGAL VALUE FOR COMPUTED GO TO AT ID NUMBER XXXXX 


The error code is always stored in the location FXCODE in FXER. Since 
this is a SYMDEF, it may be accessed by a GMAP program. 


The error code is also stored indirectly through a pointer defined in 
FXER. This pointer may be set by calls to ANYERR. If this pointer has 
been initialized to contain the address of a variable in the user's 
program via a call to ANYERR, the variable will contain the error 
code, expressed as an integer, upon return to the calling subprogram, 
after an error. 


FXOPT is an entry to .FXER. which, for a given error code, sets the 
corresponding bits in .FXSWl, .FXSW2, and .FXSW3 to the low-order bit 
of the second, third, and fourth arguments. The first argument is the 
error code. When a call is made to .FXER., the error code is used to 
shift each switch word group and set the options accordingly. 


FXALT stores the location of its argument in location FXALT]1 in .FXER. 
If the alternate error return option is used, index register 1 and the 
indicator register are restored; and a transfer is made to FXALT1 
indirect. Thus, if the alternate return is a subprogram, the RETURN 
statement transfers to the location following the call to .FXER. If no 


alternate has been supplied, a Q5 abort occurs. 
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RETURN 


A divide check, an overflow, or an underflow transfers to .FXER. via 
the program fault vector. (For a description of the fault vector, see 


the manual General Sie ea Ee Operating Supervisor.) .FXER. writes 
the error message an loads e proper values into e EAQ~-registers. 


The normal return is RET 6 (divide check) or RET 8 (overflow and 


underflow). If an alternate return is requested, the indicators and 


index register 1 are loaded from the fault vector, so that a RETURN 
statement in the alternate routine will transfer to the location 
immediately following the one that generated the fault. 


FXEM is the entry provided for error conditions detected by the user's 


. program. Error codes 61-66 are reserved for users. The statement 


CALL FXEM (ncode,msg,n) 


prints an error trace and n words of the message in the array msg. Msg 
must be an array containing character information. If either msg or n 
is omitted or is zero, no message is printed. If n is greater than 20 
words, Only 20 words are printed. 


Error codes and returns are defined in Table 3-2. 
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ocdd 


DEFAULT 
PROCEDURE 
ABORT/ 


CONTINUE 


A 


C 


FUNCTION ERROR 
Not used 
Le * J I=0,J=0 
Ley I=0,35<0 
DA** J DA=0 , J=0 
Ate A=0,0=0 
Ate A=0,J<0 
DA** I DA=0 , J<0 
BE*C B<0,C=0 
At*B A=0, B=0 
A**C A=0 ,C<0 
e**B B>88.028 
LOG (A) A=0 
LOG (B) B<0 
ARCTAN(A/B)  A=0,B=0 
SIN (A) [ajeee? 
COS (A) 
rE: B<0 
CA**K CA=0, K=0 


EXCEPTION 


RETURN 


0 —+ OR 
237 29 -SOR 
0—+EAQ 
10+ EAQ 
0 — EAQ 
0—-EAQ 
10°2- EAQ 
10°28. BAQ 


10°28 FAQ 
0—-EAQ 
0O-—+FAQ 


0—+EAQ 


Jem 


0 —AQ 


Table 3-2. Error Codes And Returns 


MESSAGE 
LINE 1 


EXPONENTIATION ERROR 
OE*FO 


EXPONENTIATION ERROR 
O** (-J) 


EXPONENTIATION ERROR 
OFF) 


EXPONENTIATION ERROR 
Q** (-~J) 


EXPONENTIATION ERROR 
(-B) **C 


EXPONENTIATION ERROR 
Ce*O 


EXPONENTIATION ERROR 
O** (-C) 


EXP (B),B GRIT THAN 88.028 
NOT ALLOWED 


LOG (0) NOT ALLOWED 
LOG (~B) NOT ALLOWED 
ATAN2 (0,0) NOT ALLOWED 


SIN OR COS ARG GRT TH 
2**27 NOT ALLOWED 


SORT (-B) NOT ALLOWED 


EXPONENTIATION ERROR 
0**0 


MESSAGE 


SET 


SET 


SET 


RESULT=0 


RESULT=2** 35=2 


RESULT=0 


RESULT=10** 38 


RESULT=0 


RESULT=0 


RESULT=10**38 


RESULT=10** 38 


RESULT=10** 32 


ae, 


RESULT=0 .C 
RESULT=0 


RESULT=0 


EVALUATE FOR +B 


SET 


RESULT=0 


OZ-€ 


0cdd 


ERROR 

CODE_ 
15 
16 
17 
18 


19 


20 
21 
22 


23 


24 


25 


26 


27 


28 


29 


30 


DEFAULT 


PROCEDURE 


ABORT/ 


CONTINUE 


C 


YO A A A 


FUNCTION 


CA**I 
DA**DB 
DA**DB 
DA**DB 
e**DA 


LOG (DA) 


LOG (DA) 


VDA 


SIN DA 
COS DA 


ARCTAN (DA/DB) 


CA/CB 
e**CA 
e**CA 
LOG (CA) 


SIN (CA) 
COS (CA) 


COS (CA) 


Table 3-2 


ERROR 


CA=0,J<0 
DA<0, DBO 
DA=0 , DB=0 
DA=0 , DB<0 
DA > 88.028 


DA=0 
DA<0 
DA<O 


|pa|> 227 


DA=0 , DB=0 


CB= (0,0) 


REAL CA>88.028 


JIMAG cal > 277 


CA= (0,0) 


|REAL (CA)| > 2’ 


IMAG (CA)>88.028 


(cont). Error 


EXCEPTION 
RETURN 


10° oe aR 
0 —+ OR 


0—-FAQ 
0—EAQ 
10° BAO 


10° 2= FAG 


-10°% rag 
O—-FAQ 


a= bar 


0-—+EAQ 


Codes And Returns 


MESSAGE 
LINE 1 


EXPONENTIATION ERROR 
O** (-J) 


EXPONENTIATION ERROR 
(-DA) **DB 


EXPONENTIATION ERROR 
OF*0 


EXPONENTIATION ERROR 
0** (-B) 


EXP (B),B GRT 88.028, 
NOT ALLOWED 


DLOG (0) NOT ALLOWED 
DLOG (-B) NOT ALLOWED 
SORT (-B) NOT. ALLOWED 


DSIN OR DCOS ARG GRT 
2**54 NOT ALLOWED 


‘DATAN2 (0,0) NOT ALLOWED 


COMPLEX 2/0 NOT ALLOWED 


EXP (Z),REAL PART GRT 
88.028 NOT ALLOWED 


EXP (Z),IMAG PART GRT 
2**27 NOT ALLOWED 


CLOG (0) NOT ALLOWED 


CSIN OR CCOS ARG WITH 
REAL PART GRT 2**27 NOT 


ALLOWED 


CSIN OR CCOS ARG WITH IM 
PART GRT 88.018 NOT 
ALLOWED 


MESSAGE 
LINE 2 


SET RESULT=(10**38,0.0) 


SET RESULT=0 
SET RESULT=0 
SET RESULT=10**38 
SET RESULT=10** 38 


SET RESULT=- (10**38) 
SET RESULT=0 
EVALUATE FOR +B 


SET RESULT=0 


SET RESULT=0 


SET RESULT=(10**38, 
LO**33) 


SET RESULT=(10**38, 
1L0** 38) 


SET RESULT=0,0 
SET RESULT 
(= ({LO**38))°5:0.40) 


SET RESULT=0 


SET RESULT=(10**38, 
LO**38) 


Tot 


0cdad 


ERROR 
CODE 
31 


32 


33 


34 


35 


36 


37 


38 


39 


40 


41 


42 


43 


44 


DEFAULT 
PROCEDURE 
ABORT/ 


CONTINUE 


(es 


FUNCTION 


BCD I/O 


BCD I/O 


LINKED 
BINARY I/O 


BCD I/O 
REWIND AND 
END FILE 
PROCESSOR 


FFFB 

FILE OPENING 
FILE OPENING 
LINKED 
BINARY I/O 
LINKED 
BINARY I/0 
LINKED 
BINARY I/O 


NAMELIST 
INPUT 


NAMELIST 
INPUT 


NAMELIST 
INPUT 


Table 3-2 


ERROR 


ILLEGAL FORMAT 
STATEMENT 


ILLEGAL CHARACTER 


IN DATA OR BAD 
FORMAT 


ATTEMPT TO READ 
OUTPUT FILE 


END-OF-FILE 


ILLEGAL REQUEST 


BACKSPACE ERROR 
FILE NOT DEFINED 
NO SPACE FOR I/0 


BUFFERS 


ILLEGAL END- 
OF-FILE 


LIST EXCEEDS 
LOGICAL RECORD 
LENGTH 


SYSOUT/FIXED 
LENGTH RECORDS 


ILLEGAL HEADING 
CARD 


ILLEGAL VARIABLE 
NAME 


ILLEGAL SUBSCRIPT 


OR ARRAY SIZE 
EXCEEDED 


(cont). Error Codes And Returns 


MESSAGE 
LINE 1_ 


EXCEPTION 
RETURN 


FORMAT AT LLLLLL, HAS 


ILLEGAL CHAR 

ILLEGAL CHAR IN DATA 
OR BAD FORMAT 

READ AFTER WRITE 

IS ILLEGAL 


END OF FILE READING 
FILE CODE FC 


REQUEST TO XXXXXX ON FC 
WAS IGNORED 
TAPE POSITIONED AT FIRST 


FILE 


FC XX 
DOES NOT EXIST 


INSUFFICIENT CORE AVAIL- 
ABLE FOR BUFFERS 


UNEXPECTED EOF OR BAD 
FORMAT 


LIST EXCEEDS LOGICAL 
RECORD LENGTH 
SYSOUT/FIXED LENGTH 
RECORDS 


ILLEGAL HEADING CARD 
BELOW 


ILLEGAL VARIABLE 
NAME BELOW 


ILLEGAL SUBSCRIPT BELOW, 
OR DATA EXCEEDS VARIABLE 


MESSAGE 
LINE 2_ 


TREAT AS END 
OF FORMAT 


TREAT ILLEGAL 
CHAR AS ZERO 


FC XX 


OPTIONAL RETURN 
NOT REQUESTED 


BACKSPACE REQ. LARGER 
THAN FILE COUNT 


FC XX 


STORE ZEROS IN 
REMAINING List 


ITEMS MAY NOT BE 
PARTITIONED-FC XX 
SCAN TERMINATED 
SKIPPING TO NEXT 
VARIABLE NAME 


SKIPPING TO NEXT 
VARIABLE NAME 


COWE 


0cdad 


ERROR 
CODE 


45 


46 


47 


48 


49 


50 


51 


52 


53 


54 


5S 


56 


58 


DEFAULT 
PROCEDURE 
ABORT/ 


CONTINUE 


C 


FUNCTION 


NAMELIST 
INPUT 


NAMELIST 
INPUT 


BACKSPACE 


_ RECORD 


NAMELIST 
INPUT 


BACKSPACE 
FILE 


BACKSPACE 
FILE 


SENSE LIGHT 
SIMULATOR 


NAMELIST 
INPUT 


SENSE SWITCH 
TEST 


FILE OPENING 


FXEP 


FILE OPENING 


BCD 1/0 


BACKSPACE 
RECORD 


Table 3-2 


ERROR 


ILLEGAL CHARACTER 


AFTER RIGHT 
PARENTHESIS 


ILLEGAL CHAR IN 
DATA 


FILE CANNOT BE 
BACKSPACED 


ILLEGAL LOGICAL 
CONSTANT 


ERRONEOUS 
END-OF-FILE 


BLOCK COUNT 
OF ZERO 
INDEX NOT 
O0<n<35 
ILLEGAL 


CHARACTER FIELD 


INDEX NOT 
l<n<6 


ATTEMPT TO WRITE 
[* 


ILLEGAL VALUE 


ATTEMPT TO READ 
pk 


ILLEGAL CHAR 
FOR L CONVERSION 


EXCEPTION 


Error Codes And Returns 


MES SAGE 
LINE 1 


ILLEGAL CHAR IN DATA 
BELOW 


ILLEGAL CHAR IN DATA 
BELOW 


FILE CODE XX, BACKSPACE 
REFUSED 


ILLEGAL LOGICAL CONSTANT 


APPEARS BELOW (OR AT 
END OF PRECEDING RECORD) 


END-OF-FILE ON READ BUT 
NOT ON PREVIOUS BACKSPACE 
OF SAME TAPE RECORD 

BLOCK COUNT = 0 


REFERENCE TO NON-EXISTENT 
SENSE LIGHT 


ILLEGAL CHARACTER 
FIELD BELOW 


NON-EXISTENT SENSE SWITCH 
TESTED 


ILLEGAL WRITE REQUEST 
ON SYSIN1 


ILLEGAL VALUE FOR COM- 


PUTED GO TO AT ID NUMBER 
XXXXX 


IT IS ILLEGAL TO READ 
FROM SYSOUT FC XX 


ILLEGAL CHAR FOR L 
CONVERSION IN DATA BELOW 


FILE NN IS CLOSED 


MESSAGE 
LINE 2 


ASSUME COMMA PRECEDES 
CHAR 


TREAT CHAR AS ZERO 


FILE IS SYSOUT OR IS 
NOT MAG TAPE, D/D 


TREAT ILLEGAL LOGICAL 
CONSTANT AS F 


DECLARED OFF IF 
TESTING IGNORED 
IF SETTING 


SKIPPING TO NEXT 
VARIABLE NAME 


SWITCH DECLARED 
OFF 


NO OPTIONAL EXIT 
EXECUTION TERMINATED 


XXXXX 


TREAT ILLEGAL 
CHARACTER AS SPACE 


Se Aas =. 


o0cdad 


ERROR 


CODE 


71 


72 


73 


74 


75 


76 


DEFAULT 


PROCEDURE 


ABORT/ 


CONTINUE _ 


Cc 


FUNCTION 


NAMELIST 


INPUT 


T**J. 


Table 3-2 (cont). Error Codes And Returns 


EXCEPTION 
RETURN 


ERROR 


EMPTY CHARACTER 
FIELD 


J>36,J3 IS EVEN 
J>36,35 IS ODD,I>0 
J>36,J IS ODD,I<0 


RESERVED FOR USERS 


FAULT 


FAULT 


FAULT 


FAULT 


FAULT 


RANDOM 
BINARY 


RANDOM 
BINARY 


RANDOM 
BINARY 


RANDOM 
BINARY 


RANDOM 
BINARY 


I/O 


70 


770 


I/O 


I/O 


EXPONENT UNDERFLOW 
INTEGER OVERFLOW 
EXPONENT OVERFLOW 


INTEGER DIVIDE 
CHECK 


FLOATING POINT 
DIVIDE CHECK 


LIST EXCEEDS LOGICAL 
RECORD LENGTH 


FILE NOT STANDARD 

SYSTEM FORMAT. ZERO 
BLOCK COUNT; BSN 

ERROR; ZERO RECORD COUNT 


NO FCB FOR FILE 


BAD RANDOM RECORD REFERENGE 


RECORD SIZE NOT SPECIFIED 
IN FCB. GIVE VIA $ FFILE 
CARD OR CALL RANSIZ 

(FC, SIZE) 


MESSAGE 
LINE 1 


EMPTY CHARACTER FIELD 


235 ~2--OR EXPONENT > 35 OR 
2° =2—-OR EXPONENTIATION OVERFLOW 
~ (2°? -2)-QR 


EXPONENT UNDERFLOW 
OVERFLOW 
EXPONENT OVERFLOW 


DIVIDE CHECK 


DIVIDE CHECK 


LIST EXCEEDS LOGICAL 
RECORD LENGTH 


FILE NOT STANDARD 
SYSTEM FORMAT 
PC xx 


LOGICAL FILE CODE XxX 
DOES NOT EXIST 


ZERO OR NEGATIVE 
RANDOM REC # 


REC SIZE NOT GIVEN 
FOR RANDOM FILE 


MESSAGE 
LINE 2 


TREAT AS BLANKS 


SET RESULT = 
#72 ((2**35)=2) 


AT LOCATION XXXXXX 
AT LOCATION XXXXXX 
AT LOCATION XXXXXX 


AT LOCATION XXXXXX 


AT LOCATION XXXXXX 


STORE ZEROS IN 
REMAINING LIST ITEMS 
FC XA 


NO OPTIONAL EXIT 
EXECUTION TERMINATED 


FC XX 


FC XX 


LL/E 


VOCE 


wocdd 


ERROR 


CODE 


77 


78 


79 


80 


81 


“82 


83 


84 


85 


86 


87 


DEFAULT 
PROCEDURE 
ABORT/ 


CONTINUE 


A 


FUNCTION 


RANDOM 
BINARY I/0 


RANDOM 
BINARY I/0 


RANDOM 
BINARY I/0 


RANDOM 
BINARY I/0 
FORMAT I/O 
ENCODE /DECODE 


FORMAT I/O 
ENCODE/DECODE 


ARCSINE 
FORMAT 1/0 
ENCODE /DECODE 
I/O 

FORMAT I/0 
ENCODE /DECODE 


I/O 


Table 3-2 (cont). 


EXCEPTION 
ERROR RETURN 


RANDOM WRITE TO LINKED 


FILE ILLEGAL. LINKED 


FILE MAY BE READ RANDOMLY 


BUT NOT WRITTEN TO. 


THE RECORD NO. GIVEN IN 
THE RANDOM READ OR WRITE 
STATEMENT IS OUTSIDE THE 


FILE LIMITS. 


LIST EXCEEDS DECLARED 
RECORD LENGTH. 


FILE IS NOT LARGE 
ENOUGH TO CONTAIN 
RECORD 


LINE EXCEEDS SIZE OF 
RECEIVING FIELD 


FIRST NON-BLANK CHAR- 
ACTER IS NOT ( 


1 ARG I> 1.0 


| INTEGERI> 2**35-1 


"“GFRC" ERROR 
ENCODE/DECODE-I/0 
MAY NOT BE USED 

RECURSIVELY 


SPACE/CORE OBTAINED 


Error Codes And Returns 


MESSAGE 
LINE 1 


RANDOM WRITE TO 
LINKED FILE ILLEGAL 


REC # OUT-OF-BOUNDS- 


LIST EXCEEDS DECLARED 
RECORD LENGTH 


FILE SPACE EXHAUSTED- 


LINE EXCEEDS SIZE OF 
RECEIVING FIELD 


FIRST NON-BLANK CHAR~ 
ACTER IS NOT ( 


1! ARGI> 1.0 

| INTEGERI> 2**35-1 
“GFRC" ERROR 
ENCODE/DECODE= 


I/O MAY 


SPACE/CORE OBTAINED 
FOR 


MESSAGE 
LINE 2: 


FC XX 


FC XX 


FC XX 


FC XX 


TREAT AS END 
OF FORMAT 


TREAT AS END 
OF FORMAT 


EVALUATE FOR 
ARG=1.0 


LIMIT TO 
2**35-1 


FC XX 
NOT BE USED 
RECURS IVELY 


LOG. FILE 
CODE XX 


LL/E 


T° PZ-€ 


Wocad 


ERROR 


CODE 


88 


89 


90 


91 


“92 


93-99 


DEFAULT 
PROCEDURE 
ABORT/ 


CONTINUE 


Cc 


Table 3-2 (cont). Error Codes And Returns 


EXCEPTION 

FUNCTION ERROR RETURN 
CALLSS END OF STRING 

CHARACTER MISSING 
EXP UNDERFLOW 
DEXP 
TAN ARG TOO LARGE 
DTAN 
ACOSH ILLEGAL ARG 
DACOSH 
ATANH ILLEGAL ARG 


NOT PRESENTLY USED 


NOTATION: I,J,K are integers 
A,B,C are real numbers 
DA,DB,DC are double-precision numbers 


MESSAGE 
LINE 1 


EXP (TOO LARGE A 
NEGATIVE NUMBER) 


LARGE ARG(71E4) 
TO TAN 


ACOSH OF NUMBER .LT. 
1.0 NOT ALLOWED 


IXi .GE. 1.0 TO 
ATANH (X) 


CA,CB,CC where CA=X,Y are complex numbers 


MESSAGE 


LINE 2 


SET RESULT 
=0.0 


MAY CAUSE LOSS 
OF PRECISION 


SET RESULT 
TO 0.0 


SET RESULT TO 
+ CR -10**38 


File Transliteration 


FUNCTION 


FMED is a file media transliteration module. 


CALLING SEQUENCE 


CALL FMEDIA (fc,media) 


fe 


i} 


FORTRAN logical file code 


media = pointer to a code representing the form of the output as follows: 


il 


, BCD NSLEW 

» BCD. cards 

, Printer 

, Time Sharing ASCII (obsolete) 

, Standard System Format ASCII (no slew) 
All others are ignored 


media 


it tt 


NWN © 


The legal combinations are as follows: 


0 to 2 3 to 0 
0 to 3 3 to 2 
0 to 5 3 to 5 
0 to 6 3 to 6 
2 to 0 6 to 0 
2 to 3 6 to 2 
2 to 5 6 to 3 
2 to 6 6 to 5 


METHOD 


FMED sets bits 18-21 of the file control block LOCSYM+5. These bits are 
used to inform the I/O edit function that each output record directed to. 
the file (fc) must be transliterated to the format represented by the media 
in bits 18-21 of FCB+5. FMED calls FOPEN to locate the file control block. 
FOPEN returns the file control block location. FMED calls GMEDIA. If media 
= 0, the nslew bit in the file control block is also set. 


Automatic file transliteration is provided and/or reformating on a_ logical 
record basis permits the following: 
dy Executing of a BCD program under time sharing. 

ae I/O can be directed to the terminal. 

b. Input files can be ASCII (media 5 or 6). 


Ce Output files can be media 0,2,3 BCD or 5,6 ASCII. 
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Execution of an ASCII program in batch. 


a. I/O can be directed to reader, printer, punch or 


b. Input files can be media 0,2,3 BCD or 5 ASCII. 


c. Output files can be media 0,2,3 BCD 


Execution of a BCD program in batch. 


or 6 ASCII. 


a. Input files can be ASCII (either media 5 or 6). 


b. Output files can be media 0,2,3 BCD 
Execution of an ASCII program under time 
a. Terminal I/O is provided. 

b. Input files can be media 5 ASCII or 


Cc. Output files can be media 0,2,3 BCD 


or 6 ASCII. 


sharing. 


0,2,3 BCD. 


or 5,6 ASCII. 


DD20 


Date and Time 


FUNCTION 


FDTM allows a user to obtain the current date and time. 


CALLING SEQUENCE 


CHARACTER A*8 
REAL B 


CALL DATIM (A,B) 


Upon return, A will contain the date in the form mm/dd/yy (with trailing 


blanks if in BCD mode); B will contain the time-of-day in hours as a 
floating point binary number. 


METHOD 


FDTM uses the MME GETIME in the batch mode and the DRL TIME in the time 


sharing mode. FDTM uses the GCVT module for transliteration on the date, if 
required. 


Cell .SETU.-6 is used to determine if the environment is batch or time 
sharing. 
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Access a Permanent File 


FUNCTION 


FTAC is used to access an existing permanent file. 


CALLING SEQUENCE 


CALL ATTACH (lgu,catfil,iprmis,mode,istat,buffer) 


3/77 


lgu = 


catfil 


iprmis = 


the FORTRAN file code (an integer expression, variable, or 
constant) 

a character constant, or variable, containing the catalog/file 
string. It must be terminated by a semicolon; embedded blanks 
are ignored. The user master catalog password is used if it 
exists; the system master catalog password is never used; 
however, subsequent passwords are required if they are part of the 
file description. 

the permissions desired. These will be ORed with any permissions 
in the catfil, 


If iprmis = 1, READ only; if iprmis = 2, WRITE only; if iprmis =3, READ and 


mode 
mode 
mode 
mode 
istat 


Hou uw uw ai 


buffer 


WRITE; otherwise, undefined and may change. 


an integer variable or constant 

0; Gets file as defined 

1; Gets file as random 

2; Gets terminal 

the status return from the file system (see the TSS System 
Programmer's Reference Manual for TSS codes), or will contain a 
status as follows: 


OK (batch mode only) 

File is currently open 

Terminal requested in batch mode (illegal) 

Additional memory needed, request denied (time sharing user is 
aborted) 

catfil all blanks 


WHR O 


4 


Wi 


null arg: get a file system buffer. 

not null arg: use this variable array as a buffer (at least 380 
words). 

Example of null arg: 


CALL ATTACH (lgu, "catfil;",iprmis,mode,istat, ) 
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METHOD 


Upon successful return from FTAC (attach), a file control block will have 
been created and the file name (or alternate name) will be in FCB -10, -9 
(in ASCII). If the file was in the available file table (AFT), it will be 
deaccessed and reaccessed with the new permissions, if necessary. 


bo 
i 

bo 

Ne) 


2DZO 


Close File, Detach Buffers, Remove from AFT 


FUNCTION 


FDTH is used to close a file and release its buffers. In time sharing, the 
file is also removed from the AFT. If more memory is needed (to deaccess 
the file) and the request is denied, the time sharing user is aborted. 


CALLING SEQUENCE 


CALL DETACH (lgu,istat ,buffer) 


lgu 
istat 


buffer 


the FORTRAN code (an integer expression, variably, or constant) 
the status return word 

istat = 0, means OK 

istat = 1, means could not get file system buffer 

null arg: get a file system buffer 

not null: use this buffer (at least 380 words) 


The following is an example of a null argument. 


CALL DETACH (lgu,istat,) 


METHOD 


See Function. 
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Attach a Temporary Mass Storage or Terminal File 


FUNCTION 


FCRA is used to create and access a temporary mass storage or terminal 
file. 


CALLING SEQUENCE: 


CALL CREATE (lgu,isize,mode,istat) 


lgu the FORTRAN file code (an integer expression, variable, or 
constant) 

the size, in words, of the temporary file wanted 

0 for a linked mass storage file 

l for a random mass storage file 

2 for a terminal file 

the status return word (see the TSS System Programmer's Reference 
Manual for TSS codes). The following codes also apply: 


isize 
mode 


istat 


successful 

1, mode is invalid 

2, file is currently open 

3, no room in AFT 

4, temporary file not available 
duplicate file name 

6, no room in PAT 

7, illegal device specified 


to 
So 
= 


tu it tu 
or 
~ 


METHOD 


If the function is successful, a FCB is created and the file code, in 
ASCII, is placed in FCB -10, -9. 
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FORTRAN DEBUGGING SYSTEM (FDS) SUBROUTINES 


Core Allocator 


FUNCTION 


FALC is a general dynamic storage allocator capable of dispensing memory 
space from a region composed of a list of disjoint free blocks. 


CALLING SEQUENCE 


The following call will allocate a block of n consecutive words: 


LDQ n,DL 
TSX1 .FDSGT 


If the request is denied, the Q register will be returned with the value 
zero. If the request is granted, bits 0-17 of the Q register will contain 
the address of the allocated block and bits 18-35 will contain the length 
of the block as designated in the call. (Currently, denial returns are 
given only in the batch mode of operation.) 


METHOD 
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The free blocks are chained together with a linked list running through 
bits 18-35 of their initial words. The last block on the chain contains 
zero in this field. Bits 0-17 of the first word of a free block contain 
the address of the last word in the block plus one. Equality of these two 
fields indicates that the block is immediately followed by another free 
block and the two may be coalesced into one. 


To locate memory space to satisfy a storage request, the chain of free 
blocks is first searched, and blocks are coalesced as the search proceeds. 
If a sufficiently large free block is found, the requested space is 
allocated and any,excess memory is returned to the free block chain. 


If no free block is found that satisfies the request, space is sought in 
the "core-hole" described by word 37 (octal) of the slave prefix. i Re 
enough memory space is available at this location, the requested space is 
allocated out of the low end and the core-hole limits are adjusted 
accordingly. 


If no memory space is available in the core-hole, a request is made to the 
operating system for additional memory. The memory space obtained in this 
operation is transformed into a new core-hole and the old core-hole is 
merged into it, if possible. 


To ensure that the core-hole does not become too small, any free blocks 


adjacent to the core-hole are merged with the core-hole and are removed 
from the free block chain. 


3=32 DD20A 


Special Entry Point 


FUNCTION 


FDBD provides a special entry point to a dummy routine. If a 
user-generated version of the FDEBUG module has been provided, control will 
be passed to the dummy routine rather than to the FDEBUG bootstrap routine. 
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FDEBUG Bootstrap 


FUNCTION 


FDEB is the bootstrap subroutine for the FDEBUG module... 


CALLING SEQUENCE 


CALL FDEBUG 
METHOD 


The FDEBUG overlay is loaded into memory (if not already loaded) and 
control is transferred to the entry point. 
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FDUMP Bootstrap 


FUNCTION 


FDUM is the bootstrap subroutine for the symbolic dump facility FDUMP. 


CALLING SEQUENCE 


CALL FDUMP 
METHOD 


The FDUMP overlay is loaded into memory (if not already loaded) and control 
is transferred to the entry point. 
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LINK/LLINK Interface 


FUNCTION 


FLKL provides an interface between LINK/LLINK and the FDEBUG module. 


CALLING SEQUENCE 


TSX1 .FDLNK 


It is assumed that the A register contains the link origin in bits 0-17 and 
the link size in bits 18-35. 


The entry point .FDLLD has the calling sequence 


CALL .FDLLD (name) 


where: 'name' is the address of the link name in character form. 


METHOD 
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Control is first passed to FDEBUG following the entry to .FDLNK so that any 
overlayed breakpoints can be removed from the breakpoint table. Control is 


then passed to .FDBDL, a _ secondary entry point in the FDEBUG bootstrap 
subroutine. 


The subroutine .FDLLD is an elaborate NOP that is capable of being 
breakpointed. The user can therefore gain control after a link has been 
loaded to install new breakpoints. 
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Delete From Wrapup List 


‘FUNCTION 


FNCL allows the user to delete subroutines from the wrapup list. 


CALLING SEQUENCE 


CALL NOCALL (subr) 


METHOD 


When a CALL NOCALL(subr) statement is included, all occurrences of the 
subroutine are deleted from the wrapup list. 
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Release Unused _ Memory 


FUNCTION 


FREL provides the capability to return a block of memory space to the free 
block list. 


CALLING SEQUENCE 


To free a block of memory, the address of the block is placed in bits 0-17 
of the Q register and the length of the block is placed in bits 18-35 (this 
is the data returned by .FDSGT). The following call is made: 


TSX1 .FDSRL 


The following call will release as much memory as possible from the high 
end of core: 


TSX1 .FDSRM 


METHOD 


Returning a block of memory space to the free list is a simple matter of 
finding the appropriate location in the chain and then adjusting the 
pointers. Blocks are not merged with the core-hole during this search. 
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Dummy Setup 


FUNCTION 


FSTP is a dummy subroutine that prevents binding of FDS modules when the 
FDS option has not been invoked, 


CALLING SEQUENCE 


TSXO .FDSET 


METHOD 


If called, this subroutine simply returns via index register zero. 
Secondary SYMDEFs are provided for FDEBUG and FDUMP that allow calls _ to 
FDEBUG and FDUMP to return to the user unless the FDS option has been 
invoked. 
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Add_to Wrapup List 


FUNCTION 


FTCL allows the user to add subroutines to the wrapup list that is 
maintained dynamically by the FORTRAN debugging system. 


CALLING SEQUENCES 


CALL ATCALL (subr) 
CALL NTCALL (subr) 


METHOD 
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The wrapup list is inspected whenever a program terminates abnormally or 
terminates when a STOP statement is executed. However, when a program 
terminates upon the execution of a CALL FTERM statement, the wrapup list is 
not inspected, thus providing a mechanism to avoid calling procedures on 
the wrapup list. 


The wrapup list is maintained in a first-in/first-out basis so that the 
subroutines entered first will be called first. If the entry is made using 
ATCALL, the subroutine will be called when the program terminates 


abnormally. If the entry is made using NTCALL, the subroutine will be 
called when the program terminates normally. 
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FF 


Timing Facility 


FUNCTION 


FTMR provides the subprogram timing measurement system of FDS, called 
FTIMER, 


CALLING SEQUENCE 


METHO 


3/77 


CALL FTIMER 


NOTE: A $ USE FTIMER control card may be included in the batch mode _ to 


force the timing measurement system to be loaded. 


D 


FTMR contains replacements for the entry points .ENTY, .RETY, and .FRETY 
that are contained in the linking subroutine FYLK. These replacements 
dynamically allocate data blocks in which timing measurement statistics are 
accumulated as the program executes. These statistics are summarized and 
printed when the program terminates (even if the termination is abnormal). 
The measurements are given only for those subprograms compiled with the FDS 
option, 
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| Wrapup and Loader 


FUNCTION 


FWRP provides the standard FORTRAN wrapup procedure as well as the overlay 
loader used by the FDEBUG and FDUMP bootstrap subroutines, 


CALLING SEQUENCE 


TSX1 .FEXIT (Termination) 
TSX1 .FDSLD (Overlay loader) 


METHOD 
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FWRP contains the standard FORTRAN wrapup procedure’ .FEXIT. When this 
procedure is called, it first determines whether the termination is normal 
or abnormal and then invokes each user-supplied wrapup routine as required. 
When all of the subroutines in the wrapup list have been called, the 
subroutine .FTERM is then invoked to provide standard termination. 


FWRP contains the overlay loader. When it is called, the A register 
contains zero (0) if an FDUMP overlay is to be loaded or two (2) if an 
FDEBUG overlay is to be loaded. 
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Linking Subroutine 


FUNCTION 


FYLK contains the subroutines .ENTY, .RETY, and .FRETY, as well as the 
setup procedure. 


CALLING SEQUENCE 


TSXO .ENTY (or .RETY or .FRETY) 
NOP eHeLie. 


METHOD 
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Subroutine .ENTY is called whenever a FORTRAN subprogram is entered. Its 
function is to save the indicator register and index register 1 in the 
error linkage area and to insert the address of the error linkage area into 
bits 0-17 of .CE.L. 


Subroutine .RETY is called whenever a RETURN statement is executed. Its 
function is to restore the indicator register and to set bits 0-17 of 
-CE.L. to the address of the caller's error linkage area. 


Subroutine . FRETY is called whenever a RETURN n statement is executed. Its 
function is to restore the indicator register, to set bits 0-17 of .CE.L. 
to the address of the caller's error linkage area, and to execute the 
return. 


The setup procedure is called by the fault processor (.FLTPR) via a TSX0 
»FDSET. Its primary function is to call FDEBUG when required prior to 
entering the main program. In the batch mode of operation, tests are 
performed to verify whether file designator 44 is present; if it is, the 
FDEBUG module is called. In the time sharing mode, bit 4 of the program 
switch word is checked to determine whether or not the FDEBUG module should 
be called. 
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SECTION IV 


MATHEMATICAL LIBRARY SUBROUTINES 


The subroutines that make up the Mathematical Library offer the FORTRAN 
user a wide range of options to help solve complex mathematical problems. These 
subroutines provide increased flexibility and computational capability to the 
FORTRAN scientific programming language. 


The purpose of the FORTRAN mathematical library subroutines is to provide a 
‘fast method of calculating the basic mathematical functions. The purpose of 
this update is to improve the accuracy and performance of the present 
subroutines and to add several new functions. The results obtained with this 
improved version of the library will in many cases differ from those obtained 
when uSing earlier versions; in these cases, answers obtained are more 
accurate. The following statements relate to the degree of accuracy for both 
single- and double-precision floating-point binary numbers. 


The single-precision subroutines were written to provide accuracy to full 
Single precision plus or minus one bit, whenever possible. This gives a 
relative error of less than 1*2**-27, 


The double-precision subroutines use polynomials with theoretical accuracy 
beyond that of the floating-point hardware (greater than 18 decimal digits). 
Any error that exists is due to rounding or truncation. Programming techniques 
were used to minimize this type of error within the routine with the result that 
most routines are accurate to about 17 or more decimal digits. 


It is assumed that the binary representation of the argument (in single or 
double precision) is the exact value, and the ensuing calculations and results 
are based on that assumption. Rounding or conversion error prior to entering 
the subroutine cannot be anticipated. However, Since argument error can cause 
serious problems in calculations, the discussion of each new or revised 
subroutine contains a brief summary of the possible effects of argument error on 
the results. This information is presented in terms of a formula that usually 
gives only the first term in the series expansion for the error term. 


The following notation is used in the error formulas: 


ACT - Correct result 

AEA - Absolute error of the argument 
AER - Absolute error of the result 
REA - Relative error of the argument 
RER - Relative error of the result 
RES - Result given by the subprogram 
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Overflow and underflow faults are not masked during the execution of the 
subroutines. In cases where the result cannot be represented as a 
floating-point number because it is too small or too large, a call is made to an 
error routine and an appropriate error message is generated. 


In certain cases, underflow may occur for valid input arguments but this 
condition will not affect the accuracy of the result. Intermediate overflow 
should not occur during the execution of the subroutines. 


MATHEMATICAL LIBRARY DESCRIPTIONS 


The FORTRAN mathematical library subroutines are composed of the function 
subprograms and are summarized in Table 4-1. The definition, calling name, and 
edit name are given for each library subroutine. The page numbers are included 
since the indexed page numbers for Section IV are no longer valid for this 
addendum, 
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Table 4-1. Mathematical Library Subroutines 
Subprogram Neiman Calling 
Function Definition Name Name No. 


1 ARCOS FASN 4-18 
cos”~ (a) 
DARCOS FDAS 4-35 
1 ACOSH FASH 4-16 
cosh”~ (a) 
DACOSH FDAH 4-34 
ARSIN FASN 4-18 
sin7! (a) 
DARSIN FDAS 4-35 
l ASINH FASH 4-16 
sinh” (a) 
DAS INH FDAH 4-34 


Absolute Value, Complex 


Arithmetic, Complex, Multiply 
and Divide 


Arccosine, Real 


Arccosine, Double 


Arccosine, Hyperbolic, Real 


Arccosine, Hyperbolic, Double 


Arcsine, Real 


Arcsine, Double 


Arcsine, Hyperbolic, 


Arcsine, Hyperbolic, Double 


Arctangent, Real 


Arctangent, Double 


Arctangent 2, Real 


Arctangent 2, Double 


Arctangent, Hyperbolic, 


e Arctangent, Hyperbolic, Double 


Cosine, Real 


Cosine, Double 


Cosine, Complex 


Cosine, Hyperbolic, Real 


cosh (a) 


Cosine, Hyperbolic, Double 


New, Subroutine 


“Revised Subroutine 
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Table 4-1 (cont). Mathematical Library Subroutines 


Subprogram Calling | Edit Page 
Function Definition Name Name No. 


Cube Root, Real CBRT FCRT 4-28 
a) 1/3 
Cube Root, Double DCBRT FDCR | 4-39 
Exponential, Real EXP FEXP 4-58 
Exponential, Double DEXP FDEX 4-40 
Exponential, Complex CE XP FCEX 4-24 
‘ EXP2 FEXP 4~58 
DEXP 2 FDEX 4-40 
EXP10 FEXP 4-58 
o* 
4-56 
4-52 
4~56 
4-52 
4-56 
4-52 
4-13 


Exponential 2, Real 


( 
2 
Exponential 2, Double 
Exponential 
1 
ee Exponential 10, Double 


e Exponential Complement, Real 


e210 


EXPC FEXC 

Exponential Complement, Double DEXPC FDXC 

EXPC2 FEXC 

DEXPC2 FDXC 

a EXPC10 FEXC 
10 -1.0 

DXPC1O FDXC 


Implied 


Exponential Complement 2, Real 


Exponential Complement 2, Double 


Exponential Complement 10, Real 


Exponential Complement 10, Double 


Exponentiation, 
Any Exponent 


Complex Base, 


Exponentiation, Complex Base, 
Integer Exponent 


Implied 


Exponentiation, Double Precision, 
Base or Exponent (or Both) 


Implied 


Exponentiation, Integer Base 
and Exponent 


Implied 


Exponentiation, Real Base 
and Exponent 


Implied 


Exponentiation, Real Base, 
Integer Exponent 


Implied 


e New Subroutine 


ee Revised Subroutine 
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Table 4-1 (cont). Mathematical Library Subroutines 


Calling] Edit Page 
Definition Name Name No. 


ALOG FALG 4-14 


Subprogram 
Function 


Logarithm, Natural, Real 


Logarithm, Natural, Double DLOG FDLG 4-42 


Logarithm, Natural, Complex CLOG FCLG 4~ 


25 
ALOG2 FALG 4-14 | 
DLOG2 FDLG 4-42 
ALOG1O | FALG 4-14 
1094 (a) 
DLOG10 | FDLG = 
POW F3XP 
aP 
DPOW FDXP 


FSIN 
sin (a) DSIN FDSN 
CSIN FCSN 
SINH F TNH 
sinh (a) 
DSINH FDPH 


2 

: TAN FTAN 4-64 
tan (a) 

DTAN FDTN 4-50 


TANH FTNH 4—66 
tanh (a) 
DTANH FDPH 4-45 


Logarithm, Base 2, Real 


Logarithm, Base 2, Double 


ee Logarithm, Common, Real 


ee Logarithm, Common, Double 


ee Power, Real (Exponentiation) 


Power, Double (Exponentiation) 


e Sine, Hyperbolic, Real 


e Sine, Hyperbolic, Double 


Square Root, Real 


Square Root, Double 


Square Root, Complex 


Real 


Tangent, 


@ Tangent, Double 


ee Tangent, Hyperbolic, Real 


e Tangent, Hyperbolic, Double 


New Subroutine 


Revised Subroutine 
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-. DEFINITIONS AND CONSIDERATIONS 


3/77 


ds 


The following descriptions show both the GMAP calling sequence and an 
example of user coding. The subroutines that do not show user coding 
are generally called by other subroutines that make use of the called 
subroutines. 


The subroutines in the math library are function subprograms and are 
not called by a CALL statement by the user. Some are called by the 
FUNCTION name (i.e., A=SQRT(X)) while others are called implicitly. 
For example, subprogram F1XP is invoked via an arithmetic statement of 
the following type: K=I**J. 


The answer is always returned in processor registers. The FORTRAN 
compiler generates the appropriate store instructions (if a 
replacement operation is involved). The registers used are determined 
by the type of the function and are as follows: 


Answer Register Processor Instruction 
Real EAQ FSTR 
Double Precision EAQ DFSTR 
Integer, Logical Q STQ 
Complex AQ STAQ 


If the answer is to be stored into a variable of a different type, the 
compiler will supply a suitable conversion algorithm. See the FORTRAN 
reference manual for the various combinations and normal type of 
results. 


When the letter Z appears in a function description, it indicates a 
complex variable, unless otherwise noted. 
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54 Errors are detected by many of the subprograms. The functional flow 
when an error is detected is as follows: 


MAIN PROGRAM 


PERFORM 
FUNCTION 


YES 


ERROR 
? 


LOAD STANDARD 
ANSWER FOR 
THIS ERROR 


CALL FXER 
TO OUTPUT 
MESSAGE 


PROCESS 
ERROR 


SET INDICATORS 
FOR RESULT 


ALTERNATE 
PROCEDURE 1 


6. Where the argument types are listed for a subprogram, the subprogram 
assumes that the arguments are of the type listed. No explicit error 
messages are issued if incorrect arguments are supplied. 


1 

It is assumed here that the alternate procedure is a FUNCTION subprogram which 
ends with a RETURN statement. This is not necessary; the alternate procedure 
could be a section of the main program (or subprogram) in which case the flow 
would not include the "SET INDICATORS---". 
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F1XP 


Exponentiation - Integer Base and Exponent 


FUNCTION 


F1XP computes I**J in an expression. 


USAGE 


K=[**J I,J are integers. 


GMAP CALLING SEQUENCE 
CALL .FXP1(I,J) 


: METHOD 
+ For I=0,J=0, chen K=0, error code 1 
2% For I=0, J<0, then K=(2**35)-2, error code 2 
3. For I=0, J>0, then K=0 
4. For I40, J=0, then K=1 
54 For I=l1, J#40, then K=1l 
6. For I=-l, J is even, K=1 
1 For I=-l, J is odd, K=-1 
8. For |Il122, J<0, K=0 
9. For 1Il22, J<36, K=I**J 
10. For II! 22, J236 
a oe If (KI>(2**35)-1, error code 60 
If J is even, K=(2**35)-2 
If J is odd, I is negative, K=-(2**35)-2 


If J is odd, I is positive, K=(2**35)-2 
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F1XP 


DD20A 


F1XP F1XP 


RETURNS 


Normal return is to the next executable statement in the calling program, 
when AQ(0-71) contains K. Normally only the last 36 bits of Q (36-71) are 
used, 

die Error code 1 if I=0, J=0, then K=0. Execution continues. 

2 Error code 2 if I=0, J<0, then K=(2**35)-2. Execution continues. 


as Error code 60 if K=I**J and [KI >(2**35)-l. Execution continues. 
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F2XP F2XP 


Exponentiation - Real Base, Integer Exponent 


FUNCTION 
F2XP computes A**K in an expression. 
USAGE 


B = A**K A is real or double precision, K is an integer. 


- GMAP CALLING SEQUENCE 


CALL .FXP2(A,K) for real A 


CALL .FDXP1(A,K) for double precision A 


Ls For A=0, K=0, then B=0, error code 3 

2% For A=0, K>0, then B=0 

o< For A=0, K<0O, then B=10**38, error code 4 
4. For A#40, K=0, then B=1.0 


5. For negative values of K, change sign of K and take the reciprocal of 
the result. 


6. K is an integer with values from -2**35 to (2**35)-1. 


A and A**K are floating-point numbers with values from -2**127 to 
(2**127) -2**64, 


de A**K is accurate to eight decimal positions for FXP2 or 16 decimal 
positions for FDXPl. 


RETURNS 
Normal return is to the next executable statement of the calling program. 


1. Error code 3 if A=0 and K=0. Then A**K=0.0. Execution continues. 


2. Error code 4 if A=0 and K<0O. Then A**K=10**38. Execution continues. 
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F3XP 


F3XP 


Exponentiation - Real Base and Exponent 


FUNCTION 


F3XP computes X**Y in an expression. 


USAGE 


C=X**Y X and Y are real. 


GMAP CALLING SEQUENCE 


CALL .FXP3 (X,Y) 


CALL POW(X,Y) 


METHOD 


The subroutine is based on the formula 
X**Y = 2** (Y*ALOG2 (X)). 

The subroutine uses the entry point .ALOG2 into the ALOG2 subroutine 

and the entry point .EXP2 into the EXP2 subroutine. (EXP2(Z) finds 

the value of 2**Z.) 


The largest known relative error is less than 1*2**-27 (one bit in 
single precision). 


EFFECT OF ARGUMENT ERROR 


RETURNS 


3/77 


The error propagation can be calculated from the guidelines provided in the 
ALOG2 and EXP2 subroutines. 


Normal return is to the next executable statement in the calling program, 
except as specified below: 


1. 


If (X.LT.0.0), the result is set to 0.0 and FXEM is called with error 
code 5. 


If (X.EQ.0.0 .AND. Y.EQ.0.0), the result is set to the largest 
possible number and FXEM is called with error code 7. 
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F3XP 


F3XP 


Appropriate overflow and underflow errors may occur in EXP2 (edit name 


3. 
FEXP) if the results cannot be expressed as a single-precision number. 
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F4XP . F4xXP 


Exponentiation - Complex Base, Any Exponent 


FUNCTION 


F4XP computes A**B in an expression. 


USAGE 


C=A**B A and B are complex. 
C=A**B A is complex; B is real or double precision. 


GMAP CALLING SEQUENCE 


CALL .FCXP2 (A,B) for complex base and exponent 
CALL .FCXP3 (A,B) for complex base and real or double-precision exponent 
METHOD 


(CA, CB are the arguments) 


CALL CLOG(CA) Complex logarithm of complex A 
STAQ AA CLOG of argument 1 

CALL .FCFMP(AA,CB) Product 

STAQ AA 

CALL CEXP (AA) Complex exponential 


For complex base and real or double-precision exponent, a complex number is 
formed with the real part equal to the argument and the imaginary part 
equal to zero. Otherwise, it uses the same method as complex base and 
exponent. 


Accuracy: 7 places + 1.0*10**-7 


RETURNS 


Normal return is to the next executable statement in the calling program. 
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FALG FALG 


Natural, Base 2, and Common Logarithms 


FUNCTION 


FALG computes logg(a) for ALOG(X), where X.GT. 0.0; log, (a) for ALOG2(X), 
the logarithm base 2; and log 49 (a) for ALOG10(X), the common logarithm or 
logarithm base 10. 


USAGE 


A=ALOG (X) 
A=ALOG2 (X) X is real. 
A=ALOG10 (X) 


GMAP CALLING SEQUENCE 


CALL ALOG(X) for loge (a) 
CALL ALOG2(X) for log» (a) 
CALL ALOG10(X) for 10949 (a) 


As The basic function calculated by the subroutine is ALOG2(X). The 
results for ALOG and ALOG1O are obtained by using the expression 
ALOG(X) = ALOG2(X) *ALOG(2.0) or ALOG10(X) = ALOG2(X) *ALOG10 (2.0). 

as The constants ALOG(2.0) and ALOG10(2.0) are stored in the subroutine; 
the necessary multiplication is performed just prior to the return. 
No multiplication is performed for ALOG2. 


cr The polynomial approximation for ALOG2 is accurate over the range 
(.707, 1.414) and is of the form 


POO+Z**2* (PO1+2Z**2* (PO2+Z**2*PO03) ) 
where: Z = (RRX-1.)/(RRX+1.). 


The range reduction of the argument utilizes the fact that: 


ALOG2 (RRX*2**IEXP) = ALOG2(RRX) + IEXP. 
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FALG FALG 


4. The following steps are performed: 
a. If X is already within the range (.707, 1.414), omit the 
remaining steps and use the polynomial directly. This prevents 
loss of precision for arguments near 1.0. 


b. EXP = exponent of X. 


om RRX = (X with exponent set to 0), so RRX is within the range (.5, 
1.0). 
d. Since the polynomial is accurate over the range (.707, 1.414), 


multiply RRX by the SQRT(2) and compensate for the multiplication 
by subtracting 0.5 (which is ALOG2(SQRT(2))) from the final 
answer. The actual multiplication is not performed; instead, the 
quotient (X-1/SQRT(2))/(X+1/SQRT(2)) is formed in finding 2 for 
use in the polynomial. Note that (SQRT(2)*X-1)/(SQRT(2)*X+1) is 
equivalent to the previous quotient. 


e. In summary: ALOG2(X) = exponent of X - 0.5 + ALOG(SQRT(2) *RRX) 
where: RRX is X with the exponent set to 0.0. 
oe The largest known relative error is less than 1*2**-27 (one bit in 
single precision). 
EFFECT OF ARGUMENT ERROR 
AER = REA, which means that the RER for arguments near 1.0 may be large. 


RETURNS 


Normal return is to the next executable statement of the calling program, 
except as specified below: 


de If (X.EQ. 0.0), the result is set to the most negative number and FXEM 
is called with error code 9. 


Zs If (X.LT.0.0), the result is set to 0.0 and FXEM is called with error 
code 10. 
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FASH FASH 


Real Hyperbolic Arcsine, Arccosine, and Arctangent 


FUNCTION 


FASH computes sinh"! (a) for ASINH(X), the value of the hyperbolic arcsine 
of xX for all inputs; cosh"! (a) for ACOSH(X), the value of the hyperbolic 
arccosine; and tanh"! (a), the value of the hyperbolic arctangent. 


USAGE 


A=ASINH (X) 
A=ACOSH (X) X is real. 


=ATANH (X) 


GMAP CALLING SEQUENCE 


CALL ASINH(X) for sinh-! (a) 
CALL ACOSH(X) for cosh“! (a) 
CALL ATANH(X) for tanh 7-1 (a) 


METHOD 


1. ASINH is calculated as ALOG(X+SORT(X**24+1.0)). 
2. ACOSH is calculated as ALOG(X+SQRT(X**2-1.0)). 


For the above two expressions, SQRT (X**2-1.) is set equal to xX for 
X.GE.1.0E5. 


36 ATANH is calculated as .5*ALOG((1.+X)/(1.-X)). 
4. These subroutines call other routines that are accurate to single 


precision. An attempt was made to minimize the accumulated error in 
calculating the mathematically correct expressions. 
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FASH FASH 


RETURNS 
Normal return is to the next executable statement of the calling program, 
except as specified below: 


i If (X.LT.1.0) for ACOSH, the result is set to 0.0 and FXEM is called 
with error code 91. 


De If (X.LE.-1.0) for ATANH, the result is set to the most negative 
number and FXEM is called with error code 92. 


3 If (X.GE.1.0) for ATANH, the result is set to the most positive number 
and FXEM is called with error code 92. 
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FASN : FASN 


Real Arcsine and Arccosine 
FUNCTION 


FASN computes sin~! (a) for ARSIN(X), the principal value of the 
trigonometric arcsine of X (in radians); and cos“! (a) for ARCOS(X), the 
principal value of the trigonometric arccosine of X (in radians). 


USAGE 
A=ARSIN (X) 
A=ARCOS (X) 
X represents any number within the range (-1.0, 1.0), meaning -1.05xX $1.0. 


The arcsine and arccosine are undefined for arguments outside of this 
range. 


GMAP CALLING SEQUENCE 


CALL ARSIN(X) for sin “! (a) 
CALL ARCOS(X) for cos 7! (a) 


a Initial range reduction is performed by setting AX=ABS(X), which 
reduces the range to (0, ARSIN(PI/2)). 


2% The range between (0, ARSIN(PI/2)) is divided into four parts, with 
the partition points of R expressed as follows: 


R(O) = 0. = SIN(0) 
R(1) = .5 = SIN(PI/6) 
R(2) = .866 = SIN(PI/3) 
R(3) = .965 = SIN(5*PI/12) 
R(4) = 1.0 = SIN(PI/2) 
3s The basic formulas used to derive the actual formulas used in the 


subroutine for the above four ranges are: 


a. T (X) = COS(N*ARCOS (X)), which means: 


T1(X) = X 
T2(X) = 2*X**2—] 
T4(X) = 8*¥X**4~8*X** 24+] 


This is the standard Chebyshev polynomial of the first kind. 


b. ARSIN(X) = (N-1)*PI/(2*N) + (1/N) *ARSIN (TN (X) ) 


where: SIN((N-2) *PI/(2*N)) <= X $1.0. 
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FASN FASN 


This form is most efficient near: 
SIN (PI/2-PI/(2*N)), where TN(X) approaches 0. 
Cs ARSIN(X) = PI/2 - 2*ARSIN(SQRT((1-X)/2)). 


d. ARCOS (X) 


PI/2 - ARSIN(X) for 0 © ARCOS(X) < PI. 


4. The following chart presents the applicable formulas and the values of 
N for the different ranges. 


Transformation Transform 

Interval Used Range 
SIN (A) < X $s SIN(B) Formula N 7RRX/ < SIN(Z) 

A B Z 

0 PI/6 None PI/6 
PI/6 PI/3 (3a) 2 PI/6 
PI/3 5*PI/12 (3b) 4 PI/6 
5*PI/12 PI/2 (3c) PI/24 


Two different minimax polynomials are used to approximate ARSIN(RRX) 
over the ranges (0, PI/6) and (0, PI/24). For the range (0, PI/6), a 
polynomial of the form 


Y*P (Y**2) /Q(Y**2) 


is used, where P is of the second order and Q is of the third order. 
For the range (0, PI/24), P and Q are both of the first order. 


oe Using all of the above formulas, the following formulas actually used 
in the subroutine can be derived. 


Range 1: 


for X (0.0, 0.5); equivalently 
X (SIN(0.0), SIN(PI/6)), or 
X (COS(PI/2), COS(5*PI/6)); 


POARG = X 

ARSIN = 0.0 + PQARG*1.0*P (PQARG** 2) /Q (PQARG**2) 

ARCOS = PI/2. - PQARG*1.0*P (PQARG** 2) /Q (PQARG** 2) 
Range 2: 


for X (0.5, .866025404); equivalently 
X (SIN(PI/6), SIN(PI/3)), or 
X (COS(5*PI/6), COS(2*PI/3)); 


PQARG = 2*xX**2-1 
ARSIN = PI/4. + PQARG*0.5*P (PQARG**2) /Q (PQARG** 2) 
ARCOS = PI/4. - PQARG*0.5*P (PQARG**2) /Q (PQARG**2) 
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FASN FASN 


Range 3: 
for X (.866025404, .965925826); equivalently 
or 


X (SIN(PI/3), SIN(5*PI/12)) 
D4 (COS (2*PI/3), COS (PI/12) ) 


v 
° 
‘ 


PQARG = 8*X**4 — 8*X**2 + 1. 
ARSIN = 3*PI/8 + PQARG*0.25*P (PQARG** 2) /Q (PQARG** 2) 
ARCOS = PI/8 - PQARG*0.25*P (PQARG**2) /Q (PQARG** 2) 


Range 4 (uses a different polynomial approximation) : 


for X (.965925826, 1.0); equivalently 
X (SIN(5*PI/12), SIN(PI/2)), or 
X (COS(PI/12), COS(0.0)); 


PQARG = (1-X)/2. 
ARSIN = PI/2 -— SQRT(PQARG) *2*P (PQARG**2) /Q (PQARG** 2) 
ARCOS = 0.0 + SQRT(PQARG) *2*P (PQARG**2) /Q (PQARG** 2) 


Now, for all ranges 
If the original argument was negative: 


ARSIN 
ARCOS 


~ARSIN 
PI - ARCOS 


Wo 


Then, return. 
For ranges 1, 2, and 3, the same formula for the polynomial is used. 
In range 4, a different formula is used and a call to SQRT must be 
made. 


6. The largest known relative error is less than 1*2**-27 (one bit in 
single precision). 


EFFECT OF ARGUMENT ERROR 


AER = AEA/(1-X**2) **.5 


If the value of X is small, RER = AEA for the arcsine. 


RETURNS 
Normal return is to the next executable statement of the calling progxam, 
except as specified below: 


1. Tf (X.LT.-1.0), FXEM is called with error code 83. Then the answer 
is found for X = -1.0. 


Zs If (X.GT.1.0), FXEM is called with error code 83. Then the answer is 
found for X = +1.0. 
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FATN FATN 


Real Arctangent 


FUNCTION 


FATN computes ‘the principal value (in radians) of the arctangent of Z for 
ATAN(Z) or the arctangent of (Y,X) for ATAN2 (Y,X), where Z = Y/X, for any 
valid input argument(s). 


USAGE 


A=ATAN (2) 
A=ATAN2 (Y,X) 


(Y,X) or 2% are any numbers, except that the expression (Y,X)=(0.,0.) is 
invalid. 


CALLING SEQUENCE 


CALL ATAN(Z) for tan-l (a) 
CALL ATAN2(Y,X) for tan-1l (a/b) 


METHOD - 


Le The algorithm is the same for both ATAN and ATAN2, except that for 
ATAN2(Z) the Z is calculated as Z = Y/X. For ATAN(Z), the answer may 
be located in either quadrant 1 or quadrant 4 (-PI/2 <ATAN < PI/2). 
For ATAN2(Y,X), the answer will be in the correct quadrant with (-PI < 
ATAN2 S PI). 


ae The ABS(Z) is found, thus reducing the range to (0.0, arctangent 
(PI/2)). 


3. The range between (0, arctangent (PI/2)) is then divided into nine 
parts with partition points R and evaluation nodes S. The partition 
points are chosen as follows: 


R(0) =0 

R(I) =tangent ((2*I-1) *PI/32.) I=1,2,...8 
R(9) =infinity 

S (I) =tangent ((2*I-2) *PI/32.) 132394052 


4. The range of Z, and thus the value of the index I, is then found 
within the range (R(I-1), R(I)). That is, R(I-1)<Z <R(I). I is found 
using a binary search to determine the range of Z. After I is found, 
the formula 

ATAN(Z) = ATAN(S(1I)) + ATAN(T) 


is used, where: T=S(I) **-1-(S(I) **-2+1)/(S (I) **-14+2Z). 
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FATN 


FATN 


ATAN(S(I)) and the values for (S(I)**-1) and (S(1I)**-2+1) are stored 
in tables. 


T is within the range (-PI/32., PI/32.) and ATAN(T) is calculated 
using a polynomial approximation of the third order which has the form 


POO+T* (PO1+T* (P02+T*P03) ) 
and provides accuracy beyond single precision. 
After the value of ATAN is found for positive 2Z (quadrant 1), the 
answer is placed in the proper quadrant by subtracting it from PI 


and/or negating it if necessary. 


The largest known relative error is less than 1*2**-27 (one bit in 
single precision). 


EFFECT OF ARGUMENT ERROR 


AER=AEA/ (1-2**2) for a lesser value of Z. 
RER=REA for the larger values of Z. 


RETURNS 


Normal return is to the next executable statement of the calling program, 
except as specified below: 


1. 


3/77 


Tf (Y.EQ.0.0.AND. X.EQ.0.0), the result is set to 0.0 and FXEM is 
called with error code ll. 
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FCAB | FCAB 


Complex Absolute Value 


FUNCTION 


FCAB computes the absolute value of the complex number X in the expression 
CABS (X). 


USAGE 


A=CABS (X) 


X represents any complex number of the form (A,B). 


GMAP CALLING SEQUENCE 


CALL CABS (X) 


METHOD 


Ls The complex argument X is treated as a number pair (A,B): 


SML = MIN(A,B) 
BIG = MAX(A,B) 
Qs CABS = ABS (BIG*SQRT(1.+(SML/BIG) **2). Note that the answer is a real 
number. 
3. The largest known relative error is less than 1*2**-27 (one bit in 


single precision). 


EFFECT OF ARGUMENT ERROR 


RER = (1./2.)*REA 


RETURNS 


Normal return is to the next executable statement of the calling program. 
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FCEX FCEX 


Complex Exponential 


FUNCTION 


FCEX computes e**Z for CEXP(Z) in an expression. 


USAGE 


A=CEXP (Z) Z is complex. 


GMAP CALLING SEQUENCE 


CALL CEXP (Z) 


METHOD 
te (Sy) 
ne ae! 
2 
B= e**X 
X = B*SIN A 
Y = B*SIN Y 
ANSWER = (X,Y) 
26 Z and e**X are complex numbers(X,Y), with 
XS88.028, !¥I<2**27, and | Yt+7/<2**27, 
2 
ee Each part (X,Y) of e**Z is accurate to seven decimal positions. 
RETURNS 


Normal return is to the next executable statement of the calling program, 
except as specified below: 


1. Error code 26 if X>88.028. Then e**Z=(10**38, 10**38). Execution 
continues. 


25 Error code 27 if IYi22**27 or if Lee / el Sots i. Then e**Z=(0.0). 
Execution continues. 
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FCLG FCLG 


Complex Natural Logarithm 


FUNCTION 


FCLG computes loggZ for CLOG(Z) in an expression. 


USAGE 


=CLOG(Z) Z is complex. 


GMAP CALLING SEQUENCE 


CALL CLOG(Z) 


METHOD 
i log e@Z=1oge (X,Y) (where Z=(X,Y) ) 
=(logeIlZl, arctan Y/X) 
2. Z and loggZ are complex numbers; values of X and Y range from 
-—2**127 to *2**127-2**100 inclusive. 
ce logeZ is accurate to seven decimal positions. 
RETURNS 


Normal return is to the next executable statement of the calling program, 
except as specified below: 


cle Error code 28 if Z2=(0,0). Then loggZ =(-(10**38), 0.0). 
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FCMP FCMP 


Complex Multiplication and Division 


FUNCTION 


FCMP computes (A,B)*(C,D) or (A,B)/(C,D) in an expression. 


USAGE 


A=B*C or A=B/C B and C are complex. 


GMAP CALLING SEQUENCE 


CALL .FCFMP(R,S) for R*S 


CALL .FCFDP(R,S) for R/S 


R = (A,B) 
S = (C,D) 
METHOD 
l. (A,B) *(C,D) = (A*C-B*D,A*D+B*C) 
2. (A,B) /(C,D) =((A,B) *(C,-D) ) / (C#*24D**2) = (A*C+B*D,B*C~A*D) /(C**24+D**2) 


3. If (A,B) = (0,0), then the quotient = (0,0). Otherwise, 


U = D/C 

¥ = (A/C) /(1+U**2) 
T = B/A 

X = Y*(T*U+1) 

Y = Y* (T-U) 


4, Before computing (A,B)/(C,D), replace the numerator by (-B,A) if IAI 
1Bl , and the denominator by (-D,C) if ICI < IDI. Adjust the quotien 
(X,Y) accordingly: 


= 
£ 


ae If |AI>IBI and ICI>IDI, then the result = (X,Y). 


b. If IAIS!IBI and ICISIDI, then the result = (X,Y). 
c. If IBIZIAl and ICI>{DI, then the result = (Y,-X). 
d. If JAI>IBI and IDIZIC!, then the result = (-Y,X). 
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FCMP FCMP 


Bie A, B, C, D, X, and Y are real numbers, with values from -2**127 to 
(2**127)-2**100 inclusive. 


6. The answer is accurate to eight decimal positions. 
RETURNS 


Normal return is to the next executable statement of the calling program, 
except as specified below: 


Ls In a division operation, error code 25 if (C,D)=(0,0). Then 
(A,B) /(C,D) =(10**38,10**38). Execution continues. 
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Real Cube Root 
FUNCTION 
FCRT computes the cube root of X in the expression CBRT(X). 
USAGE 
A=CBRT (X) 
X represents any number. 
GMAP CALLING SEQUENCE 
CALL CBRT (xX) 
METHOD 
dy X is expressed as RM*2**IE. RM is real. IE is integer. Range 
reduction is performed using the formula 
CBRT (RM*2**TE) = RM/3*CBRT(2**IE). 
If I=IE/3 and J=IE-3*1I, then 
CBRT (RM*2**IE) = CBRT (RM*2¥**J) *2**I, 
2% RRX is RM*2**J and is within the range (.5,4.). The CBRT(RRX) is 
found using an initial polynomial approximation of the form 
YO = P00+RRX* (P01+RRX*P02) 
which makes YO accurate to approximately 1.8 decimal digits. 
3. Three Newton's iterations are then used with the formula 
Yl = YO-(YO-X/Y0**2) /3 
to extend the accuracy to full single precision. 
4. The largest known relative error is less than 1*2**-27 (one bit in 
single precision). 
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EFFECT OF ARGUMENT ERROR 


RER = (1./3.) *REA-(1./9.) *REA**2... 


RETURNS 


Normal return is to the next executable statement of the calling program. 
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Complex Sine and Cosine 


FUNCTION 


FCSN computes sin Z or cos Z for CSIN(Z) or CCOS(Z) in an expression, where 
Z is in radians, 


USAGE 


A=CSIN(Z) or A=CCOS(Z) Z is complex. 


GMAP CALLING SEQUENCE 


CALL CSIN(Z) for sin Z 


CALL CCOS(Z) for cos 2 


5 Sin Z=sin (X,Y) (where Z=(X,Y)) 
=Sin X*cos(0,Y)+cos X*sin(0,Y) 
=(sin X*cosh Y,0)+(0,cos X*sinh Y) 
=(sin X*cosh Y, cos X*sinh Y) 


26 cos Z=sin (Z+7) 


3% Z, Sin Z, and cos Z are complex numbers, with [| XI<2**27, 


IX+a71<2**27, and IYI<88.028. 
2 


4. The answer is accurate to seven decimal positions. 
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RETURNS 


Normal return is to the next executable statement of the calling program, 
except as specified below: 
ve Error code 29 if: IX!I>2**27 
IX+ /21>2**27 
e**¥=0 
Then the answer is (0,0). Execution continues. 


2x Error code 30 if !1Y!>88.028. Then the answer is (10**38, 10**38). 
Execution continues. 
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Complex Square Root 


FUNCTION 

FCSQ computes the square root of Z in the expression CSORT(Z). 
USAGE 

A=CSQRT(Z) 

Z is any complex number of the form (A,B). 


GMAP CALLING SEQUENCE 


CALL 


METHOD 


RETURNS 


CSQRT(Z) 


The argument is expressed as A+B*i. 

ANSA = SQRT(/A/+CABS ((A,B))/2.). 

If (A.LT.0.AND.B.LT.0), ANSA = ~-ANSA, 

ANSB = B/(2.*ANSA). 

If (A.LT.0), return answer (ANSB, ANSA); else return (ANSA, ANSB). 


These formulas can be derived from the definitions for complex number 
multiplication. 


The largest known relative error is less than 1*2**-27 (one bit in 
single precision). 


ARGUMENT ERROR 


Using i = (-1.)**.5, express A+B*i as Rl*e** (A1*i). 
Express the answer as R2*e**(A2*i). 


RER (R2) 
RER (A2) 


-5* REA (R1) 
REA (Al) 


Normal return is to the next executable statement of the calling program. 
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Exponentiation - Complex Base, Integer Exponent 


FUNCTION 


FCXP computes A**K in an expression. 


‘USAGE 


C=A* *K A is complex. K is an integer. 


GMAP CALLING SEQUENCE 


CALL .FCXP1(A,K) 


METHOD 


dss If A=(0,0), K=0 then A**K =(0,0) error code 14 
ae If A=(0,0), K>0 then A**K =(0,0) 
Se If A=(0,0), K<0 then A**K =(10**38,0) error code 15 


4. A and C are complex numbers (X,Y) with values X and Y from 
—-2**127 to (2**127)-2**100. 


K is an integer with values from -2**35 to (2**35) <1. 
5. Each C(X,Y) is accurate to eight decimal positions. 
RETURNS 


Normal return is to the next executable statement in the calling 


bis Error code 14 if A=(0,0) and K=0. Then A**K=(0,0). 


2% Error code 15 if A=(0,0) and K<0. Then A**K=(10**38,0.0). 
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FDAH | FDAH 


Double Precision Hyperbolic Arcsine, Arccosine, and Arctangent 


FUNCTION 


FDAH computes sinh™! (a) for DASINH(X), the value of the hyperbolic arcsine 
of X for all inputs; cosh™! (a) for DACOSH(X), the value of the hyperbolic 
arccosine; and tanhwl (a) for DATANH(X), the value of the hyperbolic 
arctangent. 


USAGE 


A=DASINH (X) 
A=DACOSH (X) X is double precision. 
A=DATANH (X) 


GMAP CALLING SEQUENCE 


CALL DASINH(X) for sinh“! (a) 
. CALL DACOSH(X) for cosh“! (a) 
CALL DATANH(X) for tanh! (a) 


METHOD 


Dy DASINH is calculated as DALOG(X+DSORT(X**2+1.0)). 
2.  DACOSH is calculated as DALOG (X+DSORT(X**2-1.0)). 


For the above two expressions, DSQRT(X**2-1.) is set equal to 
X for X.GE.1.0E10. 


3. DATANH is calculated as .5*DALOG((1.+X) /(1.-X)). 
4. These subroutines call other routines that are accurate to double 


precision. An attempt was made to minimize the accumulated error in 
calculating the mathematically correct expressions. 


RETURNS 
Normal return is to the next executable statement oF the calling program, 
except as specified below: 


is If (X.LT.1.0) for DACOSH, the result is set to 0.0 and FXEM is’ called 
with error code 9l. 


2s If (X.LE.-1.0) for DATANH, the result is set to the most negative 
number and FXEM is called with error code 92. 


3s If (X.GE.1.0) for DATANH, the result is set to the most positive 
number and FXEM is called with error code 92. 
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FDAS 


Double Precision ArcsSine and Arccosine 


FUNCTION 


USAGE 


GMAP 


ME THO 


FDAS 


computes sin7l (a) for DARSIN(X), the value of the trigonometric 


arcsine of X (in radians); and cos~! (a) for DARCOS(X), the value of the 
trigonometric arccosine of X (in radians). 


A=DARSIN (X) 
A=DARCOS (X) 


X represents any number within the range (-1.0, 1.0), meaning -1.0 © X< 


1.0. 


The arcsine and arccosine are undefined for arguments outside of this 


range. 


CALLING SEQUENCE 


CALL 


DARSIN(X) for sin~! (a) 


CALL DARCOS(X) for cos 7! (a) 


D 


The method used for DARSIN is the same as that used for the 
Single-precision ARSIN subroutine (edit name FASN), except that more 
accurate minimax polynomials are used for the approximation to DARSIN. 


The polynomial used in ranges 1, 2, and 3 is of the order 6 and _ the 
polynomial used in range 4 has P of the order 3 and Q of the order 4. 


The minimax polynomial is accurate to more than 19 places. Any error 
in the result is caused only by rounding and truncation; the result 
should be accurate to approximately 17 places. 


EFFECT OF ARGUMENT ERROR 


AER=AEA/ (1-X**2) **.5 


3/77 


If the value of X is small, RER = AEA for the arcsine. 
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-RETURNS 
Normal return is to the next executable statement of the calling program, 
except as specified below: 


P If (X.LT.-1.0), FXEM is called with error code 83. Then the answer is 
found for X=-1.0. 


vam If (X.GT.1.0), FXEM is called with error code 83. Then the answer is 
found for X=+1.0. 
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FDAT FDAT 


Double Precision Arctangent 


FUNCTION 


FDAT computes the principal value (in radians) of the arctangent of Z for 
DATAN(Z) or the arctangent of (Y,X) for DATAN2(Y,X), where Z=Y/X for any 
valid input argument(s). 


USAGE 


A=DATAN (2) 
A=DATAN2 (Y,X) 


(Y,X) or Z are any numbers, except that the expression (Y,X) = (0.,0.) is” 
invalid. 


GMAP CALLING SEQUENCE 


CALL DATAN(Z) for tan 7! (a) 
CALL DATAN2(Y,X) for tan-! (a/b) 


METHOD 


1 The algorithm is the same for both DATAN and DATAN2, except that for 
DATAN2(Z) the Z is calculated as Z=Y/X. For DATAN(Z), the answer may 
be located in either quadrant 1 or quadrant 4 (-PI/2 < DATAN S$ PI/2). 
For DATAN2(Y,X), the answer will be in the correct quadrant with (-PI 
< DATAN2 S PTI). 


2 The method used for DATAN and DATAN2 is the same as that used for the 
single-precision subroutines ATAN and ATAN2 (edit name FATN), except 
that a minimax polynomial of the order seven (instead of the order 
three) is used for approximating the DATAN function. Using the same 
notation described for the ATAN function, the polynomial is 


POO+1* (PO1L+T* (P02+T* (P03+T* (PO04+T* (P05+T* (P06+T*P07)))))) 
which provides accuracy suitable for double precision. 
cr The minimax polynomial is accurate to more than 19 places. Any error 


in the result is caused only by rounding and truncation; the result 
should be accurate to approximately 17 places. 
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EFFECT OF ARGUMENT ERROR 


‘AER=AEA/ (1-2**2) for a lesser value of Z. 
RER=REA for the larger values of 2Z. 


RETURNS 


Normal return is to the next executable statement of the calling program, 
except as specified below: 


a ar If (Y.EQ.0.0.AND.X.EQ.0.0), the result is set to 0.0 and FXEM is 
called with error code 24. 


3/77 4-38 DD20A 


FDCR FDCR 


Double Precision Cube Root 


FUNCTION 


FDCR computes the cube root of X in the expression DCBRT(X). 


USAGE 


A=DCBRT (X) 


X represents any double-precision number. 


GMAP CALLING SEQUENCE 


CALL DCBRT (X) 


METHOD 


a 
Ls The method used for DCBRT is the same as that used for the 
Single-precision CBRT subroutine (edit name FCRT), except that one 
more Newton's iteration is performed (for a total of four) to provide 
accuracy to full double precision. 


EFFECT OF ARGUMENT ERROR 


RER = (1./3.)*REA~-(1./9.) *REA**2... 


RETURNS 


Normal return is to the next executable statement of the calling program. 


3/77 4-39 DD20A 


FDEX . FDEX 


Double Precision Exponential, Base 2, and Base 10 


FUNCTION 


FDEX computes e? for DEXP(X), the value of the exponential of xX for all 
valid input argument(s); 22 for DEXP2(X), the value of 2.**x: and 108 for 
DEXP10(X), the value of 10.**xX. 


USAGE 


A=DEXP (X) 
A=DEXP2 (X) 
A=DEXP10 (X) 


X represents any double-precision value so that the result will be 
presented as a double-precision number. 


GMAP CALLING SEQUENCE 


CALL DEXP (X) 
CALL DEXP2 (X) 
CALL DEXP10 (xX) 


METHOD 


Ls The basic function calculated by the subroutine is DEXP2(X). The 
results for DEXP and DEXP10 are calculated considering that 


DEXP(X) = DEXP2(X*log2(e)) and 
DEXP10(X) = DEXP2(X*log2(10)). 


The constants log2(e) and log2(10) are stored within the subroutine, 
2s Range reduction is performed considering that 
DEXP2 (RRX+INTX) = DEXP2 (RRX) *DEXP2 (INTX). 


INTX is the integer value nearest to X and RRX is xX minus’ INTX. 
Therefore, RRX is within the range (-.5,.5). 


3.  DEXP2(RRX) is found using a minimax polynomial approximation of the 
form 


(O(KP* 2) FR*P (X** 2) 7 (0 (A882 eX P (XRF 2)) 


where the polynomial P is of the second order and Q is of the third 
order. This provides results that are suitable for double precision. 
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FDEX 


FDEX 


DEXP 2 (RRX) *DEXP2(INTX) is found by adding INTX to the exponent of the 
result for DEXP2(RRX). 


The minimax polynomial is accurate to more than 19 places. Any error 
in the result is caused only by rounding and truncation; the result 
should be accurate to approximately 17 places. 


EFFECT OF ARGUMENT ERROR 


RETURNS 


3/77 


RER 
AER 


NOTE : 


ttl 


AEA 
REA* X 


Serious error amplification can occur if small errors are 
encountered in a large value of X. 


Normal return is to the next executable statement of the calling program, 
except as specified below: 


1. 


If (result.GT.2.**127), the result is set to the largest number and 
FXEM is called with error code 19. 


If (result.LT.2.**-128), the result is set to 0.0 and FXEM is called 
with error code 89. 
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Double Precision Natural, Base 2, and Common Logarithms 


FUNCTION 


FDLG computes log,(a) for DLOG(X), the natural logarithm of any input 
argument X where X.GT.0.0; Log, (a) for DLOG2(X), the logarithm base 2; and 
log 49 (a) for DLOG10(X), the common logarithm or logarithm base 10. 


USAGE 


A=DLOG (X) 
=DLOG2 (X) 
A=DLOG1O (X) 


X represents any double-precision number .GT.0.0. 


GMAP CALLING SEQUENCE 


CALL DLOG(X) for log, (a) 
CALL DLOG2(X) for log, (a) 
CALL DLOG10(X) for 109 19 (a) 


METHOD 


Le The method used in this subroutine is exactly the same as that used 
for the single-precision subroutine ALOG, ALOG2, and ALOG10 (edit name 
FALG), except that a more accurate polynomial approximation is used. 
The polynomial used for approximating DLOG2 is of the form 

Z*P(Z**2)/Q(Z**2) 
where Z is the same as that found in the single-precision subroutine. 
2. P is a third-order polynomial of the form 
POO+Z**2* (PO1+Z**2* (PO2+Z**2*POQ3) ) 
and Q is a fourth-order polynomial of the form 
QO00+Z** 2% (QO014+Z** 2* (Q024+2%**2* (Q034+2**2*Q04) ) ) 
which provides results suitable for double precision. 
36 The minimax polynomial is accurate to more than 19 places. Any error 


in the result is caused only by rounding and truncation; the result 
should be accurate to approximately 17 places. 
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EFFECT OF ARGUMENT ERROR 


AER = REA, which means that the RER for arguments near 1.0 may be large. 


RETURNS 


Normal return is to the next executable statement of the calling program, 
except as specified below: 


lL. If (X.EQ.0.0), the result is set to the most negative number and FXEM 
is called with error code 20. 


2% If (X.LT.0.0), the result is set to 0.0 and FXEM is called with error 
code 21. 
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Double Precision Remaindering 


FUNCTION 


FDMD computes A=X(mod Y) for DMOD(X,Y) in an expression. 


USAGE 


A=DMOD (X,Y) X and Y are double precision. 


GMAP CALLING SEQUENCE 


CALL DMOD (X,Y) 


METHOD 
Le If Y=0, then A=X. Otherwise, compute Z = the greatest integer < X/Y 
and give Z the same sign as that of X/Y. Then A=X- Y*Z. 


26 A, X, and Y are double=-precision numbers, with values from -2**127 to 
(2**127) -2**64 inclusive. 


3% A is accurate to 63 binary positions. 


NOTE: If X>yY*10**10, the results may be incorrect; A may be greater 
than Y. 


RETURNS 


Normal return is to the next executable statement in the calling program. 
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Double Precision Hyperbolic Sine, Cosine, and Tangent 


FUNCTION 


FDPH computes sinh(a) for DSINH(X), the value of the hyperbolic sine of xX 
for all valid input arguments; cosh(a) for DCOSH(X), the value of the 
hyperbolic cosine; and tanh(a) for DTANH(X), the value of the hyperbolic 
tangent. 


USAGE 


A=DSINH (X) 
A=DCOSH (X) 
A=DTANH (X) 


X represents any double-precision value so that the result will be 
presented as a double-precision number. 


GMAP CALLING SEQUENCE 


CALL DSINH(X) for sinh (a) 
CALL DCOSH(X) for cosh (a) 
CALL DTANH(X) for tanh(a) 


METHOD 


Ay DCOSH(X) is calculated as (DEXP(X)+1./DEXP(X))/2.0.° 
2. DSINH(X) is calculated as follows: 


If (/X/.GE.1.44269), DSINH 


It 


(DEXP (X) -1. /DEXP (X))/2.0. 


If (/X/.LT.1.44269), DSINH (DEXPC (X) / (DEXPC (X)+1.0)+DEXPC(X))/2.0. 
3. DTANH(X) is calculated as follows: 

If (X.GE.22.0), DTANH = 1.0. 

If (X.LE.-22.0), DTANH = -1.0. 

If (/X/.LT.1.44269), DTANH = DEXPC(2*X) /(2.+DEXPC(2*X)). 


If (/X/.GT.1.44269.AND./X/.LT.22.0), 
DTANH = (DEXP(X)~-1./DEXP (X) ) / (DEXP(X)+1./DEXP(X)). 


4. These subroutines call other routines that are accurate to double 
precision. The complementary forms of the exponential functions were 
used to help minimize accumulated errors in calculating the 
mathematically correct expressions. 
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EFFECT OF ARGUMENT ERROR 


Lx For DSINH: AER 
RER 


AEA*COSH (X) 
AEA*COSH(X) = REA 


2's For DCOSH: AER 
RER 


AEA*DSINH (X) 
REA*TANH (X) = AEA 


Se For DTANH: AER 
RER 


(1.-TANH (X) **2) *AEA 
2. *AEA/SINH (2*X) 


For a small value of X, RER = REA; for a large value of xX, RER is 
affected less and less by REA. 


RETURNS 


No error handling is performed by this subroutine, but potential 
underflow/overflow conditions will be displayed by the DEXP or DEXPC 
routines that are called by this subroutine. 
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FDSN 


Double Precision Sine and Cosine 


FUNCTION 


FDSN 
for 


computes sin(a) for DSIN(X), the value of the trigonometric sine of xX 
any valid input argument (in radians); and cos(a) for DCOS(X), the 


value of the trigonometric cosine of X for any valid input argument (in 
radians). 


USAGE 


GMAP 


A=DSI 


N(X) 


A=DCOS (X) 


X represents any double-precision number, but accuracy may be affected in 
finding the sine or cosine of large input arguments due to the periodicity 
of the trigonometric functions. 


CALLI 


CALL 
CALL 


METHOD 


NG SEQUENCE 


DSIN(X) for sin(a) 
DCOS(X) for cos (a) 


The method used for this subroutine is equivalent to that used for the 
single-precision sine and cosine subroutines (edit name FSIN), except 
that a more accurate minimax polynomial approximation is used. The 
Minimax polynomial is of the ninth order and has the form 


_ P(Y) =P00+Y* (PO1+Y* (P02+Y* (PO3+Y* (P04+Y* (PO5+Y¥* (P06+Y* (PO7+Y* (P08 
+¥*P09)))))))). 


The minimax polynomial is accurate to more than 19 places. Any error 
in the result is caused only by rounding and truncation; the result 
should be accurate to approximately 17 places. 


EFFECT OF ARGUMENT ERROR 


3/77 


AER = 


NOTE : 


AEA 


The relative error may increase for larger arguments outside of the 
principal range (-PI, +PI). 
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RETURNS 


Normal return is to the next executable statement of the calling program, 
except as specified below: 


1. If (/X/.GT.2**27), FXEM is called with error code 23 to issue a 
warning. Then the sine or cosine is found with as much accuracy as 
possible. | 
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FDSQ FDSQ 


Double Precision Square Root 


FUNCTION 


FDSQ computes the square root of xX in the expression DSQRT(X), where 
X.GE. 0.0. 


USAGE 

A=DSQRT (X) 

X represents any double-precision number .GE.0.0. 
GMAP CALLING SEQUENCE 


CALL DSQRT (xX) 


METHOD 


ls The method used for DSQRT is exactly the same as that used for the 
Ssingle-precision SORT subroutine (edit name FSQR), except that one 
more Newton's iteration is performed (for a total of three). 


EFFECT OF ARGUMENT ERROR 


RER = .5*REA 


RETURNS 


Normal return is to the next executable statement of the calling program, 
except as specified below: 


Le If (X.LT.0.0), FXEM is called with eérror code 22. Then find 
DSQRT (/X/) . 
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Double Precision Tangent 


FUNCTION 


FDTN computes tan(a) for DTAN(X), the value of the trigonometric tangent of 
X for any valid input argument (in radians). 


USAGE 
A=DTAN (X) 


X represents any double~precision number, but accuracy may be affected in 
finding the tangent of large input arguments due to the periodicity of the 
trigonometric functions. 


GMAP CALLING SEQUENCE 


CALL DTAN(X) for tan (a) 


1. The algorithm used to find DTAN(X) is identical to that used for the 
Single-precision tangent subroutine (edit name FTAN), except that a 
more accurate minimax polynomial is used. In this case, the 
polynomials P and Q are of the fourth order; P has the form 

POO+Y* (PO01+Y* (P02+Y* (P03+Y*P04) ) ) 
and Q has the form 
QO0+Y* (QO01+Y* (Q02+Y* (Q03+Y*004))). 
2. The minimax polynomial is accurate to more than 19 places. Any error 


in the result is caused only by rounding and truncation; the result 
should be accurate to approximately 17 places. 


EFFECT OF ARGUMENT ERROR 


AER 
RER 


AEA/ (COS (X) **2) 
2./SIN (2*X) 


NOTE: Near singularities of the expression X=(K+.5)*PI, large errors in 
values may occur. 
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FDTN 


RETURNS 


3/77 


Normal return is to the next executable statement of the calling program, 
except as specified below: 


1. If (/X/.GT.2**27), FXEM is called with error code 90 to issue a 
warning. Then the tangent is found with as much accuracy as 
possible. 
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Double Precision Exponential Complement 


FUNCTION 


FDXC computes e®-1.0 for DEXPC(X), the value of the exponential complement 
of xX (e**X-1.0) for all valid input arguments; 22-1.0 for DEXPC2(X), the 
value of 2.**X-1.0; and 107-1.0 for DXPC10(X), the value of 10.**X-1.0. 


These functions are useful for calculating the above expressions when the 
arguments are near the value zero. However, if the standard exponential 
subroutines were used for calculating an expression such as 2.**X-1.0, 
accuracy would be lost for arguments that are near 0.0 because of the 
subtraction of two almost identical numbers (1.0 and -1.0). This loss of 
accuracy is avoided by using these subroutines. 


USAGE 


A=DEXPC (X) 
A=DEXPC 2 (X) 
=DXPC 10 (X) 


X represents any double-precision value so that the result will be 
presented as a double-precision number. 


GMAP CALLING SEQUENCE 


CALL DEXPC (X) 
CALL DEXPC2 (X) 
CALL DXPC10 (xX) 


‘METHOD 


Lis The basic function calculated by the subroutine is DEXPC2(X). The 
results for DEXPC and DXPC10 are calculated considering that 


DEXPC (X) = DEXPC2(X*1lo0g2 (e) ) and 
DXPC10(X) = DEXPC2(X*log2(10)). 


The constants log2(e) and log2(10) are stored within the subroutine. 


2. For /X/.LT.0.5, the function is calculated using a minimax polynomial 
approximation of the form 


EXPC2 = 2*X*P(X**2)/(Q(X**2) —X*P (X**2) ) 


where: P(Y) 
Q(Y) 


POO+Y* (PO1+Y*P02) and 
QOO+Y* (QO01+Y* (Q02+¥*Q03)). 
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FDXC 


For /X/.GE.0.5, the function DEXPC2 is simply calculated as 
DEXP2 (X)=-1.0. That is, a call to DEXP2 is made followed by a 
subtraction of 1.0. 


The minimax polynomial is accurate to more than 19 places. Any error 
in the result is caused only by rounding and truncation; the result 
should be accurate to approximately 17 places. 


EFFECT OF ARGUMENT ERROR 


RER 
AER 


NOTE : 


RETURNS 


No 


3/77 


AEA 
REA*X 


Serious error amplification can occur if small errors are 
encountered in a large value of xX. 


error handling is performed by this subroutine. However, for large 
positive or negative arguments the subroutine DEXP2 is called and it will 
generate any appropriate error messages. 
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FDXP 


Double Precision Exponentiation, Real Base or Exponent (or Both) 


FUNCTION 


FDXP has five entry points that compute the value of D1l**D2 for any valid 
double-precision numbers Dl and D2. For entry points other than .FDXP2, 


the arguments are converted to double precision. The entry points and 
computation are as follows: 


DPOW(D1,D2) or .FDXP2(D1,D2) for D1**D2 
-FDXP3(D1,R2) for D1**R2 


~FXP4 (R1,D2) 
~FXP5 (I1,R2) 
~FXP6(I1,D2) 


where: DL 
D2 
Il 
I2 
Rl 
R2 


USAGE 


X = D1**D2 


for 
for 
for 


any 
any 
any 
any 
any 
any 


GMAP CALLING SEQUENCE 


CALL DPOW(D1,D2) 
CALL .FDXP2(D1,D2) 
CALL .FDXP3(D1,R2) Notation is described above 


CALL .FXP4(R1,D2) 
CALL .FXP5(1I1,R2) 
CALL .FXP6(I1,D2) 


METHOD 


R1**D2 
I1**R2 
T1**D2 


double-precision number .GE.0.0 
double-precision number 

integer .GE.0 

integer 

single-precision number .GE.0.0 
single-precision number 


l. The subroutine is based on the formula 


X*¥*Y = 2** (Y*DLOG2 (X)). 


26 The subroutine uses the entry point .DLOG2 into the DLOG2 subroutine 


and the entry point .DEXP2 into the DEXP2 subroutine. (DEXP2(Z) finds 
the value of 2**Z.) 


ce The minimax polynomial used is accurate to more than 19 places. Any 
error in the result is caused only by rounding and truncation; the 
result should be accurate to approximately 17 places. 
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EFFECT OF ARGUMENT ERROR 


The error propagation can be calculated from the guidelines provided in the 
DLOG2 and DEXP2 subroutines. 


RETURNS 


Normal return is to the next executable statement in the calling program, 
except as specified below: 


g Ie If (D1.LT.0.0), the result is set to 0.0 and FXEM is called with error 
code 16. 


2% If (D1.EQ.0.0.AND.D2.EQ.0.0), the result is set to 0.0 and FXEM is 
called with error code 17. 


ce If (D1.EQ.0.0.AND.D2.LT.0.0), the result is set to the largest 
positive number and FXEM is called with error code 18. 


4. Appropriate overflow and underflow errors may occur in DEXP2 (edit 


name FDEX) if the results cannot be expressed as a double-precision 
number. 
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_FEXC | FEXC 


Real Exponential Complement 


FUNCTION 


FEXC computes e*-1.0 for EXPC(X), the value of the exponential complement 
of xX (e**X=1.0) for all valid input arguments; 27-1.0 for EXPC2(X), the 
value of 2.**xX-1.0; and 107-1.0 for EXPC10(X), the value of 10.**xX-1.0. 


These functions are useful for calculating the above expressions when the 
arguments are near the value zero. However, if the standard exponential 
subroutines were used for calculating an expression such as 2.**X-1.0, 
accuracy would be lost for arguments that are near 0.0 because of the 
subtraction of two almost identical numbers (1.0 and -1.0). This loss of 
accuracy is avoided by using these subroutines. 


USAGE 


A=EXPC (X) 
A=EXPC2 (X) 
A=EXPC1O (X) 


X represents any value so that the result will be presented as a 
single-precision number. 


GMAP CALLING SEQUENCE 


CALL EXPC (X) 
CALL EXPC2 (X) 
CALL EXPC10 (xX) 


1. The basic function calculated by the subroutine is EXPC2(X). The 
results for EXPC and EXPC10 are calculated considering that 


EXPC(X) = EXPC2 (X*1log2 (e) ) and 
EXPC10(X) = EXPC2(X*log2(10)). 


The constants log2(e) and log2(10) are stored within the subroutine. 


Ze For /X/.LT.1.0, the function is calculated using a minimax polynomial 
approximation of the form 


PO0O+X* (PO1+X* (P02+X* (P03+X* (P04+X* (P05+X* (P06+X*P07)))))). 
Se For /X/.GE.1.0, the function EXPC2 is simply calculated as 


EXP2 (X)-1.0. That is, a call to EXP2 is made followed by a 
subtraction of 1.0. 
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4. The largest known relative error is less than 1*2**-27 (one bit in 
Single precision). 


EFFECT OF ARGUMENT ERROR 


AEA 
REA*X 


RER 
AER 


to 


NOTE: Serious error amplification can occur if small errors are 
encountered in a large value of X. 


' RETURNS 


No error handling is performed by this subroutine. However, for large 
positive or negative arguments, the subroutine EXP2 is called and it will 
generate any appropriate error messages. 
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FEXP FEXP 


Real Exponential, Base 2, and Base 10 


FUNCTION 


FEXP computes e8 for EXP(X), the value of the exponential of X for all 
valid input argument(s); 232 for EXP2(X), the value of 2.0**xX; and 104 for 
EXP10(X), the value of 10**xX, 


USAGE 


A=EXP (X) 
A=EXP 2 (X) 
A=EXP 10 (X) 


X represents any value so that the result will be presented as a 
Single-precision number. 


GMAP CALLING SEQUENCE 


CALL EXP (X) 
CALL EXP2 (X) 
CALL EXP10(X) 


1. The basic function calculated by the subroutine is EXP2(X). The 
results for EXP and EXP10 are calculated considering that 


EXP(X) = EXP2(X*log2(e) ) and 
EXP10(X) = EXP2(X*1log2(10)). 


The constants log2(e) and log2(10) are stored within the subroutine. 
2. Range reduction is performed considering that 
EXP2 (RRX+INTX) = EXP2(RRX) *EXP2 (INTX). 


INTX is the largest integer less than X and RRX is xX minus’ INTX,. 
Therefore, RRX is within the range (0.,1.). 
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3. EXP2(RRX) is found using a minimax polynomial approximation of the 
form 


P00+RRX* (PO01L+RRX* (P02+RRX* (P03+RRX* (P04+RRX* (PO05+RRX* 
(P06+RRX*P07)))))) 


which provides results suitable for single precision. 


4. EXP2 (RRX) *EXP2(INTX) is found by adding INTX to the exponent of the 
result for EXP2(RRX). 


5. The largest known relative error is less than 1*2**-27 (one bit in 
Single precision). 


EFFECT OF ARGUMENT ERROR 


RER 
AER 


AEA 
REA*X 


NOTE: Serious error amplification can occur if small errors are 
encountered in a large value of X. 


RETURNS 
Normal return is to the next executable statement of the calling program, 
except as specified below: 


l. If (result .GT.2.**127), the result is set to the largest number and 
FXEM is called with error code 8. 


2. If (result .LT.2.**-128), the result is set to 0.0 and FXEM is’ called 
with error code 89. 
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FSIN FSIN 


Real Sine and Cosine 


FUNCTION 


FSIN computes sin(a) for SIN(X), the value of the trigonometric sine of X 
for any valid input argument (in radians); and cos(a) for COS(X), the value 
of the trigonometric cosine of X for any valid input argument (in radians). 


USAGE 
A=SIN (X) 
A=COS (X) 
X represents any single-precision number, but accuracy may be affected in 


finding the sine or cosine of large input arguments due to the periodicity 
of the trigonometric functions. 


GMAP CALLING SEQUENCE 


CALL SIN(X) for sin(a) 
CALL COS(X) for cos (a) 


METHOD 


nA The first procedure to be followed in calculating this function is to 
reduce the range of the argument to (-.5,+.5) by employing the 
periodicity of the sine or cosine function. The steps are: 
ae AX = X*(1./PTI) 


b. If the entry was for COS(X), then AX=Ax+.5, since COS(X) = 
SIN (PI/2+X). 


Cx INTAX = the nearest integer to AX. 


d. If INTX is odd, RRX = INTAX minus AX. If INTX is even, RRX = AX 
minus INTAX. 


2. Since RRX is now within the range (-.5,+.5), the polynomial must be 
accurate within the range (-PI/2,+PI/2). 


3. A polynomial approximation is then used to find SIN(RRX) and is of the 
form 


RRX*2*P (4*RRX**2) 


where: P(Y) has the form 


POO+Y* (POL+Y* (P02+Y* (P03+Y* (PO4+Y¥*P05)))). 
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FSIN 


oA 
Zz 


4. The range reduction described in step ld above is where the error for 
large values of the argument may occur. 


5 The largest known relative error is less than 1*2**-27 (one bit in 
single precision). 


EFFECT OF ARGUMENT ERROR 


‘NOTE: “The relative error may increase for large arguments outside of the 
principal range (-PI,+PI). 


RETURNS 
Normal return is to the next executable statement of the calling program, 
except as specified below: 
de If (/X/.GT.2**27), FXEM is called with error code 12 to issue a 


warning. Then the sine or cosine is found with as much accuracy as 
possible. 
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-FSOR : FSOR 


| Real Square Root 


FUNCTION 


FSQR computes the square root of xX in the expression SQRT(X), where 
X.GE. 0.0. 


USAGE 


A=SQRT (X) 


X represents any single-precision number .GE.0.0. 


GMAP CALLING SEQUENCE 


CALL SQRT (X) 


by X is expressed as M*2**E, where 


M = the mantissa of X(.5,1.0) and 
E = the exponent of X. 


2s Range reduction is performed using formulas based on 
E even: SQRT(X)=(2*SQRT(.25*M) ) *2**E/2 
E odd: SORT (X) =(2*SQRT (.5*M) ) *2** (E-1) /2 


3. The SQRT(.25*M) or SQRT(.5*M) is expressed as SQRT(RRX) and E/2 or 
(E-1)/2 is expressed as EXPNT. 


4. SQRT(RRX), where RRX is (.125,.5), is found using an initial 
polynomial approximation accurate to more than two places. Two 
Newton's iterations are then used to find 2*SQRT(RRX) : 


YO = P00+RRX* (P01+P02*RRX) 
Yl = .5*(Y0O+RRX/YO0) 
Y2 = 2*SQRT (RRX) =2*.5 (Y1+RRX/Y0) 


The result is then expressed as 
SQRT(X) = Y2*2**EXPNT 


bs The largest known relative error is less than 1*2**-27 (one bit in 
Single precision). 
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EFFECT OF ARGUMENT ERROR 


RER = .5*REA 


RETURNS 


Normal return is to the next executable statement of the calling program, 
except as specified below: 


Ts If (X.LT.0.0), FXEM is called with error code 13. Then SQRT(/X/) is 
found with as much accuracy as possible, 


3/77 | 4-63 _ _ DD20A 


Real Tangent 


FUNCTION 


FTAN 


FTAN computes tan(a) for TAN(X), the value of the trigonometric tangent of 
X for any valid input argument (in radians). 


USAGE 


A=TAN (X) 


X represents 


finding the tangent of large input arguments due to the periodicity of 
trigonometric functions. 


GMAP CALLING SEQUENCE 


CALL TAN(X) for tan(a) 


3/77 


The first procedure to be followed in calculating this function is 
reduce the range of the argument to (-.5,+.5) by employing 
periodicity of the tangent function. The steps are: 

ae AX=X*(2/PI). The value 2/PI is stored as a constant. 

b. INTAX = the nearest integer to AX. 

Cc. RRX = INTAX minus AX, 


Since RRX is now within the range (-.5,+.5), the polynomial must 
accurate within the range (-PI/4,+PI/4). 


A minimax polynomial is used of the form 
RRX*P (RRX) /Q(RRX) 


where the polynomials P and Q are of the second order; P(Y) has 
form 


PO00+Y* (PO01+Y*P02) 
and Q has the form 
Q00+Y* (Q01+Y¥*Q02). 


If the original argument was in octant l, 2, 5, or 6, the result 
Q(RRX) / (P (RRX) *RRX) . 


any single-precision number, but accuracy may be affected in 


the 


to 


be 


the 


is 
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FTAN FTAN 


If the original argument was in octant 0, 3, 4, or 7, the result is 
-RRX*P (RRX) /Q(RRX). 


These final adjustments are required because of the symmetry of the 
tangent function. 


Da The largest known relative error is less than 1*2**-27 (one bit in 
single precision). 


EFFECT OF ARGUMENT ERROR 


AEA/ (COS (X) **2) 
2./SIN (2*X) 


AER 
RER 


NOTE: Near singularities of the expression X=(K+.5)*PI, large errors in 
values may occur. 


RETURNS © 


Normal return is to the next executable statement of the calling program, 
except as specified below: 


Le If (/X/.GT.2**27), FXEM is called with error code 90 to issue a 
warning. Then the tangent is found with as much accuracy as possible. 
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FTNH 


Real Hyperbolic Sine, Cosine, and Tangent 


FUNCTION 


FTNH computes sinh(a) for SINH(X), the value of the hyperbolic sine of X 
for all valid input arguments; cosh(a) for COSH(X), the value of the 


hyperbolic cosine; and tanh(a) for TANH(X), the value of the hyperbolic 
tangent. 


USAGE 


A=S INH (X) 
A=COSH (X) 
A=TANH (X) 


X represents any value so that the result will be presented as a 
single-precision number. 


GMAP CALLING SEQUENCE 


CALL SINH(X) for sinh(a) 
CALL COSH(X) for cosh(a) 
CALL TANH(X) for tanh(a) 


Ls COSH(X) is calculated as (EXP(X)+1l1./EXP(X))/2.0. 
2% SINH(X) is calculated as follows: 


If (/X/.GE.1.44269), SINH 


(EXP (X) -1./EXP(X))/2.0. 
If (/X/.LT.1.44269), SINH 


(EXPC (X) / (EXPC (X)+1.0) + EXPC(X))/2.0. 
3. TANH (X) is calculated as follows: 


If (X.GE.22.0), TANH = 1.0. 
If (X.LE.-22.0), TANH = ~1.0. 
If (/X/.LT.1.44269), TANH = EXPC(2*X) /(2.+EXPC(2*X)). 


If (/X/.GT.1.44269.AND./X/.LT.22.0), TANH = (EXP(X)-1./EXP(X)) 
/ (EXP (X)+1./EXP(X)). 


4. These subroutines call other routines that are accurate to single 
precision. The complementary forms of the exponential functions were 
used to help minimize accumulated errors in calculating the 
mathematically correct expressions. 
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EFFECT OF ARGUMENT ERROR 


Ls For SINH: AER = AEA*COSH(X) 

RER = AEA*COSH(X) = REA 
2. For COSH: AER = AEA*SINH(X) 

RER = REA*TANH(X) = AEA 


3% For TANH: AER 
RER 


(L.~TANH (X) **2) *ABA 
2. *AEA/SINH ( 2*X) 


For a small value of xX, RER = REA; for a large value of X, RER is 
affected less and less by REA. 


RETURNS 


No error handling is performed by this’ subroutine, but potential 
underflow/overflow conditions will be displayed by the EXP or EXPC routines 
that are called by this subroutine. 
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This section describes library subroutines that are utilized 


SECTION V 


NONMATHEMATICAL LIBRARY SUBROUTINES 


Manner as those described in Section IV. The 
section are nonmathematical in nature. 


subroutines described 


Table 5-1 lists the nonmathematical library subroutines. 


FMDE 


PCAT 


FSRT 


FTSF 


FTSG 


FRSW 


FSHF 


FRND 


FFER 


Table 5-1. Nonmathematical Library Subroutine 


Function 


Mode determination by FORTRAN 
Character string manipulation 
Sort character array of data 
Memory usage 

User time and identification 
Set/reset switch word 
Shift/rotate word contents 


Random number generator 


File and Record Control I/O error 
recovery . 


CD600E1.075 
CbD600E1.069 
CD600E1.070 
CD600E1.072 
CD600E1.071 
CD600E1.074 
CD600E1.073 
CD600D3.009 


CD600E1.011 


in the 


in 


same 
this 
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Mode Determination by FORTRAN 


' FUNCTION 


FMDE is used to provide the FORTRAN user with means to determine whether 
the environment is batch or time sharing and whether it is BCD or ASCII. 


CALLING SEQUENCE: 


j = MODE(i) 


for i 1, j = 0 if batch, j = 1 if time sharing 


for i= 2, j 


Mt 


0 if BCD, j = 1 if ASCII 


If iis neither 1 nor 2, j is set to -l. 


METHOD 


The FORTRAN compiler and library provide flags for determination of mode 
and character set. These are interrogated to give information to the user. 
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“Character String Manipulation 


FUNCTION 


FCAT is used to provide the FORTRAN user with ability to move or compare abs 
character substring of arbitrary length and position within a string. 


CALLING SEQUENCE: 


CALL CONCAT (a,n,b,m,f) 


where: 
a = string to be replaced 
n = initial character of a (n=l implies first character) 
b = replacement string 
m = initial character of b (m=1 implies first character) 
f = number of characters to be replaced; if f£ is not given, 1 is 
assumed 
Causes characters n through (n + f -1) of a to be replaced with 
characters m through (m + f -1) of b. 
Example: 
0010 CHARACTER A*20/"FIFTEEN WERE THERE "/ 
0020 CHARACTER B*20/"SIXTEEN WERE ABSENT "/ 
0030 PRINT A,B 
0040 CALL CONCAT (A,1,B,1,3) 
0050 PRINT,A,B 
0060 STOP;END 
READY 
* RUN 


FIFTEEN WERE THERE SIXTEEN WERE ABSENT 
SIXTEEN WERE THERE SIXTEEN WERE ABSENT 


i = KOMPCH (a,n,b,m,f) is function where substring of b is 
compared character by character to substring of a. 


if b 
at: 5 


ai b 


METHOD 


=a, I= 0 
>a, l=l 


<a, Il=-e-1 


The move is done by the way of the SC modification after tally words are built 
using parameters in the calling sequence. 
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Sort Array of Data 


FUNCTION 


FSRT is used to sort an array of data in ascending or descending order. 
SORT is for ascending sort; SORTD is for descending sort. 


CALLING SEQUENCE: 


CALL SORT/SORTD (array, nrec, lrs, key, ees key) 


array = the name of the array to be sorted. 

nrec = the number of items, or logical records in the array. 

lrs = the logical record size, or the size in words of each item in the 
array. 

key = the relative word number of the ith sort key in each logical record 
and must be such that 0< key, <lrs. Record comparisons are made 
starting with key and either’ progress through to key,, or until a 
non-equal comparison is made. Any number of sort keys may be 
specified; however, at least one must always be specified. 

METHOD 


Implementation is by use of the Shell method as described in ACM Computer 
Survey, December 1971. 
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Memory Usage 


FUNCTION 


FTSF is used to obtain memory usage of a program, the terminal © station 
code, and to move data from/to a 10-word memory file. 


CALLING SEQUENCE: 


CALL MEMSIZ (3) 
j is returned as the number of 1024-word blocks of memory currently . 
allocated this job. The return is a 3-character variable. MEMSIZ can 
also be used as a function; e.g., IF (MEMSIZ (j) .GE.n) ... 


CALL TERMNO (a) 


where a is a type character variable capable of receiving the 
2-character station code. 


For batch, the call returns blanks. 


CALL CORFL (loc,i,j,k) 


loc = first word address of area to or from which data is to be moved 
i = number of words to be moved such that 1< i <10 
= relative location in the 10-word area at which transfer is to begin 
k = 0, data to be transferred into the 10-word area 
k = 1, data to be transferred out of the 10-word area 


This call is ignored in batch. 


METHOD 


Memory size is determined from the BAR. TERMNO and CORFL use_ services 
provided by the time sharing system. 
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User Time 


. FUNCTION 


FTSG 


And Identification 


is used to obtain user time and user identification. 


CALLING SEQUENCE 


CALL 


CALL 


CALL 


CALL 


METHOD 


CORSEC (a) 

a, a real variable, is returned as the product of 1024-word blocks 
currently allocated and processor time in seconds. CORSEC can also be 
used as a function. For example, 


IF (CORSEC (A) .GE.B) 


TERMTM (a) 


a, a real variable, is returnec as hours since log-on time; ignored in 
batch. 


USRCOD (s) 


Ss is character type variable of at least 12 characters into which the 
user identification will be returned. (Time sharing only; ignored in 
batch. ) 


PTIME (a) 


a, a real variable, is returned as elapsed processor time in hours. 
PTIME can also be used as a function. 


Service routines provided by the time sharing system are used to obtain 


data. 
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Set/Reset Switch Word 


FUNCTION 


FRSW sets or resets bits in the program switch word. 


‘CALLING SEQUENCE 


j = ISETSW(i) for setting switch word. 
j = IRETSW(i) for resetting switch word. 


i contains an octal value that specifies one or more bits that are to 
be set or reset. 


HW 


1 set bit on 
0 = leave bit as is 


j is resultant switch word setting in the respective positions of the 
program switch word. 


J 
0 


bit is on 
bit is off 


In time sharing, bits 0-17 cannot be changed. 


METHOD 


MME GESETS/GERETS or DRL SETSWH/RSTSWH is used. 
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Shift/Rotate Word Contents 


FUNCTION 


FSHF shifts/rotates word contents. 


CALLING SEQUENCE: 


j = ILS (k,p) k left shift by p bit positions 

j = IRS(k,p) k right shift by p bit scateaons (fill vacated positions 
with contents of bit position zero) 

j = ILR(k,p) k left rotate by p bit positions 

j = ILR(k,p) k right logical shift by p bit positions (fill vacated 
positions with zeros) 
j = new word 
k = word to shift 
Pp = number of bits 

METHOD 


k is shifted as indicated by p bit positions and stored into j. k and p are 
assumed to be integers; if j is not an integer, the compiler provides a 
conversion prior to storage. 
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- Random Number Generator 


FUNCTION 


FRND provides four different calls to generate random numbers from a 
uniform (rectangular) distribution. 


USAGE 


A = RAND (range) (seed is always 1) 

Use of this version of the function results in the same sequence of random 
numbers each time the program is executed. 

A = RANDT (range) (seed = current time register value) 


Use of this version results in a different sequence of random numbers’ each 
time the program is executed. 


For the above two uses, the algorithm used is the same. A is a _ positive 
real number such that 0O<A <range. 


A = FLAT (seed) 

This version allows the seed to be varied, but the range is constant 
(0 A l) 

A = UNIFM2 (seed, mean, width) 


This version allows the seed and the range to be varied. The range will be 
such that (mean - width/2 <A <mean + width/2). 


For example: 

A = UNIFM2 (6.5, O65: 150) 

will generate a set of random numbers between 0.0 and 1.0. 
For the above two uses, the algorithm used is the same. 


NOTE: For all calls, the value of the initial argument (seed) initializes 


the algorithm for generation of the sequence of random numbers. For 
all subsequent calis to the function, during execution of the same 
program unit, the value of the argument is ignored. 
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File and Record Control I/O Error Recovery 


FUNCTION 


FFER provides the user with some control on File and Record Control erro s 
Only (not FORTRAN). The subroutine allows the user to set his error routi se 
address into the file control block; a default address is provided if the 
user does not provide the address. 


File control block (FCB) word -5 contains an address for user. error 
recovery. When a file is first used, this word (FCB -5) is set with a 
pointer to FOPE and the user-specified address, if one is supplied, is 
placed in the lower half of FCB -15. At any time, the user can place his 
recovery address in FCB -15 by the following call. 


CALLING SEQUENCE 


CALL FLGRC (lgu,ptr) 


lgu = numeric file code 


ptr address of the recovery routine 


METHOD 


See Function. 
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INDEX 


.FBDT. 
~FBDT. 


sF BLT. 
-FBLT. 


~FBST. 
CALL .FBST. 


~FCFDP 
CALL .FCFDP 


-FCFMP 
CALL .FCFMP 


~FCNV. 
~FCNV. 


~FCNVC. 
-FCNVC. 


-FCNVD. 
-FCNVD. 


-FCNVI. 
«PCNVI: 


-FCNVL. 
e FCNVL e 


-FCNVR. 
-FCNVR. 


~FCXP1 
CALL .FCXP1 


oFCXP2 
CALL .FCXP2 


oFCXP3 
CALL .FCXP3 


»FDXP1 
CALL .FDXP1 


 .FDXP2 
CALL .FDXP2 


»-FDXP3 
CALL .FDXP3 


~FEFT. 
CALL .FEFT. 
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° FEOF e 


CALL .FEOF. 2-36 
.-FEXIT 

CALL .FEXIT 2-50 
» FGTFB 

CALL .FGTFB 2-37 
.FLTPR | 

CALL .FLTPR 2-31 
~FOPEN 

CALL .FOPEN 2-37 
-  FRWT. 

CALL .FRWT. 2-34 
-FSLEW 

CALL .FSLEW 2-39 
e FXEM ° 

CALL .FXEM. 3-14 
~FXP1 

CALL .FXP1 4-12 
. FXP2 | 

CALL .FXP2 4-14 
~FXP3 

CALL .FXP3 4-16 
~FXP4 

CALL .FXP4 4-21 
~FXP5 

CALL .FXP5 4-21 
- FXP6 

CALL .FXP6 4-21 
»FXSW1 

»FXSW1 3-13 
«FXSW2 

~FXSW2 3-14 
-FXSW3 

~-FXSW3 3-14 
« RCOV 

CALL .RCOV 3-6 
-SETU. 

-SETU. 2-28 
ABSOLUTE 

Complex Absolute Value . 4-10 
ACCESS 

Access a Permanent File 3-28 
AFT 

Close File, Detach Buffers, Remove from AFT 3-30 
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ALOG 


CALL ALOG 4-39 
ALOG10 

CALL ALOG10 4-39 
ARCOS 

CALL ARCOS 4-60 
ARCTANGENT 

Double Precision Arctangent 4-31 

Real Arctangent 4-42 
ARITHMETIC 

Arithmetic Fault Processor 2-31 
ARRAY 

Sort Array of Data 5-4 
ARSIN 

CALL ARSIN 4-60 
ASCII/BCD 

ASCII/BCD Indicators 3-7 
ASSIGNMENT | 

Character String Assignment 2-23 
ATAN 

CALL ATAN 4-42 
ATAN2 

CALL ATAN2 4-42 
ATTACH 

Attach a Temporary Mass Storage or Terminal File 3-31 

CALL ATTACH 3-28 
BACKS PACE 

Backspace Record 2-32 

File Forwardspace and Backspace 2-51 
BINARY 

Linked Binary Input/Output Interface 2-4 

Random Binary Input/Output 2-10 

Specify Record Size, Random Binary File 2-54 
BLOCK 

File Control Block 2-44 
BUFFERS 

Close File, Detach Buffers, Remove from AFT 3-30 
CABS . 

CALL CABS 4-10 
CALL 

CALL .FBST. 2-32 

CALL .FCFDP 4-7 

CALL .~FCFMP 4-7 

CALL .FCXP1 4-19 

CALL .FCXP2 4-18 

CALL .FCXP3 4-18 

CALL .FDXP1 4-14 

CALL .FDXP2. 4-21 

CALL .FDXP3 4-21 
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CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
‘CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALE 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CAEL 
CALL 
CARL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 


~FEFT. 


FEOF. 
~PEXIT 
.FGTFB 
~FLTPR 
- FOPEN 
-FRWT. 
-F SLEW 
~FXEM. 
~FXPL 
~FXP2 
~FXP3 
~FXP4 
~FXP5 
~FXP6 
-RCOV 
ALOG 
ALOG10 
ARCOS 
ARSIN 
ATAN 
ATAN2 
ATTACH 
CABS 
CALLSS 
CCOS 
CEXP 
CLOG 
CONCAT 
CORFL 
CORSEC 
COS 
CREATE 
CSIN 
CSORT 
DATAN 
DATAN2 
DATIM 
DCOS 
DEFIL 
DETACH 
DEXP 
DLOG 
DLOG10 
DMOD 
DSIN 
DSQRT 
DUMP. 
DUMP 
DUMPA 
DUMPA 
DVCHK 
EXIT 
EXP 
FCLOSE 
FILBSP 
FILFSP 
FLGEOF 
FLGERR 
FLGFRC 
FLGRC 
FMEDIA 
FPARAM 
FXALT 
FXEM 
FXOPT 
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CALL LINK 
CALL LLINK 
CALL MEMSIZ 
CALL OVERFL 
CALL PDUMP 
CALL PDUMP 
CALL PDUMPA 
CALL PDUMPA 
CALL PTIME 
CALL RANSIZ 
CALL RANSIZ 
CALL SETBUF 
CALL SETFCB 
CALL SETLGT 
CALL SIN 
CALL SLITE 
CALL SLITET 
CALL SORT/SORTD 
CALL SQRT 
CALL SSWTCH 
CALL TANH 
CALL TERMNO 
CALL TERMTM 
CALL USRCOD 
Call TSS Subsystem 


CALLGT 
CALLGT 


CALLSS 
CALL CALLSS 


CARRIAGE 
Carriage Control Simulator 


CCOS 
CALL CCOS 


CEXP 
CALL CEXP 


CHARACTER 
Character String Assignment 
Character String Compare 
Character String Manipulation 


CLOG 
CALL CLOG 


CLOSE 
Close File, Detach Buffers, Remove from AFT 


COMMON 


Double Precision Natural and Common Logarithms 
Natural and Common Logarithms 


COMPARE 
Character String Compare 
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COMPLEX : 
~ Complex Absolute Value 
Complex Base and Double Precision Exponent 
Complex Base and Exponent 
Complex Base and Real Exponent 
Complex Exponential 
Complex Multiplication and Division 
Complex Natural Logarithm 
Complex Sine and Cosine 
Complex Square Root 
Exponentiation - Complex Base, Integer Exponent 


CONCAT 
CALL CONCAT 


CONSOLE 
Console Communication 


CORFL 
CALL CORFL 


CORSEC 
CALL CORSEC 


COS 
CALL COS 


COSINE 
Complex Sine and Cosine 
Double Precision Sine and Cosine 
Real Sine and Cosine 


CREATE 
CALL CREATE 
Create TSS Temporary File 


CSIN 
CALL CSIN 


CSORT 
CALL CSQRT 


DATA~ERROR 
Initialization of Data-Error Processing 


DATAN 
CALL DATAN 


DATAN2 
CALL DATAN2 


DATE 
Date and Time 


DATIM 
CALL DATIM 


DCOS 
CALL DCOS 


DEBUG 
Object Time Debug Processor 


DEFIL 
CALL DEFIL 
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DETACH 
CALL DETACH 
Close File, Detach Buffers, Remove from AFT 


DEXP 
CALL DEXP 


DIVIDE 
Exponent Register Overflow and Divide Check Tests 


DIVISION 
Complex Multiplication and Division 


DLOG 
CALL DLOG 


DLOG10 
CALL DLOG10 


DMOD 
CALL DMOD 


DOUBLE 
Complex Base and Double Precision Exponent 
Double Precision Arctangent 
Double Precision Base and Exponent 
Double Precision Base and Real Exponent 
Double Precision Exponential 
Double Precision Modulus 
Double Precision Natural and Common Logarithms 
Double Precision Powers of Ten Table 
Double Precision Sine and Cosine 
Double Precision Square Root 
Integer Base and Real or Double Precision Exponent 
Real Base and Double Precision Exponent 


DSIN 
CALL DSIN 


DSQRT 
CALL DSQRT 


DUMP 
CALL DUMP 
CALL DUMP 
Memory Dump 
Namelist and Dump Output 


DUMPA 
CALL DUMPA 
CALL DUMPA 


DVCHK 
CALL DVCHK 


END-OF-FILE 
End-of-File (On Input) Processor 
Initialization of End-of-File Processing 


ENDFILE 
Rewind and Endfile Processor 


ERROR 
Execution Error Monitor 
I/O Error Recovery 
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EXIT 


CALL EXIT 


EXP 


CALL EXP 


EXPONENT 


Complex Base and 
Complex Base and 
Complex Base and 
Double Precision 
Double Precision 
Exponent Register Overflow and Divide Check Tests 


Double Precision Exponent 
Exponent 

Real Exponent 

Base and Exponent 

Base and Real Exponent 


Exponentiation - Complex Base, Integer Exponent 
Exponentiation - Floating-Point Base Integer Exponent 
Exponentiation - Integer Base and Exponent 
Exponentiation - Real Base and Exponent 


Integer Base and 


Real or Double Precision Exponent 


Real Base and Double Precision Exponent 


EXPONENTIAL 


Complex Exponential 
Double Precision Exponential 
Natural Exponential 


EXPONENTIATION 


Exponentiation - Complex Base, Integer Exponent 
Exponentiation Floating-Point Base Integer. Exponent 
Exponentiation - Integer Base and Exponent 


Exponentiation 


F1XP 
F1XP 


F2XP 
F2XP 


F3XP 
F3XP 


F4XP 
F4XP 


FALG 
FALG 


FASC 


FASC 


‘PASN 
FASN 


FATN 
FATN 


FAULT 


Real Base and Exponent 


Arithmetic Fault Processor 


FBCD 
FBCD 


FBLO 
FBLO 
FBLO 
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BST 
FBST 


FCAB 
FCAB 


FCAL 
FCAL 


{CAT 
FCAT 


FCEX 
FCEX 


FCHA 
FCHA 


FCLG 
FCLG 


FCLO 
FCLO 


FCLOSE 
CALL 


FCMP 
FCMP 


FCOM 
FCOM 


FCRA 
FCRA 


FCSL 
FCSL 


FCSN 
FCSN 

FPCSO 
FCSQ 


FCXP 
FCXP 


FDAT 
FDAT 


FDBG 
FDBG 


FDEF 
FDEF 

FDEX 
FDEX 


FDIO 
FDIO 


FDLG 
FDLG 


FDMD 
FDMD 


FCLOSE 
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' FDMP 
FDMP 


FDPT 
FDPT 


FDSN 
, FDSN 


FDSQ 
FDSQ 


FDTH 
FDTH 


-FDTM 
FDTM 


FDXP 
FDXP 


FEFT 
FEFT 


FEOF 
FEOF 


FEXP 
FEXP 


FFEE 
FFEE 


FFER 
FFER 
FFER 


FFFB 
FFFB 


FFLT 
FFLT 


“FIFA 
FIFA 


FILBSP 
CALL FILBSP 


FILE 
Access a Permanent File 


Attach a Temporary Mass Storage or Terminal File 


Close File, Detach Buffers, Remove from AFT 


Create TSS Temporary File 


File Closing 
File Control Block 


File Forwardspace and Backspace 


File Opening 
File Transliteration 


Specify Record Size, Random Binary File 


FILFSP 
CALL FILFSP 


FLAT 
FLAT 
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FLGEOF 
CALL 


FLGERR 
CALL 


FLGFRC 
CALL 


FLGRC 
CALL 


FLHS 
FLHS 


FLIT 
FLIT 


FLNK 
FLNK 


FLGEOF 


FLGERR > 


FLGFRC 


FLGRC 


FLOATING~POINT 


Exponentiation - Floating-Point Base Integer Exponent 


FMDE 
FMDE 


FMED 
FMED 


FMEDIA 
CALL 


FNLI 
FNLI 


FNLO 
FNLO 


FOFA 
FOFA 


FOPE 
FOPE 


FORMAT 


Format Controlled Sequential Input/Output 


FORTRAN 


Mode Determination by FORTRAN 


FMEDIA 


FORWARDSPACE 


File Forwardspace and Backspace 


FPARAM 
CALL 


'| FPAW 


FPAW 


FRBC 
FRBC 


FRCV 
FRCV 


FRDA 
FRDA 


FPARAM 
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FRDB 
_ . FRDB 


FRDD 
FRDD 


“FRND 
FRND 


FRRD 
FRRD 
FRRD 
FRRD 


FRSW 
FRSW 


FSET 
FSET 


FSHF 
FSHF 


FSIN 
FSIN 


FSLI 
FSLI 
FSLI 


FSLO 
FSLO 
FSLO 


FSLW 
FSLW 


FSOR 
FSQR 


FSRT 
FSRT 


FSTU 


FSTU - 


FSTU 


FSWI 
FSWI 


FTAC 
FTAC 


FTGF 
FTGF 


FTLK 
FTLK 


FTNH 
FTNH 


FTSF 
FTSF 


FTSG 
FTSG 
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‘FTSU 
FTSU 


FVFI 
FVFI 


FVFO 
FVFO 


FXALT 
CALL 


FXEM 
CALL 
FXEM 


FXER 
FXER 


FXIT 
PAIT 


FXMA 
FXMA 


FXOPT 
CALL 


FXALT 


FXEM 


FXOPT 


GENERATOR 
Random Number Generator 


H* 

Restore Link - H* 
HYPERBOLIC 

Real Hyperbolic Tangent 


1/0 


I/O Error Recovery 
I/O Parameters of Run-Time Library 


IDENTIFICATION 


User 


ILR 
ILR 


ILS 
ILS 


Time And Identification 


INDICATORS 
ASCII/BCD Indicators 


INITIALIZER 
Pre-execution Initializer 
Pre-execution Initializer 


INPUT 


End-of-File (On Input) Processor 
Namelist Input 
Terminal Input Recovery 


DD20 


TNEULCOUTEUL 
Format Controlled Sequential Input/Output 
INPUT/OUTPUT LIBRARY 
INPUT/OUTPUT LIBRARY SUBROUTINES 
Linked Binary Input/Output Interface 
Random Binary Input/Output 
Short List Input/Output Processor 


INTEGER 
Exponentiation - Complex Base, Integer Exponent 


Exponentiation - Floating-Point Base Integer Exponent 


Exponentiation - Integer Base and Exponent 


Integer Base and Real or Double Precision Exponent 


INVERSE 
Inverse Sine/Cosine 


IRETSW 
_ IRETSW 


‘IRS 
IRS 


ISETSW 
ISETSW 


LIBRARY 
I/O Parameters of Run-Time Library 
INPUT/OUTPUT LIBRARY 
INPUT/OUTPUT LIBRARY SUBROUTINES 
LIBRARY CALLS 
MATHEMATICAL LIBRARY 


LINK 
CALL LINK 
Restore Link - H* 
Linked Binary Input/Output Interface 
Restore Links During Execution (Batch) 
Restore Links During Execution (Time Sharing) 


LIST 
Short List Input/Output Processor 


LLINK 
CALL LLINK 


LOGARITHM 
Complex Natural Logarithm 
Double Precision Natural and Common Logarithms 
Natural and Common Logarithms 


“LOGICAL 
DEFINE LOGICAL UNIT TABLE 
Logical Unit Table Routines 


MANIPULATION 
Character String Manipulation 


MASS 
Attach a Temporary Mass Storage or Terminal File 


MATHEMATICAL 
MATHEMATICAL LIBRARY 
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MEMORY 


Memory Dump 2-43 

Memory Usage 5-5 
MEMSTZ 

7 CALL MEMSIZ 5-5 

MODE : 

MODE 5-2 

Mode Determination by FORTRAN 5 (9 
MODULUS 

Double Precision Modulus 4=5 
MONITOR 

Execution Error Monitor 379 
MULTIPLICATION 

Complex Multiplication and Division 4-7 
NAMELIST 

Namelist and Dump Output 2-21 

Namelist Input 2-17 
NATURAL 

Complex Natural Logarithm 4-54 

Double Precision Natural and Common Logarithms 4-34 

Natural and Common Logarithms 4-39 

Natural Exponential 4-37 
NUMBER 

Random Number Generator 5-9 
OBJECT 

Object Time Debug Processor 2-27 
OUTPUT 

Namelist and Dump Output 2-21 

Output Stop and Pause Information 2-25 
OVERFL 

CALL OVERFL 3-8 
OVERFLOW 

Exponent Register Overflow and Divide Check Tests 3-8 
PARAMETERS 

I/O Parameters of Run-Time Library 2-55 
PAUSE 

Output Stop and Pause Information 2-25 
PDUMP 

CALL PDUMP 2-43 

CALL PDUMP 2-21 
PDUMPA 

CALL PDUMPA 2-21 

CALL PDUMPA 2-43 
PERMANENT | 

Access a Permanent File 3-28 
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PRE-~EXECUTION 
Pre-execution Initializer 
Pre-execution Initializer 


PRECISION 
Complex Base and Double Precision Exponent 
Double Precision Arctangent 
Double Precision Base and Exponent 
Double Precision Base and Real Exponent 
Double Precision Exponential 
Double Precision Modulus 
Double Precision Natural and Common Logarithms 
Double Precision Powers of Ten Table 
Double Precision Sine and Cosine 
Double Precision Square. Root 
Integer Base and Real or Double Precision Exponent 
Real Base and Double Precision Exponent 


PROCESSING 
Initialization of Data-Error Processing 
Initialization of End-of-File Processing 


PTIME 
CALL PTIME 


RAND 
RAND 


RANDOM 
Random Binary Input/Output 
Random Number Generator 
Specify Record Size, Random Binary File 


RANDOM RECORD SIZE 
Random Record Size 


RANDT 
RANDT 


RANSIZ 
CALL RANSTZ 
CALL RANSI2Z 


REAL 

Complex Base and Real Exponent 

Double Precision Base and Real Exponent 
Exponentiation - Real Base and Exponent 


Integer Base and Real or Double Precision Exponent 
Real Arctangent 

Real Base and Double Precision Exponent 

Real Hyperbolic Tangent 

Real Sine and Cosine 

Real Square Root 


RECORD 
Backspace Record 
Specify Record Size, Random Binary File 


RECOVERY 
_I/O Error Recovery 
Terminal Input Recovery 


REGISTER 
Exponent Register Overflow and Divide Check Tests 


REMOVE 
Close File, Detach Buffers, Remove from AFT 
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RESTORE 
Restore Link - H* 
Restore Links During Execution (Batch) 


Restore Links During Execution (Time Sharing) 


REWIND 
Rewind and Endfile Processor 


RUN-TIME 
I/O Parameters of Run-Time Library 


SENSE 
Sense Light Simulator 
Sense Switch Test 


SEQUENTIAL 
Format Controlled Sequential Input/Output 


SET/RESET 
Set/Reset Switch Word 


_SETBUF 
CALL SETBUF 


SETFCB 
CALL SETFCB 


SETLGT 
CALL SETLGT 


SHIFT/ROTATE 
Shift/Rotate Word Contents 


SIN 
CALL SIN 


SINE 
Complex Sine and Cosine 
Double Precision Sine and Cosine 
Real Sine and Cosine 


SINE/COSINE 
Inverse Sine/Cosine 


SIZE 
Specify Record Size, Random Binary File 


SLITE 
CALL SLITE 


SLITET 
CALL SLITET 


SORT 
Sort Array of Data 


SORT/SORTD 
CALL SORT/SORTD 


SPECIFY 
Specify Record Size, Random Binary File 


SQRT 
CALL SQRT 
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SQUARE ROOT 
Complex Square Root 
Double Precision Square Root 
Real Square Root 


SSWTCH 
CALL SSWTCH 


STOP 
Output Stop and Pause Information 


STORAGE 
Attach a Temporary Mass Storage or Terminal File 


STRING 

Character String Assignment 
Character String Compare 
Character String Manipulation 


SUBSYSTEM 
Call TSS Subsystem 


TANGENT 
Real Hyperbolic Tangent 


TANH 
CALL TANH 


TEMPORARY 
Attach a Temporary Mass Storage or Terminal File 
Create TSS Temporary File 


TERMINAL 
Attach a Temporary Mass Storage or Terminal File 
Terminal Input Recovery 


TERMINATION 
Job Termination 


TERMNO 
CALL TERMNO 


TERMTM 
CALL TERMTM 


TRANSLITERATION 
File Transliteration 


TSS 
Call TSS Subsystem 
Create TSS Temporary File 


UNIFM2 
UNIFM2 


USER 
User Time And Identification 


USRCOD_ 
CALL USRCOD 


VALUE 
Complex Absolute Value 
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